Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 2,171,429 Δ5,958 bp IS3‑mediated [gatR]gatY [gatR], gatD, gatB, gatA, gatZ, gatY

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 2170375–2171428 2177386 5959–7012 2 [1] [0] 31 [insE5]–gatY [insE5],insF5,gatR,gatD,gatB,gatA,gatZ,gatY

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 2171428NA (NA)31 (0.960) 28/488 0.2 100% noncoding (1258/1258 nt) IS3 repeat region
?NC_000913 2177387 = 0 (0.000)intergenic (‑183/+125) gatY/fbaB tagatose‑1,6‑bisphosphate aldolase 2/fructose‑bisphosphate aldolase class I

CGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTG  >  NC_000913/2177156‑2177630
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
tgctgcaacgtgaatgcctgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCt                                                                                                                                                                                                                                   <  2:114349‑M2/19‑1 (MQ=255)
ggctgctccgataatgcctgcggggaaagcgtctttcctgctccgaaagggcattgtatccatggagaacactttatcagccgggaaataatgcgtgcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtttggcggcctcagtccggaacaatgtgaaaacaagaacctcgcttaggcctgtgtccatattatgtgggtaggatcaAAACATACGTGTCAGTGCTa                                                                                                                                                                                                                                  <  2:227460‑M2/20‑1 (MQ=255)
               gcctgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAg                                                                                                                                                                                                                   <  2:221760‑M2/35‑1 (MQ=255)
               gcctgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAg                                                                                                                                                                                                                   <  2:17590‑M2/35‑1 (MQ=255)
               gcctgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGGGCTATGTAGGCATGATAAg                                                                                                                                                                                                                   <  1:268955‑M2/35‑1 (MQ=255)
                cctgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgggggtaggataaAAACATCCGTGGCAGGGCTATGTAGGCATGAtaaa                                                                                                                                                                                                                   >  2:100136‑M2/216‑249 (MQ=255)
                  tgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagcacctcgcttaggcctgtgtccatattacgtgggaaggatcaAACCCTACGTGGCAGTGCTAAGTAGGCATGAtaagccg                                                                                                                                                                                                                >  1:188508‑M2/214‑248 (MQ=255)
                  tgcgtggaaagcttctttcattcgctgaaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACg                                                                                                                                                                                                                >  1:202530‑M2/214‑251 (MQ=255)
                                             aaagtggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTg                                                                                                                                                                                     <  2:227665‑M2/65‑1 (MQ=255)
                                                 tggaatgtatccatggagaacactttatcagccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                                                                                                                                                                                 >  1:325550‑M2/183‑251 (MQ=255)
                                                                              agccgggaaataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTAc                                                                                                                                                    >  2:79557‑M2/154‑251 (MQ=255)
                                                                                       ataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                                                                                                                                                                    >  1:36063‑M2/145‑226 (MQ=255)
                                                                                       ataatgcgggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc                                                                                                                                                                    <  2:36063‑M2/82‑1 (MQ=255)
                                                                                         aatgctggcaacggtgtttaattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaaccgcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGtgtt                                                                                                                                         <  2:11140‑M2/109‑1 (MQ=255)
                                                                                                            aattatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCa                                                                                                                      <  1:149509‑M2/128‑1 (MQ=255)
                                                                                                              ttatatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAgg                                                                                                                    >  1:92528‑M2/122‑251 (MQ=255)
                                                                                                                 tatcgaatgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGcgg                                                                                                                 >  1:155609‑M2/119‑250 (MQ=255)
                                                                                                                       atgtgattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTaa                                                                                                           <  1:8206‑M2/139‑1 (MQ=255)
                                                                                                                           gattacaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAAttt                                                                                                        <  1:14809‑M2/142‑1 (MQ=255)
                                                                                                                               acaatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcccttcggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTCCGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCt                                                                                                    >  1:119645‑M2/105‑250 (MQ=255)
                                                                                                                                 aatcggtggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATc                                                                                                 <  1:186384‑M2/149‑1 (MQ=255)
                                                                                                                                       gggcggcggcacagttggtgtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCATACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATa                                                                                           <  2:85679‑M2/155‑1 (MQ=255)
                                                                                                                                                          gtggcggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGtt                                                                         <  1:198814‑M2/173‑1 (MQ=255)
                                                                                                                                                               ggcctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAg                                                                   <  1:223019‑M2/179‑1 (MQ=255)
                                                                                                                                                                 cctcagtccggaacaatttgaaaacaagaacctcgcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGtt                                                                 >  2:87717‑M2/71‑251 (MQ=255)
                                                                                                                                                                                                   gcttaggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCttt                               <  1:49121‑M2/215‑1 (MQ=255)
                                                                                                                                                                                                      taggcctgtgtccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTAACGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTa                              <  1:79557‑M2/216‑1 (MQ=255)
                                                                                                                                                                                                                 ccatattacgtgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATc                  <  2:8567‑M2/228‑1 (MQ=255)
                                                                                                                                                                                                                           tgggtaggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGAGGCTTCAGCGTCGCATCAGGCATAGGCTGCCGGATGCCGCGTAAACGCCTTATCCGTCCTACGGGGTGGTCTTGTTCGCATAAGGCGTATCAGGAGATAGTAATTTTGCTATCGAGATAAACGTCCAGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTcc       >  1:304240‑M2/13‑251 (MQ=255)
                                                                                                                                                                                                                               taggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCt   >  1:227919‑M2/9‑251 (MQ=255)
                                                                                                                                                                                                                                aggatcaAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTAAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTg  >  1:302014‑M2/8‑251 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                   
CGCCGCGCTGTGCGTTGTTCAGCATCTGCTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTCCTACGGGGTGGTGTTGTTCGTATAAGGCGTATCAGGCGATAGTAATTTTGCTATCGAGATAAACGTCCTGCACGGCGTTAATCAGTTTCACGCCGTCAGCCATCGATTTCTTGAACGCTTTACGTCCAAGAATCAGCCCCATTCCGCCTG  >  NC_000913/2177156‑2177630

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