Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,468 A→T 5.3% M20L (ATG→TTG)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,4680AT 5.3% 237.6 / 2.8 75M20L (ATG→TTG) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base T (2/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

CGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGAT  >  NC_000913/869321‑869616
                                                                                                                                                   |                                                                                                                                                    
cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                                                                                                                    <  2:185516/150‑1 (MQ=255)
cGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATg                                                                                                                                                    <  2:793044/150‑1 (MQ=255)
    gCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCAt                                                                                                                                                >  2:1731393/1‑150 (MQ=255)
      aaTGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCt                                                                                                                                              >  1:1904567/1‑150 (MQ=255)
             cGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGGAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCCTCTAGTGcct                                                                                                                                       >  2:948503/1‑148 (MQ=255)
              gggCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCAtt                                                                                                                                      >  2:444179/1‑150 (MQ=255)
               ggCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTc                                                                                                                                     >  1:980546/1‑150 (MQ=255)
                gCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCg                                                                                                                                    <  1:362920/150‑1 (MQ=255)
                     ggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGcc                                                                                                                               >  2:1761381/1‑150 (MQ=255)
                     ggCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGcc                                                                                                                               >  1:1678499/1‑150 (MQ=255)
                      gCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCa                                                                                                                              >  2:185055/1‑150 (MQ=255)
                         ttACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGAAAGAGCTGTACACCGTAGTGGGTTAGTGGCACTGGGCATTGCGACAGCGATGATGGCATCTTGTGAATTCGCTGCCAAAg                                                                                                                           >  1:877390/1‑150 (MQ=255)
                           aCGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGCATAAAATGGCAAGCGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGAt                                                                                                                         >  2:883678/1‑150 (MQ=255)
                            cGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGAtg                                                                                                                        <  2:1904567/150‑1 (MQ=255)
                              tCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCTGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTTGCGCTGGGCATTGCGACAGCGTTGATGGCATATTGTGCATTAGCTGCCAAAGAtttg                                                                                                                      >  2:813558/1‑150 (MQ=255)
                               cGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATgtgg                                                                                                                     <  1:722174/150‑1 (MQ=255)
                                           aTCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGTTGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAg                                                                                                         >  2:1425309/1‑150 (MQ=255)
                                           aTCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGTTGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAg                                                                                                         >  2:1432371/1‑150 (MQ=255)
                                             ttatctACGCATGAATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGcc                                                                                                                               <  2:636415/121‑1 (MQ=255)
                                             ttatctACGCATGAATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGcc                                                                                                                               >  1:636415/6‑126 (MQ=255)
                                             cAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGa                                                                                                       >  2:23098/1‑150 (MQ=255)
                                                  tACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGaa                                                                                                  >  1:597548/1‑150 (MQ=255)
                                                   aCGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAAt                                                                                                 <  1:1247641/150‑1 (MQ=255)
                                                    cGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAAtt                                                                                                >  1:933123/1‑150 (MQ=255)
                                                      cATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTc                                                                                              >  1:1805844/1‑150 (MQ=255)
                                                          cATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTcacca                                                                                          <  1:112070/150‑1 (MQ=255)
                                                                tAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCg                                                                                    <  2:389045/150‑1 (MQ=255)
                                                                  gATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGAt                                                                                  >  1:1577744/1‑150 (MQ=255)
                                                                   aTAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATc                                                                                 <  1:444179/150‑1 (MQ=255)
                                                                      aCATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGt                                                                              <  2:898718/150‑1 (MQ=255)
                                                                      aCATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGt                                                                              >  1:802745/1‑150 (MQ=255)
                                                                       cATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGGGGCGGTATTATCGAATTTCACCACTCTCGATCCGTa                                                                             <  2:1610526/150‑1 (MQ=255)
                                                                           cAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGAATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGAc                                                                         >  2:1525205/1‑150 (MQ=255)
                                                                              gCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGGTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCa                                                                      <  2:1841070/150‑1 (MQ=255)
                                                                                 gagaATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAAt                                                                   >  1:903434/1‑150 (MQ=255)
                                                                                     aTAAAATGGCAAGAGCTGTACACCGTAGTGGGTAAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGaca                                                               >  1:262845/1‑150 (MQ=255)
                                                                                      tAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGacac                                                              >  1:1914662/1‑150 (MQ=255)
                                                                                         aaTGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGtt                                                           >  2:864366/1‑150 (MQ=255)
                                                                                         aaTGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGtt                                                           >  1:663613/1‑150 (MQ=255)
                                                                                         aaTGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGtt                                                           >  2:307640/1‑150 (MQ=255)
                                                                                          aTGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTa                                                          <  1:23098/150‑1 (MQ=255)
                                                                                              cAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTAtctc                                                      <  1:206636/150‑1 (MQ=255)
                                                                                                   gCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGcc                                                 >  2:126604/1‑150 (MQ=255)
                                                                                                    cTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCg                                                <  1:1761381/150‑1 (MQ=255)
                                                                                                    cTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCg                                                >  1:415219/1‑150 (MQ=255)
                                                                                                           ccGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGaa                                         >  1:1606493/1‑150 (MQ=255)
                                                                                                                 tGGGTTAGTGGCGGTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGTTCGAATTTCTCGACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGtt                                   <  1:1671276/150‑1 (MQ=255)
                                                                                                                  gggTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGttt                                  >  1:493758/1‑150 (MQ=255)
                                                                                                                  gggTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGttt                                  <  1:1731393/150‑1 (MQ=255)
                                                                                                                  gggTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGttt                                  <  2:933123/150‑1 (MQ=255)
                                                                                                                  gggTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGttt                                  <  2:1805844/150‑1 (MQ=255)
                                                                                                                  gggTTAGTGGCGCGTGGCATTGCGACAGCTTTGATGGCATCTTGTTCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATGGAATTTCACCACGCTCTATCCGTATGACGAAAATGACACGTTATCTCAGGCCGTAGCGAAATCGttt                                  <  2:597548/150‑1 (MQ=255)
                                                                                                                   ggTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGtttt                                 <  1:126604/150‑1 (MQ=255)
                                                                                                                           gcgcTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAgggg                         >  1:1656043/1‑150 (MQ=255)
                                                                                                                             gcTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCt                       >  1:462803/1‑150 (MQ=255)
                                                                                                                               tGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGt                     <  1:669345/150‑1 (MQ=255)
                                                                                                                                  gCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCg                  >  1:1081800/1‑150 (MQ=255)
                                                                                                                                   cATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCgg                 >  1:975034/1‑150 (MQ=255)
                                                                                                                                   aaTTGCGACATCGTTGATGTCATCTTGTGCATTCTCTGCCAAAGATGTGGTGGTGGCGGTAGGTTCTAATTTCACGACGCTCGATTCGTATTACGCAAATGACAGTTTATCTCAGGCCTTAGCGAAATCGTTTTTGCAGGGGCTGTTCgg                 <  2:262845/149‑1 (MQ=255)
                                                                                                                                     ttGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTc               <  1:307640/150‑1 (MQ=255)
                                                                                                                                      tGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCt              >  2:1838178/1‑150 (MQ=255)
                                                                                                                                      tGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCt              <  2:493758/150‑1 (MQ=255)
                                                                                                                                       gCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTg             <  1:1404848/150‑1 (MQ=255)
                                                                                                                                       gCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTg             <  2:1577744/150‑1 (MQ=255)
                                                                                                                                        cGACAGCGTTGTTGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTGACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTgg            <  1:1432371/150‑1 (MQ=255)
                                                                                                                                        cGACAGCGTTGTTGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTgg            <  1:1425309/150‑1 (MQ=255)
                                                                                                                                         gACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGa           <  2:663613/150‑1 (MQ=255)
                                                                                                                                         gACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGa           <  2:1678499/150‑1 (MQ=255)
                                                                                                                                         gACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGAGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGACTGGa           >  2:1596669/1‑150 (MQ=255)
                                                                                                                                          aCAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGAt          >  2:836138/1‑150 (MQ=255)
                                                                                                                                           cAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATa         <  2:1606493/150‑1 (MQ=255)
                                                                                                                                             gCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATaaa       <  1:813558/150‑1 (MQ=255)
                                                                                                                                             gCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATaaa       <  2:980546/150‑1 (MQ=255)
                                                                                                                                              cgatcaTGGCATCTTGTGCATTCGCTCCCAAAGATGTCGTGGTAGCGGTAGGATCCAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCGGTAGCGCAATCGTTTTACCACGGGCTGTTCCGTCTGGAGAAag      >  2:543151/6‑150 (MQ=255)
                                                                                                                                               gTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAaga     >  2:1763584/1‑150 (MQ=255)
                                                                                                                                                  gATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGAt  >  2:1878593/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTTGCAATGCGTCGGGCCGGGGCATTACGTCGCACAACCACAATCAGAATACGCATTCATGCGTAGATAACATTCAGGCGGAGAATAAAATGGCAAGAGCTGTACACCGTAGTGGGTTAGTGGCGCTGGGCATTGCGACAGCGTTGATGGCATCTTGTGCATTCGCTGCCAAAGATGTGGTGGTGGCGGTAGGATCGAATTTCACCACGCTCGATCCGTATGACGCAAATGACACGTTATCTCAGGCCGTAGCGAAATCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGAT  >  NC_000913/869321‑869616

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