Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 2,034,358 Δ19,491 bp 100% [insB‑14][wcaJ] 20 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 2034358–2035096 2053848 18753–19491 22 [0] [0] 22 [insB‑14]–[wcaJ] [insB‑14],insA‑14,wbbD,wbbC,wzy,wbbB,wbbA,vioB,vioA,wzx,rmlC,rfbA,rfbD,rfbB,galF,wcaM,wcaL,wcaK,wzxC,[wcaJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 20343570 (0.000)18 (0.670) 16/278 NT 100% noncoding (737/768 nt) IS1 repeat region
?REL606 2053849 = 0 (0.000)coding (823/1395 nt) wcaJ predicted UDP‑glucose lipid carrier transferase

GTACCATTACCACCACTTAATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGC  >  REL606/2034209‑2034506
                                                                                                                                                    |                                                                                                                                                     
gTACCATTACCACCACTTAATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAg                                                                                                                                                      <  2:92472‑M1/150‑2 (MQ=255)
  aCCATTACCACCACTTAATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtt                                                                                                                                                    <  2:442956‑M1/150‑4 (MQ=255)
  aCCATTACCACCACTTAATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtt                                                                                                                                                    <  2:243884‑M1/150‑4 (MQ=255)
              acTTAATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCAAGACAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggt                                                                                                                                        >  1:417331‑M1/1‑135 (MQ=255)
                   atatAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaacc                                                                                                                                   >  2:186321‑M1/1‑130 (MQ=255)
                   atatAAACTGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATCTATCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaacc                                                                                                                                   <  1:177379‑M1/150‑21 (MQ=255)
                          cGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaa                                                                                                                            <  2:487071‑M1/150‑28 (MQ=255)
                             aaTAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcg                                                                                                                         >  1:137990‑M1/1‑120 (MQ=255)
                               tAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTATAACAAATATCCGGTAATGACTACAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggcaagcggc                                                                                                                       >  1:352365‑M1/1‑118 (MQ=255)
                                        aGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatac                                                                                                              <  2:372319‑M1/150‑42 (MQ=255)
                                        aGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatac                                                                                                              <  1:56709‑M1/150‑42 (MQ=255)
                                          cTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacag                                                                                                            <  1:353495‑M1/150‑44 (MQ=255)
                                                                   tttATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatc                                                                                   <  1:125018‑M1/150‑69 (MQ=255)
                                                                     tatatatataGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctc                                                                                 <  2:550672‑M1/150‑71 (MQ=255)
                                                                                                      gTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaaggttgaaggtaaa                                                >  1:557197‑M1/1‑47 (MQ=255)
                                                                                                         tAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaatgttgaaggtaaagac                                             <  1:358765‑M1/150‑107 (MQ=255)
                                                                                                         gacagaatatccGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaatgttgaaggtaaagac                                             <  1:526508‑M1/145‑107 (MQ=255)
                                                                                                              aaTATCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaatgttgaaggtaaagacctgtc                                        >  2:528024‑M1/1‑39 (MQ=255)
                                                                                                                 atCCGGTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagtggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaatgttgaaggtaaagacatcgggga                                     <  2:394828‑M1/150‑115 (MQ=255)
                                                                                                                      gTAATGACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaatgttgaaggtaaagacatcggggatcagc                                <  2:138384‑M1/150‑120 (MQ=255)
                                                                                                                           gACTCCAACTTATTGATAGTGTTTTAgtttgagcaggcggttaaccccggaaagcggcgtgtcatacagcggcaccaccggtacgccattcatctcttcgaggcgtgaatggagaatgttgaaggtaaagacatcggggaacagcagcac                           <  2:461904‑M1/150‑125 (MQ=255)
                                                                                                                                                    agtttgcgcaggaggttaaccccggaacgctgcgtgtcatacagcggcaccaccggtacgccattcatctcttggaggcgtgaatggagaatgtggaaggtaaagacatcggggatcagcagcaccgcacaggtggtctccgccatttta  >  1:302942‑M1/1‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTACCATTACCACCACTTAATATAACCGGAATAAGAGGTGAGCTCATAAAAATCTCCTGATAACTTTTTTATATATATAGCCACGAAAACTAGTGGGAAGAAGTCTAACAAATATCCGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGC  >  REL606/2034209‑2034506

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