Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,773,173 Δ1 bp coding (893/1914 nt) mtlA → mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,773,1730G.94.0% 192.4 / 4.5 51coding (893/1914 nt)mtlAmannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base . (30/17);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATC  >  NC_000913/3772940‑3773420
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
cccAGGTCCAGGTATGGGCGTGCTGCTGGCGTTCATTTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCTGTGCGGCAATCATCCACTTCTTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTa                                                                                                                                                                                                                                          <  2:208607/249‑1 (MQ=255)
   aggtccaggtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCATGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGCATATCCTGGGCGGGGGTCTGGGTTTGCCGGCGTCTCCGGTTTTATCCTTTCTTTACTGGc                                                                                                                                                                                                                                     >  2:234766/1‑251 (MQ=255)
           gtATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTTGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGacac                                                                                                                                                                                                                             >  1:321322/1‑251 (MQ=255)
                        gctgGCGTACATGTTCTTTGGTCGTGGTAGCGTTAAACAGTCTGCTGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCtt                                                                                                                                                                                                                 <  2:82240/250‑1 (MQ=255)
                        gctgGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGt                                                                                                                                                                                                                                           >  1:392922/1‑224 (MQ=255)
                        gcgggCGTACATGTTCTTTGGGCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGt                                                                                                                                                                                                                                           <  2:392922/221‑1 (MQ=255)
                        actgGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCCCCGGTTTTATCCTTGCTGTACTGGCGCTGGCACCAAAAggggact                                                                                                                                                                                                                 >  2:166656/2‑245 (MQ=255)
                                aCATGTCCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGTGCGGTGCGGCAATCATCCACTCCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCGTGCTTTACTGGCGATGACACCAAAAGGTGCTTACTTCGCt                                                                                                                                                                                                         <  2:246364/250‑1 (MQ=255)
                                       atttGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTCTGACAGGCGTGTTCACGCTGACTGTCCTGGGCGGGGGGATGGTTTCTCCGGCATCTTCGGTTTTATACTTGCTGTGCTGGCGATGTCCCCCAAAGGTGCTTGCTTCGCTAACATcgc                                                                                                                                                                                                 >  2:100001/2‑251 (MQ=255)
                                                                  tGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGCCTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCCGGGTTCTCCGGCATCTCCGGTTTTATCCTTGATGTACTGGCGATGGCACCAACAGGGGCTTACTTCGCTACCATCGCGGGTGGGTGTGGGGCGCTGGCTGtctc                                                                                                                                                                      >  2:442253/1‑251 (MQ=255)
                                                                   gcgggcggTGCGGCAAGCATCCACTTCCTGGGGGGTATCCACGAAATCTACTGCCCGTATGAGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCGCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGtc                                                                                                                                                                        <  2:110673/248‑1 (MQ=255)
                                                                     gggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTgg                                                                                                                                                                                                                                      <  1:337974/184‑1 (MQ=255)
                                                                     gggcggTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTgg                                                                                                                                                                                                                                      >  2:337974/1‑184 (MQ=255)
                                                                                    catcCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTATTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAt                                                                                                                                                   >  1:179399/1‑251 (MQ=255)
                                                                                       ccttctCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGAGCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTAtttt                                                                                                                                                <  2:350124/246‑1 (MQ=255)
                                                                                            tCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCt                                                                                                                                             <  1:260252/249‑1 (MQ=255)
                                                                                                  ggggTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTAGCAGTCATCCTCGGCGGTATGACTGGGGTGTTCAAGCTGACTCTCCTGGGCGGGGGGCTGGGGTCTCCGGCATCTCCGGTTTTATCCTTGCTGGACTGGAGCTGACACCAAAAGGTGGGTACTTGGCTAACATCGCGGGGGTGTGTGAGGCGATGGGTGTCTCCTTCGGTGTCTCGGCCTATTTGCTGAAAAcc                                                                                                                                      >  2:140975/1‑250 (MQ=255)
                                                                                                        tccacgccaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCCAAAGGTGCTTACTTCTCTAACCTCGCGGGTGTGGGTGCGGCGCTGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAACACCAGCaa                                                                                                                                 >  2:206423/9‑249 (MQ=255)
                                                                                                               aaTCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCCTCTCCGGTTTTATCCTTGCTGTACTGGCGATGACCCCAAAAGGGGCTTACTTCGCTCACATCGCGGGGGTGTGTGCGGCGATGGCTGTCTCCTTAGTTTTCTCTGCTATTTTGCTGAAAAACAGCAAAGGGaa                                                                                                                           >  2:386036/1‑248 (MQ=255)
                                                                                                                 tCTACATCCCGAATGTGATGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGACGGCTATCCTGGGCGGTGGTCTGGTTTCTCCAGCATTTCCGGATTTATCCTTGCTGGACTGGGGTTAACACCACAAGGAGATTACTTCGCTAACATCGAGGGTGTGTGAGCGGCGATGGATGTTTCCTTGGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGGGaaaaaat                                                                                                                      >  2:314473/1‑247 (MQ=255)
                                                                                                                    aCTTCCCGGTTGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGAGTGTTCACGCTGACTATCCTGGGCTGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaagaagaag                                                                                                                   <  2:270734/251‑1 (MQ=255)
                                                                                                                       tCCCGTAGGTGCTGATGAATCCGCGTCTGATCCGCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAaagaagaag                                                                                                                   <  2:326476/248‑1 (MQ=255)
                                                                                                                        cccGTATGGGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAgatg                                                                                                                <  1:235357/250‑1 (MQ=255)
                                                                                                                            tATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCCTCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGGGCTTACTTCGCTAACATCGCGGGTGGGGGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTTCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATAtt                                                                                                           >  2:357469/1‑251 (MQ=255)
                                                                                                                               gtgCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGCTGACTCCACAAGGTGCTTTCTTCGCTCACATCGCGGGTGTGGTTGCGGCGGTGGCTGTCTCCTTCGTTGTCTCTGCTATTTTTCTGAAAACCAGCCAAGGGAAAGCAGAAGATGATAttta                                                                                                         >  2:160835/1‑248 (MQ=255)
                                                                                                                                     atgaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAagcagc                                                                                                   >  1:428873/1‑250 (MQ=255)
                                                                                                                                       gaATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGGGGGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTTAagccgaa                                                                                                  >  2:254660/1‑246 (MQ=255)
                                                                                                                                         aTCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtc                                                                                              >  1:151406/1‑251 (MQ=255)
                                                                                                                                             gcgTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACtcgtcg                                                                                          <  1:323168/251‑1 (MQ=255)
                                                                                                                                               gTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGAATATCCTGGGCGGTGGGCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACCCCAAAAGGGGCTTACTTCGCTAACCTCACGGGGGTGTGTGCGGCGATGGGTGTATCCTTTGTTGTCTCTGCTTTTTTTCTGTAAACCAGACAAGGGAACGAAGAAGTGGATATTGAAACAACAAAtcgtcct                                                                                         >  2:343449/1‑248 (MQ=255)
                                                                                                                                                 cTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATg                                                                                      >  1:301549/1‑251 (MQ=255)
                                                                                                                                                  tGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTTTTTTGCTGAAAACCCGCCAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGc                                                                                     >  2:195050/1‑251 (MQ=255)
                                                                                                                                                         cGCAGTTATCCTCGGCGGTATGACTGGCGTGGTCACGCTGACTATCCTGGGCGGTGGACTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACTTCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCTTATGCAGGACa                                                                               <  2:145817/250‑1 (MQ=255)
                                                                                                                                                          gCAGTCATCCGCGGCGGTATGACTGGCCTGTTCACGCTGACTATCCTGGGCGGTTGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATGTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATg                                                                             <  2:321322/251‑1 (MQ=255)
                                                                                                                                                              tCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGGGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGaa                                                                           >  1:80811/1‑249 (MQ=255)
                                                                                                                                                                        cggTATGCCTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTa                                                                >  2:386624/1‑250 (MQ=255)
                                                                                                                                                                                cTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAt                                                       >  2:330085/1‑251 (MQ=255)
                                                                                                                                                                                    cGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTcc                                                   >  1:247379/1‑251 (MQ=255)
                                                                                                                                                                                    cGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCAGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTATCTGCTATTTTGCTGCAAACAAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCTGGACATGCAAGCTGAGTCTAAAGGCGCATCTcc                                                   >  2:350827/1‑251 (MQ=255)
                                                                                                                                                                                        ttCACGCTGAGTATCCGGGGCGGTGGTCTGGCTTCTCCGGCATCTCAGGTTTTATCCTTGCTGTACTGGCGCTGACACCCCAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGACAACCAGCAAAGTGAAGGTAGAAGTTGATATTGACGCCGCAACTCGGCGTATGCAGGACTTGAAAGCTGAGTCTAAAGGCGCATCTCCGCTg                                               >  2:45347/1‑251 (MQ=255)
                                                                                                                                                                                          cACGCTGACTATCCGGTGCGGTTGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGTTGCTTACTTCGCTAACATCGCGGGTTTGTGTGCTGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGt                                              <  1:103497/250‑1 (MQ=255)
                                                                                                                                                                                                   tATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTATGCTGGCg                                     >  1:403138/1‑250 (MQ=255)
                                                                                                                                                                                                       ctgggaagtggaCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCAGTAATGGCGATGCCACCAAACGGTGCTTACTTCGCTAACATCGCGGGTGTGAGTGCGGCGATGGATGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAACGAAGCAGATTATATTGAAGCAGCAACTCGTCGTATGCAGGCCATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGAtg                                  >  2:310865/8‑249 (MQ=255)
                                                                                                                                                                                                        tGGGCGGTGGTCTGGTTTCTCCGGCATTTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGAtgtg                                >  1:349436/1‑250 (MQ=255)
                                                                                                                                                                                                              gtggtCTGGTTTCTCCGGCATATCCGGTTTTATCCTTGCTGTACTGGCGATGACACCACAAGGCGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCCGGACATGAAAGCTGAGTCTAAAGGCGCCTCTCCGCTGGCTGCTGGCGCTGGGGCTAACg                         >  2:82856/1‑251 (MQ=255)
                                                                                                                                                                                                                  tCTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCt                     >  1:216372/1‑251 (MQ=255)
                                                                                                                                                                                                                   cTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAt                                                       >  1:13275/1‑216 (MQ=255)
                                                                                                                                                                                                                   cTGGTTTCTCCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACGTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCGCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCAt                                                       <  2:13275/216‑1 (MQ=255)
                                                                                                                                                                                                                          ctcCGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCt                     <  1:274804/243‑1 (MQ=255)
                                                                                                                                                                                                                          ctcCGGCATCTCCGGTTTTATACTTGCTGTACTGGCGATGACACCAAACGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCt                     >  2:274804/1‑243 (MQ=255)
                                                                                                                                                                                                                             cGGCATCTCCGGTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCacgtac          >  1:148931/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGCAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAat   >  1:80641/4‑251 (MQ=255)
                                                                                                                                                                                                                                      ccggTTTTATCCTTGCTGTTCTGGCGATGACACCAAAAGTTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTaaaa    <  2:368273/247‑1 (MQ=255)
                                                                                                                                                                                                                                      ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGGAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTaaa     <  1:461946/246‑1 (MQ=255)
                                                                                                                                                                                                                                      ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTaaa     <  1:386036/246‑1 (MQ=255)
                                                                                                                                                                                                                                      ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTaa      <  1:330085/245‑1 (MQ=255)
                                                                                                                                                                                                                                      ccggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAat   <  1:372819/248‑1 (MQ=255)
                                                                                                                                                                                                                                       cggTTTTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAatc  >  1:175895/2‑250 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        
CCCAGGTCCAGGTATGGGCGTGCTGCTGGCGTACATGTTCTTTGGTCGTGGTAGCGCTAAACAGTCTGCGGGCGGTGCGGCAATCATCCACTTCCTGGGGGGTATCCACGAAATCTACTTCCCGTATGTGCTGATGAATCCGCGTCTGATCCTCGCAGTCATCCTCGGCGGTATGACTGGCGTGTTCACGCTGACTATCCTGGGCGGTGGTCTGGTTTCTCCGGCATCTCCGGGTTCTATCCTTGCTGTACTGGCGATGACACCAAAAGGTGCTTACTTCGCTAACATCGCGGGTGTGTGTGCGGCGATGGCTGTCTCCTTCGTTGTCTCTGCTATTTTGCTGAAAACCAGCAAAGTGAAAGAAGAAGATGATATTGAAGCAGCAACTCGTCGTATGCAGGACATGAAAGCTGAGTCTAAAGGCGCATCTCCGCTGTCTGCTGGCGATGTGACTAACGACCTGAGCCACGTACGTAAAATC  >  NC_000913/3772940‑3773420

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