Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,196,367 A→T D129E (GAT→GAA hldE ← heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,196,3670AT97.4% 431.4 / 2.5 113D129E (GAT→GAAhldEheptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base T (38/72);  total (40/73)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

Only 100 of 113 total aligned reads displayed.

TGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCG  >  NC_000913/3196220‑3196514
                                                                                                                                                   |                                                                                                                                                   
tGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCGTGATTAATCCGCTCGTGCAGCGGCTGCGGTTc                                                                                                                                                   <  2:1565930/150‑1 (MQ=255)
tGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTc                                                                                                                                                   <  2:70506/150‑1 (MQ=255)
aggATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATc                                                                                                                                                   >  2:1450685/2‑150 (MQ=255)
 ggATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCaa                                                                                                                                                 <  2:1008329/151‑1 (MQ=255)
    tcaatcaGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAacac                                                                                                                                              >  1:433457/1‑151 (MQ=255)
        tcaGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTc                                                                                                                                          >  1:50372/1‑151 (MQ=255)
          aGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGa                                                                                                                                        <  2:956257/151‑1 (MQ=255)
           gCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGaa                                                                                                                                       <  1:1649495/151‑1 (MQ=255)
           gCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCCGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGa                                                                                                                                        <  2:278489/150‑1 (MQ=255)
            cACCGGAACACCCGCCTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGaa                                                                                                                                       <  2:867392/150‑1 (MQ=255)
              ccGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAAcc                                                                                                                                    <  1:377182/151‑1 (MQ=255)
                 gAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACcttc                                                                                                                                 <  1:24569/151‑1 (MQ=255)
                     acCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCttcttc                                                                                                                              >  1:998356/1‑150 (MQ=255)
                       ccGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCaaa                                                                                                                           <  2:1122394/151‑1 (MQ=255)
                        cGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCaaaa                                                                                                                          <  1:562184/151‑1 (MQ=255)
                        cGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCaaaa                                                                                                                          <  1:4916/151‑1 (MQ=255)
                         gCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACTCCTTCGAAACCTTCTTCaaaa                                                                                                                          <  2:462105/150‑1 (MQ=255)
                         gCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAAt                                                                                                                         >  1:1250493/1‑151 (MQ=255)
                         gCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAAt                                                                                                                         >  1:1050017/1‑151 (MQ=255)
                            ttACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATcc                                                                                                                       <  2:820589/150‑1 (MQ=255)
                            ttACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATcc                                                                                                                       >  1:409712/1‑150 (MQ=255)
                              aCGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCACCACCTTCGAAACCTTCTTCAAAATCCAGa                                                                                                                    <  2:433426/151‑1 (MQ=255)
                              aCGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGa                                                                                                                    <  2:1042953/151‑1 (MQ=255)
                                 cgcCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACgg                                                                                                                 <  2:117137/151‑1 (MQ=255)
                                  gcCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGa                                                                                                                <  1:152538/151‑1 (MQ=255)
                                   ccAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGAt                                                                                                               <  2:1119897/151‑1 (MQ=255)
                                    cAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATc                                                                                                              >  1:450356/1‑151 (MQ=255)
                                     aGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCa                                                                                                             <  2:7148/151‑1 (MQ=255)
                                     aGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCa                                                                                                             <  1:668330/151‑1 (MQ=255)
                                      gTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATc                                                                                                              >  1:1060265/1‑149 (MQ=255)
                                             catcTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTgttg                                                                                                      <  1:1127428/150‑1 (MQ=255)
                                                  gctgTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTgttggtt                                                                                                   <  1:1232494/148‑1 (MQ=255)
                                                  gctgTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGc                                                                                                 <  2:1371944/150‑1 (MQ=255)
                                                  gctgTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCg                                                                                                <  2:1505321/151‑1 (MQ=255)
                                                   ctgTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCgg                                                                                               >  1:1654367/1‑151 (MQ=255)
                                                         gCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGt                                                                                         >  1:1059642/1‑151 (MQ=255)
                                                           ttGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTAc                                                                                       >  1:919637/1‑151 (MQ=255)
                                                           ttGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTAc                                                                                       >  2:1623898/1‑151 (MQ=255)
                                                            tGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTAc                                                                                       <  1:263573/150‑1 (MQ=255)
                                                             gCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTAccc                                                                                     <  1:481734/151‑1 (MQ=255)
                                                              ccAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCg                                                                                    <  1:1337108/151‑1 (MQ=255)
                                                               cAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGt                                                                                   >  1:1350551/1‑151 (MQ=255)
                                                               cAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGt                                                                                   <  2:1645199/151‑1 (MQ=255)
                                                               cAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGt                                                                                   >  1:1200963/1‑151 (MQ=255)
                                                                  cgcACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTCCCCGTaa                                                                                 <  1:274477/150‑1 (MQ=255)
                                                                  cgcACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAAt                                                                                >  2:731606/1‑151 (MQ=255)
                                                                  cgcACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAAGCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAAt                                                                                <  1:262958/151‑1 (MQ=255)
                                                                    cACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAAttt                                                                              >  1:514372/1‑151 (MQ=255)
                                                                         ttGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGt                                                                           <  2:408073/149‑1 (MQ=255)
                                                                         ttGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTa                                                                          <  1:1066799/150‑1 (MQ=255)
                                                                         ttGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTa                                                                          >  2:74369/1‑150 (MQ=255)
                                                                         ttGGCGTAGTCAGAAAGCACCAGCGCGCAAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTa                                                                          >  1:1538383/1‑150 (MQ=255)
                                                                          tGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTaa                                                                         <  2:1313909/150‑1 (MQ=255)
                                                                          tGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAAt                                                                        >  1:1106005/1‑151 (MQ=255)
                                                                          tGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAAt                                                                        >  1:1553765/1‑151 (MQ=255)
                                                                            gCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATgg                                                                      <  1:182057/151‑1 (MQ=255)
                                                                                aGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCgg                                                                  >  2:711913/1‑151 (MQ=255)
                                                                                aGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCg                                                                   <  2:1553861/150‑1 (MQ=255)
                                                                                    aGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGc                                                              <  2:1348817/151‑1 (MQ=255)
                                                                                     gAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGc                                                              <  1:1712222/150‑1 (MQ=255)
                                                                                     gAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCg                                                             >  1:1204375/1‑151 (MQ=255)
                                                                                        aGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCg                                                          <  2:761493/151‑1 (MQ=255)
                                                                                         gCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCgg                                                         <  2:1067090/151‑1 (MQ=255)
                                                                                         gCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCgg                                                         <  2:83997/151‑1 (MQ=255)
                                                                                            ccAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTAc                                                      <  2:1133731/151‑1 (MQ=255)
                                                                                            ccAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTAc                                                      >  2:581479/1‑151 (MQ=255)
                                                                                             cAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTAc                                                      >  1:1215456/1‑150 (MQ=255)
                                                                                             cAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACa                                                     <  2:1108608/151‑1 (MQ=255)
                                                                                             cAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACa                                                     <  2:1330675/151‑1 (MQ=255)
                                                                                             cAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACa                                                     <  1:1080003/151‑1 (MQ=255)
                                                                                              aGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTAc                                                      >  1:607487/1‑149 (MQ=255)
                                                                                                cgcgcCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGaa                                                  <  1:939873/151‑1 (MQ=255)
                                                                                                   gcCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACg                                               >  1:209494/1‑151 (MQ=255)
                                                                                                         cGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTc                                          <  2:959888/150‑1 (MQ=255)
                                                                                                         cGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTc                                          <  2:1250590/150‑1 (MQ=255)
                                                                                                         cGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCg                                         <  2:191310/151‑1 (MQ=255)
                                                                                                          gAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGc                                        >  1:590564/1‑151 (MQ=255)
                                                                                                           aaCTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCa                                       <  1:1666216/151‑1 (MQ=255)
                                                                                                            aCTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCAt                                      <  2:148608/151‑1 (MQ=255)
                                                                                                            aCTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCAt                                      >  2:754369/1‑151 (MQ=255)
                                                                                                               cAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCAttt                                    <  1:1654342/150‑1 (MQ=255)
                                                                                                                   gcCTGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATttgacgt                               <  2:930404/151‑1 (MQ=255)
                                                                                                                      tGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATttgacgttg                             >  2:1070071/1‑150 (MQ=255)
                                                                                                                      tGATTAATCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATttgacgtt                              >  1:1480725/1‑149 (MQ=255)
                                                                                                                           aaTCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTc                        <  2:1313602/150‑1 (MQ=255)
                                                                                                                           aaTCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCg                       <  1:924285/151‑1 (MQ=255)
                                                                                                                            aTCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTc                        <  1:1219465/149‑1 (MQ=255)
                                                                                                                            aTCCGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCgg                      <  1:754928/151‑1 (MQ=255)
                                                                                                                              ccGCTCGTGCAGCGGCTGCGGTTCAACACCTTTGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGc                     <  2:948086/150‑1 (MQ=255)
                                                                                                                              ccGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGc                     <  1:55005/150‑1 (MQ=255)
                                                                                                                               cGCTCGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCa                   <  2:1278661/151‑1 (MQ=255)
                                                                                                                                   cGTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTGCCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCaga                 <  2:274939/149‑1 (MQ=255)
                                                                                                                                    gTGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCag                  >  1:1094606/1‑147 (MQ=255)
                                                                                                                                     tGCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGAtt             >  1:540745/1‑151 (MQ=255)
                                                                                                                                      gCAGCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGAttt            <  1:971343/151‑1 (MQ=255)
                                                                                                                                         gCGGCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACt         >  2:1494833/1‑151 (MQ=255)
                                                                                                                                           ggCTGCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCa       <  2:186138/151‑1 (MQ=255)
                                                                                                                                               gCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAgcgc   <  1:711855/151‑1 (MQ=255)
                                                                                                                                               gCGGTTCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAgc     <  2:409744/149‑1 (MQ=255)
                                                                                                                                                cGGTTCAACACCTTCGAAACCTTCTTCAACATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGCCGGCCAGAGATTTACTCAgcgcg  >  2:1063651/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGGATCAATCAGCACCGGAACACCCGCTTTACGCGCCAGTTGGATCATCTGCTGTACGCTTGCCAGCGCACCTTTGGCGTAGTCAGAAAGCACCAGCGCGCCAATCGAACTCAGCGCCTGATTAATCCGCTCGTGCAGCGGCTGCGGATCAACACCTTCGAAACCTTCTTCAAAATCCAGACGGATCAGCTGTTGGTTGCGGGAAAGTACCCGTAATTTGGTAATGGTCGGATGCGTCGGTACAGAAACGAAGTCGCATTTGACGTTGACGTCGGCCAGAGATTTACTCAGCGCG  >  NC_000913/3196220‑3196514

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