Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,361,226 Δ9,811 bp [gltF][nanE] [gltF], yhcA, yhcD, yhcE, insH10, yhcE, yhcF, yhcG, nanQ, nanK, [nanE]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3361226 3371036 9811 32 [0] [0] 33 [gltF]–[nanE] [gltF],yhcA,yhcD,yhcE,insH10,yhcE,yhcF,yhcG,nanQ,nanK,[nanE]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 33612250 (0.000)32 (0.560) 31/490 0.5 100% coding (50/765 nt) gltF periplasmic protein GltF
?NC_000913 3371037 = 0 (0.000)intergenic (‑1/+47) nanE/nanT putative N‑acetylmannosamine‑6‑phosphate 2‑epimerase/N‑acetylneuraminate:H(+) symporter

TCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCAGAT  >  NC_000913/3370801‑3371286
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         
caccacaaaccacctcaccagcaccttattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGc                                                                                                                                                                                                                                             <  1:78433‑M2/15‑1 (MQ=255)
               caccagcaccttattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTAc                                                                                                                                                                                                                               <  2:373900‑M2/29‑1 (MQ=255)
                          tattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                                                                                                                                                                                    >  2:436854‑M2/211‑250 (MQ=255)
                          tattattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                                                                                                                                                                                    <  2:243506‑M2/40‑1 (MQ=255)
                             tattcacatgatcgtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAAc                                                                                                                                                                                                                 >  2:250652‑M2/208‑250 (MQ=255)
                                          gtaactataaaattcattaaataaatacatcccagcaagttggctatttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTg                                                                                                                                                                                                   <  2:287924‑M2/57‑1 (MQ=255)
                                                                                        tttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGTACCGc                                                                                                                                                     >  1:450443‑M2/149‑251 (MQ=255)
                                                                                         ttctgaacgagttatttatttacattcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGc                                                                                                                                                     <  2:122225‑M2/103‑1 (MQ=255)
                                                                                                                  tcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGccc                                                                                                                                                                                                                                    <  1:102235‑M2/24‑1 (MQ=255)
                                                                                                                  tcacaaacagatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGccc                                                                                                                                                                                                                                    >  2:102235‑M2/123‑146 (MQ=255)
                                                                                                                           gatatatatactgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTAATGAGTCCGc                                                                                                                  >  2:140899‑M2/114‑251 (MQ=255)
                                                                                                                                  atactgcttcagctttttattttaaaaaataaatgtgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCtt                                                                                                           <  1:69543‑M2/145‑1 (MQ=255)
                                                                                                                                     ctgcttcagctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCg                                                                                                        >  1:176158‑M2/104‑251 (MQ=255)
                                                                                                                                            agctttttattttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGca                                                                                                  >  1:437085‑M2/97‑250 (MQ=255)
                                                                                                                                                    attttaaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGc                                                                                          <  1:436854‑M2/162‑1 (MQ=255)
                                                                                                                                                         aaaaaataaatttgttattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAc                                                                                     <  2:176158‑M2/167‑1 (MQ=255)
                                                                                                                                                                        tattaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGTACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATc                                                                      <  2:450443‑M2/182‑1 (MQ=255)
                                                                                                                                                                         attaaggatatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGNAACCAACGCTGACCGCGAGAAGGCATATcc                                                                     >  2:210795‑M2/68‑250 (MQ=255)
                                                                                                                                                                            aaggatatgatcatatgtttttcaaaaagaaactcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATACCCAGcc                                                                 >  2:424667‑M2/65‑251 (MQ=255)
                                                                                                                                                                                atatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGCATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAAt                                                             >  2:144839‑M2/61‑251 (MQ=255)
                                                                                                                                                                                 tatgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGGACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGCGAAGGCATATACAGCCCAATc                                                            >  1:156745‑M2/60‑251 (MQ=255)
                                                                                                                                                                                  atgttcatatgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATc                                                            >  1:379953‑M2/59‑250 (MQ=255)
                                                                                                                                                                                          atgtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGAt                                                    <  1:392353‑M2/200‑1 (MQ=255)
                                                                                                                                                                                            gtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATcac                                                 >  1:5796‑M2/49‑251 (MQ=255)
                                                                                                                                                                                            gtttttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATcac                                                 <  2:372845‑M2/203‑1 (MQ=255)
                                                                                                                                                                                                ttcaaaaagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATcaccac                                              <  1:375553‑M2/206‑1 (MQ=255)
                                                                                                                                                                                                      aagaacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTc                                       >  1:378319‑M2/39‑251 (MQ=255)
                                                                                                                                                                                                        gaacctaacaacagctgctatttgtgaagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGTTTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGACCGTc                                       >  2:408324‑M2/37‑249 (MQ=255)
                                                                                                                                                                                                         aacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTc                                       >  2:242599‑M2/36‑248 (MQ=255)
                                                                                                                                                                                                         aacctcacaacagctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTc                                       >  1:112047‑M2/36‑248 (MQ=255)
                                                                                                                                                                                                                      ctgctatttgtgcagcactttcACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCg                       >  2:177910‑M2/23‑251 (MQ=255)
                                                                                                                                                                                                                                           cACAGCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTATCTTTTTGTTTTGACTAAATCGTTTTTGGCGCAGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGACATGAGTACGCAACGCTTCCGGGCGGAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCAGAt  <  1:290861‑M2/250‑1 (MQ=255)
                                                                                                                                                                                                                                            acatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTGACTTTTGGTTTTGACTAAATCGTTTTTGGCGCAGCCAAACGGCACGGCAAAGCTGAATGGTTTACCGTCGATAGCGTCATGAGGACGCAACGCTTCCGGGCGCAACCAACGCTGAAAGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGCGAGAGATGACAGCGCAGTAccc      >  1:254022/1‑246 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         
TCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTGGCTTGCAGATTTGCCACACCTTCAATGCGAATGGCAACCGCGCCCGCCTGTTCTGCCGCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATATCCAGCCCAATCAGCAGGATCACCACGAACGTCAGACTGAACGAGAGCGATGCCAGCGCAGTACCCAGAT  >  NC_000913/3370801‑3371286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ 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.