Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,972,673 A→C 8.8% E51A (GAG→GCG) ‡ rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6730AC 8.8% 324.9 / 17.4 131E51A (GAG→GCG) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base A (56/62);  new base C (9/3);  total (66/65)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01

Only 100 of 138 total aligned reads displayed.

GAGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAACC  >  NC_000913/3972524‑3972821
                                                                                                                                                     |                                                                                                                                                    
gagaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcga                                                                                                                                                      >  1:206700/1‑150 (MQ=255)
  gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    >  2:658238/1‑150 (MQ=255)
  gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    >  1:188598/1‑150 (MQ=255)
  gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    >  2:304568/1‑150 (MQ=255)
     aaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCt                                                                                                                                                 <  1:787178/150‑1 (MQ=255)
       aTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCttt                                                                                                                                               >  1:998686/1‑150 (MQ=255)
        ttCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTg                                                                                                                                              <  2:982857/150‑1 (MQ=255)
         tCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGc                                                                                                                                             <  2:1569359/150‑1 (MQ=255)
          cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGCCGCGGTGGTAGTGGTCGTTAAGCTGACCTCCGCCGGAAACCTGATGTCGGTGGCACCGGTCGCGCCAAGCGCGCGGTTTGcc                                                                                                                                            >  1:1661368/1‑150 (MQ=255)
          cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                            <  1:115454/150‑1 (MQ=255)
          cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAACCAGCCCTTTGcc                                                                                                                                            >  1:1426059/1‑150 (MQ=255)
          cTGATAACAGGGGGGGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                            <  1:434249/150‑1 (MQ=255)
           tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAACCGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcgc                                                                                                                                           <  2:1522382/150‑3 (MQ=255)
           tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcgc                                                                                                                                           <  2:1361274/150‑3 (MQ=255)
           tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCGTTGCCt                                                                                                                                           >  1:321820/1‑150 (MQ=255)
           tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGAGCGTTATATCATCAACGAATCGAGCGACGCGGTGGTAGTGGTCGATAAGCTGCCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                           >  2:269168/1‑150 (MQ=255)
             aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGAGGGGGTGGTAGTGGTCGATAAGCTGACCTAAGCCGGAAGCCTGATGTCGCTGGCTCCGGTCGCGCAAAGAGAGAGCTTTGCCttt                                                                                                                                         >  1:928007/1‑150 (MQ=255)
              tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                        <  2:699618/150‑1 (MQ=255)
               accAGGTGGTGCCGGGTGTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAACCGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                       <  1:1647657/148‑1 (MQ=255)
               aaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                       >  1:1173204/1‑150 (MQ=255)
                  aGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAACGACCGCATTGACTTTGAGaa                                                                                                                                    >  2:1677130/1‑150 (MQ=255)
                     tggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCTAAGCGAGCGCGTTGCCTTTGAGAAAGt                                                                                                                                 >  1:1548407/1‑150 (MQ=255)
                     tggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGt                                                                                                                                 >  1:65889/1‑150 (MQ=255)
                      ggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGCGCGACGCGGGGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGGCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGtt                                                                                                                                <  2:646786/150‑1 (MQ=255)
                      ggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGtt                                                                                                                                >  2:1181199/1‑150 (MQ=255)
                      ggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGtt                                                                                                                                <  1:1425834/150‑1 (MQ=255)
                       gtgCCGGGTTTATTGGCTCGGCGCTTGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTg                                                                                                                               <  1:1484828/150‑1 (MQ=255)
                       gtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTg                                                                                                                               <  1:1491550/150‑1 (MQ=255)
                        tgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGa                                                                                                                              <  1:181781/150‑1 (MQ=255)
                         gCCGGGTTTATTGGCTCGGCGCTGGGGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGat                                                                                                                             <  1:1284235/150‑1 (MQ=255)
                          ccGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGata                                                                                                                            <  1:539382/150‑1 (MQ=255)
                           cGGGTTTATTGGCTCGGCGATGGTGCGTTATATCATCAACGCAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGGCGCTGGCACCGGTCGCGCAAATCGAGCGCGGTGGCGATTAGAAACTTGatat                                                                                                                           >  2:509202/1‑150 (MQ=255)
                               tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAACCGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGGCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGc                                                                                                                       <  2:1096550/150‑1 (MQ=255)
                               tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTAAGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGc                                                                                                                       <  2:1322546/150‑1 (MQ=255)
                                  aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGAt                                                                                                                    >  1:283654/1‑150 (MQ=255)
                                       cTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGc                                                                                                               >  1:91672/1‑150 (MQ=255)
                                         cGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAg                                                                                                             <  2:1713823/150‑1 (MQ=255)
                                         cGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAg                                                                                                             <  1:304568/150‑1 (MQ=255)
                                          ggCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCCAAGCGGGCGGGTTGGCGGTGTGAAAGTTTTTTTATGCGATCGGGccgc                                                                                                            >  1:1304406/1‑147 (MQ=255)
                                          agcgcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTACTCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGa                                                                                                            <  2:488514/149‑1 (MQ=255)
                                           gcgcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCCACCGCTATGCCTTAGAGAAAGTTGTTAACTGCCCTAGGGGAGaa                                                                                                           >  1:374797/1‑150 (MQ=255)
                                                  tGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAc                                                                                                    <  1:1124663/150‑1 (MQ=255)
                                                     gTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                 <  1:1170137/150‑1 (MQ=255)
                                                          ttcatcaACGAAACGCGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTc                                                                                            <  1:298247/149‑1 (MQ=255)
                                                          atcatcTAAGAAACGAGCGACGCGGCGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTc                                                                                            <  2:174513/150‑1 (MQ=255)
                                                          atcatcAACGAAACGAGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTACGCCGGAACCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTc                                                                                            <  2:1268707/150‑1 (MQ=255)
                                                            catcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc                                                                                          >  2:411182/1‑150 (MQ=255)
                                                             atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACt                                                                                         >  1:1146828/1‑150 (MQ=255)
                                                             atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACt                                                                                         <  1:1369465/150‑1 (MQ=255)
                                                             atcaACGAAACGAGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACt                                                                                         <  2:593757/150‑1 (MQ=255)
                                                                   gaaacgaGCGACGCGGTGGTAGTGGGCGATAAGCTGACCTGCTCCGGAAACCTGATGACGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCAt                                                                                   <  2:852818/150‑1 (MQ=255)
                                                                    aaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCCTc                                                                                  >  2:291029/1‑150 (MQ=255)
                                                                     aacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCa                                                                                 >  1:1023159/1‑150 (MQ=255)
                                                                       cgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATGTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGc                                                                               >  1:241640/1‑150 (MQ=255)
                                                                        gtgCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGATACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCTCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGATCTGGCACGCGTATTCACTGAGCATCAGcc                                                                              <  1:855544/148‑1 (MQ=255)
                                                                         aGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCCGCCa                                                                             >  2:862546/1‑150 (MQ=255)
                                                                          gCGACGCGGTGGTAGTGGACGATAAGCTGACCTAAGCCGGAAACCTGATGTCGCTGGCTCCGGTCGCGCAAAGCGAGCGCTTAGCCTTTGAGAAAGTTGATCTCTGCGGTCGGGCAGAACTGGCACGCGTATTAAGAGAGCATCAGcgat                                                                            >  2:267109/1‑147 (MQ=255)
                                                                           cGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCCAGGGCTTTGGCGTTGATTAATGTTATATATTCGATAAGGCCGAAAAGGCCCGCGGATTAAATGAGCATCAGCCacc                                                                           >  2:1675147/1‑148 (MQ=255)
                                                                            gACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGAc                                                                          >  2:392408/1‑150 (MQ=255)
                                                                              cgcgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtg                                                                        >  2:225877/1‑150 (MQ=255)
                                                                                cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCTGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACTCGTATTCACTGAGCATCAGCGAGACtgtg                                                                      <  1:829933/150‑1 (MQ=255)
                                                                                   tggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCa                                                                   <  2:998686/150‑1 (MQ=255)
                                                                                    ggtAGTGGTCGATACGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAt                                                                  <  2:108624/150‑1 (MQ=255)
                                                                                      tAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGCGAAAGTTGATATCTGCGATCGGGCAGCACTGGCACGCGTATTCAGTGAGCATCAGCCAGCCTGTGTCCTGc                                                                >  1:48605/1‑150 (MQ=255)
                                                                                          ggTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATATGCGATAGGGCAGAACTCGCACGCGTATTCACTGAGCATCAGCCAGCAGGTGTAATGCAACt                                                            >  2:434277/1‑150 (MQ=255)
                                                                                           gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGGCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTg                                                           >  1:1597571/1‑150 (MQ=255)
                                                                                           gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTg                                                           >  2:1495276/1‑150 (MQ=255)
                                                                                             cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCCCGCGTATTCCATGAGCCTCAGCCCGACTGTGTCATGCAAATGGc                                                         >  2:1546828/1‑150 (MQ=255)
                                                                                                  aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACGGGCCCGCGTATTCACTGAGCATCCGCCAGACTGTGTGATGCATCTGGCAGCCg                                                    >  2:828541/1‑150 (MQ=255)
                                                                                                   gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCCGCCGa                                                   >  2:269391/1‑150 (MQ=255)
                                                                                                    cTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaa                                                  >  1:223589/1‑150 (MQ=255)
                                                                                                       aCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGc                                               <  2:206700/150‑1 (MQ=255)
                                                                                                          tACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                            <  2:1698785/150‑1 (MQ=255)
                                                                                                          tACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                            <  2:1269760/150‑1 (MQ=255)
                                                                                                           aCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCCGCCGAAAGCCATg                                           >  2:334784/1‑150 (MQ=255)
                                                                                                           aCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATg                                           >  2:1180387/1‑150 (MQ=255)
                                                                                                           aCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATg                                           <  1:1005904/150‑1 (MQ=255)
                                                                                                              ccGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTg                                        <  2:807937/150‑1 (MQ=255)
                                                                                                                gttaaCCTGATGTCGCTGGCAGCGGTCGCGCAATGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAAGTTGCACGCGTATTCACTGAGCATCAGCGAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGAc                                      <  2:191688/147‑1 (MQ=255)
                                                                                                                     ccTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                 <  2:1016017/150‑1 (MQ=255)
                                                                                                                        gATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAt                              <  2:197479/150‑1 (MQ=255)
                                                                                                                         aTGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCCGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAtt                             >  1:1010115/1‑150 (MQ=255)
                                                                                                                            tCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGAc                          <  1:392408/150‑1 (MQ=255)
                                                                                                                             cGCTGGCACCGGTCGCGCAAAGCGCGCGGGTTGCCAGTGAGAAAGTTGATATCTGCGATCGGGCAGAACAGGCACGCGGATTAACATCGCCTAAGCCAGCCAGTGGCATGCAACTGGCAGCCCAAACCCCTGTTGACCCGTCTATTGAcc                         >  2:871947/1‑149 (MQ=255)
                                                                                                                             cGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACg                         <  2:223589/150‑1 (MQ=255)
                                                                                                                                    aCCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGc                  <  1:6650/150‑1 (MQ=255)
                                                                                                                                       ggTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGc               >  2:561478/1‑150 (MQ=255)
                                                                                                                                       ggTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGGTGATATCTGCGATCGGGCAGAACAGGCACGCGTATTGACTGAGCATCAGCCAGCCACTGACAAGCATCTGTCAGACCAACTCCATATTGAACGTTCTATTCACGGCCCCGAAGc               >  2:873524/1‑150 (MQ=255)
                                                                                                                                        gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCg              >  1:1536376/1‑150 (MQ=255)
                                                                                                                                        gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAAGCATCTGGCATCCGAAAGCCATGTTCACCGTTCTATTGACGGCCCGGCAGCg              >  1:1506457/1‑150 (MQ=255)
                                                                                                                                        gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACAGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCg              <  2:746425/150‑1 (MQ=255)
                                                                                                                                        caaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCg              >  2:1572427/6‑150 (MQ=255)
                                                                                                                                          aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCCCTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt            >  2:1332255/4‑150 (MQ=255)
                                                                                                                                          aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCTTCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt            >  2:1579036/4‑150 (MQ=255)
                                                                                                                                           aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGCCCGTTCTATTGACGGCCCGGCAGCGttt           >  2:1162762/3‑150 (MQ=255)
                                                                                                                                           aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGttt           >  1:491257/3‑150 (MQ=255)
                                                                                                                                             gcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAt         >  1:27041/1‑150 (MQ=255)
                                                                                                                                             gcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAt         <  2:928007/150‑1 (MQ=255)
                                                                                                                                                aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGTCCCGGCAGCGTTTATTGa      >  1:996287/1‑150 (MQ=255)
                                                                                                                                                aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGa      >  1:1477379/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAACC  >  NC_000913/3972524‑3972821

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