Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_001416 41,601 (AAC)2→1 coding (521‑523/873 nt) NinC → NinC protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_001416 41601 41603 3 29 [0] [0] 28 NinC NinC protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_001416 = 416000 (0.000)23 (1.000) 19/90 0.0 100% coding (520/873 nt) NinC NinC protein
?NC_001416 41604 = 0 (0.000)coding (524/873 nt) NinC NinC protein

TGAACTGTCAGACTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACAACCATTCGATTTGCAAATACCGGAACATCTCGGTAACTGCATATTCTG  >  NC_001416/41551‑41649
                                                  |                                                 
tGAACTGTCAGACTTTGAGAAGGAAGATATTCTCGCATGGTGGAAGCAAC                                                   >  1:727‑M1/1‑50 (MQ=255)
    cTGTCAGACTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACcatt                                               >  1:2919‑M1/1‑46 (MQ=255)
      gTCAGACTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACcattcg                                             >  1:2683‑M1/1‑44 (MQ=255)
        cAGACTTTGGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACcattcga                                            <  1:2944‑M1/50‑8 (MQ=255)
           aCTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACcattcgatttg                                        >  1:1852‑M1/1‑39 (MQ=255)
            cTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACcattcgatttgc                                       <  1:1426‑M1/50‑13 (MQ=255)
            cTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACcattcgatttgc                                       >  1:3385‑M1/1‑38 (MQ=255)
                     aggaagATATCCTCGCATGGTGGAAGCAACcattcgatttgcaaataccg                               <  1:4972‑M1/50‑21 (MQ=255)
                      ggtagataTCCTCGCATGGTGGAAGCAACcattcgatttgcaaataccgg                              <  1:3597‑M1/47‑22 (MQ=255)
                         agATATCCTCGCATGGTGGAAGCAACcattcgatttgcaaataccggaac                           <  1:2288‑M1/50‑25 (MQ=255)
                          gATACCCTCGCATGGTGGAAGCAACcattcgatttgcaaataccggaaca                          <  1:4834‑M1/50‑26 (MQ=255)
                           atatCCTCGCATGGTGGAAGCAACcattcgatttgcaaatactggaacat                         >  1:1043‑M1/1‑24 (MQ=255)
                           atatCCTCGCATGTGGAAGCAACcattcgatttgcaaatacggaacatct                        >  1:604‑M1/1‑23 (MQ=255)
                             atCCTCGCATGGTGGAAGCAACcattcgatttgcaaataccggaacatct                       >  1:3580‑M1/1‑22 (MQ=255)
                                cTCGCATGGTGGAAGCAACcattcgatttgcaaataccggaacatctcgg                    >  1:3982‑M1/1‑19 (MQ=255)
                                  cGCATGGTGGAAGCAACcattcgatttgcaaataccggaacatctcggta                  >  1:1076‑M1/1‑17 (MQ=255)
                                   gCATGGTGGAAGCAACcattcgatttgcaaataccggaacatctcggtaa                 >  1:4349‑M1/1‑16 (MQ=255)
                                   gCATGGTGGAAGCAACcattcgatttgcaaataccggaaaatctcggtaa                 >  1:419‑M1/1‑16 (MQ=255)
                                     aTGGTGGAAGCAACcattcgatttgcaaataccggaacatctcggtaact               <  1:648‑M1/50‑37 (MQ=255)
                                      tggtggCAGCAACcattcgatttgcaaataccggaacatctcggtaactg              <  1:2513‑M1/50‑38 (MQ=255)
                                      tggtggAAGCAACcattcgatttgcaaataccggaacatctcggtaactg              <  1:3874‑M1/50‑38 (MQ=255)
                                           gAAGCAACcattcgatttgcaaataccggaacatctcggtaactgcatat         >  1:1257‑M1/1‑8 (MQ=255)
                                            aaGCAACcattcgatttgcaaataccggaacatctcggtaactgcgtatt        >  1:4347‑M1/1‑7 (MQ=255)
                                             aGCAACcattcgatttgcaaataccggaacatctcggtaactgcatattc       >  1:901‑M1/1‑6 (MQ=255)
                                               cAACcattcgatttgcaaataccggaacgtctcggtaactgcatattctg     <  1:2162‑M1/50‑47 (MQ=255)
                                               cAACcattcgatttgcaaataccggaaacatctcggtaactgcatattct     >  1:1668‑M1/1‑4 (MQ=255)
                                                 aCcattcgatttgcaaataccggaacatctcggtaactgcatattctgca   >  1:3156‑M1/1‑2 (MQ=255)
                                                  ccattcgatttgcaaataccggaacatctcggtaactgcatattctgcat  <  1:4063‑M1/50‑50 (MQ=255)
                                                  ccattcgatttgcaaataccggaacatctcggtaactgcatattctgcat  <  1:849‑M1/50‑50 (MQ=255)
                                                  |                                                 
TGAACTGTCAGACTTTGAGAAGGAAGATATCCTCGCATGGTGGAAGCAACAACCATTCGATTTGCAAATACCGGAACATCTCGGTAACTGCATATTCTG  >  NC_001416/41551‑41649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.