Predicted mutation
evidence seq id position mutation annotation gene description
MC JC AM260480 2,912,452 Δ39 bp intergenic (+243/–) h16_B2579 → /  cold shock protein, DNA‑binding/–

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ AM260480 2912452 2912490 39 20 [0] [0] NA h16_B2579/– cold shock protein, DNA‑binding/–

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AM260480 1 =0 (0.000)13 (0.690) 10/130 0.6 100% intergenic (–/‑1344) –/pA2 –/replication initiation protein RepA
?AM260480 = 2912451 0 (0.000)intergenic (+242/–) h16_B2579/– cold shock protein, DNA‑binding/–

ACCGCCAGGCGCTCCGCTCCAGCCCCCGCATGCCCATGCCGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCTGCTGCACCGCACGAGGTCAGGCACGGCCCGTGAGC  >  AM260480/2912330‑2912490
                                                                                                                         |                                                                                                         
aCCGCCAGGCGCTCCGCTCCAGCCCCCGCATGCCCATGCCGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcacc                                                                                               <  1:349795‑M2/134‑13 (MQ=255)
aCCGCCAGGCGCTCCGCTCCAGCCCCCGCATGCCCATGCCGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcac                                                                                                <  2:119425‑M2/133‑12 (MQ=255)
aCCGCCAGGCGCTCCGCTCCAGCCCCCGCATGCCCATGCCGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcac                                                                                                >  1:119425‑M2/1‑122 (MQ=255)
                            catgcccatgccGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcc                                                                   <  1:962343‑M2/134‑41 (MQ=255)
                              tgcccatgccGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctc                                                                 <  2:64794‑M2/134‑43 (MQ=255)
                                  catgccGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcag                                                             >  2:962343‑M2/1‑88 (MQ=255)
                                      gcaaggccagccgGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaat                                                         >  1:1176345‑M2/6‑84 (MQ=255)
                                           gccAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggta                                                    >  1:730458‑M2/1‑79 (MQ=255)
                                                        ccTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggtacggtctgtcctgg                                       >  1:64794‑M2/1‑66 (MQ=255)
                                                                                            cAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtca                                                                                    >  2:1172857‑M2/1‑30 (MQ=255)
                                                                                            cAGCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtca                                                                                    <  1:1172857‑M2/53‑24 (MQ=255)
                                                                                              gCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggtacggtctgtcctggcaggcctgcttgcgctgggcagcgggtacgcagaaa   <  2:910623‑M2/132‑105 (MQ=255)
                                                                                              gCCTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggtacggtctgtcctggcaggcctgcttgcgctgggcagcgggtacgcagaaa   <  1:1003149‑M2/132‑105 (MQ=255)
                                                                                               ccTTGCGCTCCACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggtacggtctgtcctggcaggcctgcttgcgctgggcagcgggtacgcagaaac  >  2:1003149‑M2/1‑27 (MQ=255)
                                                                                                        ccACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggta                                                    >  1:860635‑M2/1‑18 (MQ=255)
                                                                                                        ccACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcctcgcagaaatgggta                                                    <  2:860635‑M2/73‑56 (MQ=255)
                                                                                                        ccACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcct                                                                  >  1:360897‑M2/1‑18 (MQ=255)
                                                                                                        ccACGCAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagcct                                                                  <  2:360897‑M2/59‑42 (MQ=255)
                                                                                                             cAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagc                                                                    <  1:745324‑M2/52‑40 (MQ=255)
                                                                                                             cAGAATCGGGTACgcctgctgcaccgcacgaggtcaggcacggcccgtgagc                                                                    >  2:745324‑M2/1‑13 (MQ=255)
                                                                                                                         |                                                                                                         
ACCGCCAGGCGCTCCGCTCCAGCCCCCGCATGCCCATGCCGCCGCCAGCCGGCACGCCTGGCCTTGCGCGGCGCACTTTCCCGCTCACCATCCAGCCTTGCGCTCCACGCAGAATCGGGTAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCCTGCTGCACCGCACGAGGTCAGGCACGGCCCGTGAGC  >  AM260480/2912330‑2912490

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