Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 828,746 T→C 5.2% H79R (CAC→CGC)  ybhF ← putative ABC exporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913828,7460TC 5.2% 296.3 / 4.2 96H79R (CAC→CGC) ybhFputative ABC exporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (35/56);  new base C (3/2);  total (38/58)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

AGAGACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATG  >  NC_000913/828599‑828895
                                                                                                                                                   |                                                                                                                                                     
agagACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGc                                                                                                                                                     <  1:646142/150‑1 (MQ=255)
agagACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGc                                                                                                                                                     >  1:1347062/1‑150 (MQ=255)
 gagaCGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCa                                                                                                                                                    <  1:441449/150‑1 (MQ=255)
   gaCGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAgc                                                                                                                                                  >  1:613702/1‑150 (MQ=255)
      gTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAgcgcg                                                                                                                                               >  1:789235/1‑150 (MQ=255)
       tAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGCGCTgcacgc                                                                                                                                              >  2:1091859/1‑150 (MQ=255)
         aaCTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCcg                                                                                                                                            <  2:1008858/150‑1 (MQ=255)
           cTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCcgtc                                                                                                                                          <  1:867367/150‑1 (MQ=255)
                gcagGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGtttt                                                                                                                                     <  1:652145/150‑1 (MQ=255)
                gcagGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGtttt                                                                                                                                     <  1:1173194/150‑1 (MQ=255)
                gcagGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGCGCTGCACGCcgtgctcgg                                                                                                                                     <  1:395273/150‑7 (MQ=255)
                 cagGCGATCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCAATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGttttt                                                                                                                                    <  1:306533/150‑1 (MQ=255)
                 cagGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGttttt                                                                                                                                    <  2:1241078/150‑1 (MQ=255)
                  agGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTg                                                                                                                                   <  2:760440/150‑1 (MQ=255)
                     cGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATc                                                                                                                                <  1:56859/150‑1 (MQ=255)
                     cGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATc                                                                                                                                >  1:659015/1‑150 (MQ=255)
                       aGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCgg                                                                                                                              <  2:659015/150‑1 (MQ=255)
                        gCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGa                                                                                                                             >  1:279491/1‑150 (MQ=255)
                              gTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAg                                                                                                                       <  2:1243573/150‑1 (MQ=255)
                               tttGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGCGCTGCACGCCGTCGTTTTTGATCGGATCAAAGc                                                                                                                      <  1:1091859/150‑1 (MQ=255)
                                ttGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGTATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGcc                                                                                                                     <  1:1002983/150‑1 (MQ=255)
                                   cTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAAt                                                                                                                  >  1:585225/1‑150 (MQ=255)
                                    ttACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATGACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATc                                                                                                                 <  1:150534/150‑1 (MQ=255)
                                     tACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATAGACGTTCTCCATGACCGTGAGATCTTCATACAGACGAAATTTCTGCGGCAAATAACCGAGCACGGCGTCCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCaaacc                                                                                                                >  1:285677/1‑147 (MQ=255)
                                     tACGAGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCACCGTCGTTTTTGATCGGATCAAAGCCAATCa                                                                                                                >  2:114895/1‑150 (MQ=255)
                                          gCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTg                                                                                                           <  1:756084/150‑1 (MQ=255)
                                          gCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTg                                                                                                           >  1:896950/1‑150 (MQ=255)
                                          gCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTg                                                                                                           <  2:229655/150‑1 (MQ=255)
                                          gCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTg                                                                                                           <  2:981789/150‑1 (MQ=255)
                                                  ggTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCCGATCAAAGCCAATCACCGTGGCACTGcc                                                                                                   >  2:494805/1‑150 (MQ=255)
                                                      aCGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCCAAGCCAATCACCGTGGCACTGCCGCTg                                                                                               >  2:1570910/1‑150 (MQ=255)
                                                      aCGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTg                                                                                               >  2:1576402/1‑150 (MQ=255)
                                                       cGCTGCGCAGATCCGCGTACAGATTGAGGTTCTCCATCCCCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGt                                                                                              >  2:610706/1‑150 (MQ=255)
                                                            ggcAAATCGGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACGAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGGGCCGTCGTTTTTGATGGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGt                                                                                         <  2:333330/149‑1 (MQ=255)
                                                            cgcAAATCCGCGTACCGCATGAGGCTCTCCATCCCCGTGACATCTTCATCCAGACCAAATTGCTGCGGCATATAACCGACCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCgggg                                                                                         <  2:1559988/150‑2 (MQ=255)
                                                            cgcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGt                                                                                         <  2:138225/150‑1 (MQ=255)
                                                            cgcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGt                                                                                         >  2:1483033/1‑150 (MQ=255)
                                                            cgcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGt                                                                                         >  1:882723/1‑150 (MQ=255)
                                                             gcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGtt                                                                                        <  2:1347062/150‑1 (MQ=255)
                                                             gcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGtt                                                                                        <  2:227256/150‑1 (MQ=255)
                                                               aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                      <  1:1155704/150‑1 (MQ=255)
                                                               aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                      <  1:1369298/150‑1 (MQ=255)
                                                               aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                      <  1:687046/150‑1 (MQ=255)
                                                               aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                      <  1:545856/150‑1 (MQ=255)
                                                                 aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGTGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                    <  2:887115/150‑1 (MQ=255)
                                                                 aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                    <  2:585225/150‑1 (MQ=255)
                                                                  tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                   <  2:1273105/150‑1 (MQ=255)
                                                                  tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                   <  1:92402/150‑1 (MQ=255)
                                                                  tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                   <  1:969522/150‑1 (MQ=255)
                                                                  tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                   >  1:319342/1‑150 (MQ=255)
                                                                            gATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCa                                                                         <  2:789235/150‑1 (MQ=255)
                                                                             aTTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCaa                                                                        >  1:829369/1‑150 (MQ=255)
                                                                             aTTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCaa                                                                        <  2:936564/150‑1 (MQ=255)
                                                                             aTTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCaa                                                                        >  1:657431/1‑150 (MQ=255)
                                                                              ttGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAAc                                                                       >  1:279535/1‑150 (MQ=255)
                                                                               tGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACa                                                                      >  2:748473/1‑150 (MQ=255)
                                                                                gAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACAt                                                                     >  2:1364174/1‑150 (MQ=255)
                                                                                gAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACAt                                                                     >  1:1607485/1‑150 (MQ=255)
                                                                                  ggTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATcc                                                                   >  2:1429400/1‑150 (MQ=255)
                                                                                  ggTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATcc                                                                   >  1:281699/1‑150 (MQ=255)
                                                                                   gTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCg                                                                  >  1:9468/1‑150 (MQ=255)
                                                                                     tttcCATCACCGTGAGATCTTCATACAGACCAAATTTCCGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGTATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCa                                                                <  2:1453626/148‑1 (MQ=255)
                                                                                      ctcCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCAt                                                               <  2:1077473/150‑1 (MQ=255)
                                                                                       tcCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATc                                                              <  2:1553295/150‑1 (MQ=255)
                                                                                        ccATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACTGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCa                                                             <  2:195623/150‑1 (MQ=255)
                                                                                         cATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAg                                                            <  1:896895/150‑1 (MQ=255)
                                                                                         cATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAg                                                            >  2:1631213/1‑150 (MQ=255)
                                                                                             aCCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTg                                                        <  2:260049/150‑1 (MQ=255)
                                                                                               cGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGt                                                      <  2:518786/150‑1 (MQ=255)
                                                                                               cGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCACCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAAGATCCGCATCAGCGTGGt                                                      <  1:114895/150‑1 (MQ=255)
                                                                                                gTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGtt                                                     <  2:1280034/150‑1 (MQ=255)
                                                                                                     atcatcATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGAGTGCAGCGCGGCGTCGTTTTTGATTGTTTCAAAGCCAATGAGCGTGTCAATGACTATTTCGGGTTTCAGTAATACCGCCAACTTCGGAATCAGCTTTGTTTTAcc                                                <  2:1456365/150‑1 (MQ=255)
                                                                                                          cATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCAc                                           <  1:216108/150‑1 (MQ=255)
                                                                                                           aTACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCAcc                                          >  1:1048152/1‑150 (MQ=255)
                                                                                                            tACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCg                                         >  1:927627/1‑150 (MQ=255)
                                                                                                            tACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCg                                         >  1:160260/1‑150 (MQ=255)
                                                                                                                  ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                   <  1:809881/150‑1 (MQ=255)
                                                                                                                  ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                   >  2:430533/1‑150 (MQ=255)
                                                                                                                  ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                   <  2:1363997/150‑1 (MQ=255)
                                                                                                                   cAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGcc                                  <  1:169637/150‑1 (MQ=255)
                                                                                                                   cAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGcc                                  >  2:785003/1‑150 (MQ=255)
                                                                                                                     aaTTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGcccc                                <  2:732335/150‑1 (MQ=255)
                                                                                                                      aTTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCcca                               <  2:1137946/150‑1 (MQ=255)
                                                                                                                         tCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACGTCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCccacca                            <  2:923197/150‑1 (MQ=255)
                                                                                                                           tGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTCCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAAc                          <  2:896950/150‑1 (MQ=255)
                                                                                                                                  tataACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTca                   <  1:1631213/150‑1 (MQ=255)
                                                                                                                                  tataACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTca                   <  2:160260/150‑1 (MQ=255)
                                                                                                                                   ataACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCTGGCCCCACCAACCCCGTcac                  >  2:43973/1‑150 (MQ=255)
                                                                                                                                    taACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTcaca                 >  2:643567/1‑150 (MQ=255)
                                                                                                                                     aaCCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACAt                <  2:279491/150‑1 (MQ=255)
                                                                                                                                      aaaaacgaCGGCGCGCTGCACGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATa               >  1:1407933/8‑150 (MQ=255)
                                                                                                                                         gAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAAcc            >  2:1198716/1‑150 (MQ=255)
                                                                                                                                         aacgaCGGCGCGCTGCCCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAAcc            >  1:98766/5‑150 (MQ=255)
                                                                                                                                            cACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGc         <  1:1429400/150‑1 (MQ=255)
                                                                                                                                             aCGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCg        >  1:1001858/1‑150 (MQ=255)
                                                                                                                                                 cGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGaa    <  1:819256/150‑1 (MQ=255)
                                                                                                                                                   tGCAGCGCGCGGTCGTTTTTGTTCGGATGAAAGCCAATCTCCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCTGCGTGGTTTTTCCTGCTCCGTGCGGGCCCACCAACCCCGTCAGATAACCGGCGTGAATg  <  2:1219469/150‑1 (MQ=255)
                                                                                                                                                   tGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATg  <  2:657431/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGAGACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATG  >  NC_000913/828599‑828895

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: