Predicted mutation
evidence seq id position mutation annotation gene description
MC JC AM260479 1,464,441 Δ37 bp coding (125‑161/552 nt) h16_A1354 → conserved hypothetical protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ AM260479 1464441 1464477 37 23 [1] [0] 22 h16_A1354 conserved hypothetical protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260479 = 14644401 (0.040)22 (0.810) 18/212 0.4 97.8% coding (124/552 nt) h16_A1354 conserved hypothetical protein
?AM260479 1464478 = 0 (0.000)coding (162/552 nt) h16_A1354 conserved hypothetical protein

GATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC  >  AM260479/1464311‑1464567
                                                                                                                                 |                                                                                                                               
gaTGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctg                                                                                                                            <  2:285393‑M1/134‑6 (MQ=255)
    ccATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacg                                                                                                                        >  1:512028‑M1/1‑125 (MQ=255)
    ccATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacg                                                                                                                        >  1:525485‑M1/1‑125 (MQ=255)
                 tgACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcc                                                                                                            <  1:256727‑M1/133‑22 (MQ=255)
                     cTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                                       <  2:388757‑M1/134‑27 (MQ=255)
                                   ggtggGCGAACCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacgg                                                                                         >  1:506153‑M1/1‑94 (MQ=255)
                                           aaCCGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggt                                                                                 <  1:477985‑M1/134‑49 (MQ=255)
                                              cGTCCAGCGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                              <  1:357321‑M1/134‑52 (MQ=255)
                                                          ccgccACACGAGCACCGGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaata                                                                    <  1:546828‑M1/133‑62 (MQ=255)
                                                             ccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgg                                                                                                                       >  1:130006‑M1/1‑69 (MQ=255)
                                                             ccACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgg                                                                                                                       <  2:130006‑M1/79‑11 (MQ=255)
                                                                               gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCCAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                              >  2:97506‑M1/1‑51 (MQ=255)
                                                                               gCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgag                                                                              <  1:97506‑M1/102‑52 (MQ=255)
                                                                                    tGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggt                                           <  2:512028‑M1/132‑87 (MQ=255)
                                                                                     ggCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcga                                        >  2:59726‑M1/1‑45 (MQ=255)
                                                                                       cATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggttagccgccgaatatgtggcgaatgggcccgcgcccggtcgact                                      <  2:1055177‑M1/134‑92 (MQ=255)
                                                                                              gCCATGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctga                               >  1:413323‑M1/1‑36 (MQ=255)
                                                                                                  tGAGCAACCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggat                                                                                                       <  2:1013410‑M1/58‑27 (MQ=255)
                                                                                                  tGAGCAACCGGTGCAAGGTCGACAGCCGCTGGCGCTGCACGGGCTGCGGCGCCTGGAt                                                                                                       >  1:1013410/1‑58 (MQ=21)
                                                                                                        aCCGGTGCAAGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggc                     >  1:32572‑M1/1‑26 (MQ=255)
                                                                                                                 aGGTCGACAGCCGCTGGcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcagg             >  2:676943‑M1/1‑17 (MQ=255)
                                                                                                                         agccgctggcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcaggcatggccgc    <  2:413323‑M1/130‑126 (MQ=255)
                                                                                                                           ccgctggcgctgcacgggctgcggcgcctggatgtaccagagcacggcattgggtgagccgccgaatatgtggcgaatgggcccgcgcccggtcgactggcctgataatgacggctgagcaggcatggccgcgc  >  2:1112376‑M1/3‑7 (MQ=255)
                                                                                                                                 |                                                                                                                               
GATGCCATGCACGCCTGTGACCTGTGTCATCAGCTGGTGGGCGAACCGTCCAGCAGTGCCGCCACACGAGCACCTGGCTGCCTCTGGCATTGCGGCCATGAGCAACCGGTGCAAGGTCGACAACTGCTGGTGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGCTGCACGGGCTGCGGCGCCTGGATGTACCAGAGCACGGCATTGGGTGAGCCGCCGAATATGTGGCGAATGGGCCCGCGCCCGGTCGAC  >  AM260479/1464311‑1464567

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: