Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 828,749 A→T 5.2% L78Q (CTG→CAG)  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,7490AT 5.2% 297.3 / 4.9 97L78Q (CTG→CAG) ybhFputative ABC exporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (35/57);  new base T (3/2);  total (38/59)
Fisher's exact test for biased strand distribution p-value = 3.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

GAGACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGTA  >  NC_000913/828600‑828898
                                                                                                                                                     |                                                                                                                                                     
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                                                                                            >  1:882723/1‑150 (MQ=255)
                                                           cgcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGt                                                                                            >  2:1483033/1‑150 (MQ=255)
                                                            gcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGtt                                                                                           <  2:227256/150‑1 (MQ=255)
                                                            gcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGtt                                                                                           <  2:1347062/150‑1 (MQ=255)
                                                              aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                         <  1:687046/150‑1 (MQ=255)
                                                              aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                         <  1:545856/150‑1 (MQ=255)
                                                              aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                         <  1:1369298/150‑1 (MQ=255)
                                                              aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                         <  1:1155704/150‑1 (MQ=255)
                                                                aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGTGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                       <  2:887115/150‑1 (MQ=255)
                                                                aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                       <  2:585225/150‑1 (MQ=255)
                                                                 tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                      <  1:969522/150‑1 (MQ=255)
                                                                 tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                      <  2:1273105/150‑1 (MQ=255)
                                                                 tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                      <  1:92402/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                                                                           >  1:657431/1‑150 (MQ=255)
                                                                            aTTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCaa                                                                           <  2:936564/150‑1 (MQ=255)
                                                                             ttGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAAc                                                                          >  1:279535/1‑150 (MQ=255)
                                                                              tGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACa                                                                         >  2:748473/1‑150 (MQ=255)
                                                                               gAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACAt                                                                        >  1:1607485/1‑150 (MQ=255)
                                                                               gAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACAt                                                                        >  2:1364174/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                                                               >  2:1631213/1‑150 (MQ=255)
                                                                                        cATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAg                                                               <  1:896895/150‑1 (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                                      >  2:430533/1‑150 (MQ=255)
                                                                                                                 ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                      <  2:1363997/150‑1 (MQ=255)
                                                                                                                 ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                      <  1:809881/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)
                                                                                                                                                   gCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATgg    >  1:1281399/1‑150 (MQ=255)
                                                                                                                                                   gCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATgg    <  2:9468/150‑1 (MQ=255)
                                                                                                                                                     aGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGTa  <  1:1475641/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGTA  >  NC_000913/828600‑828898

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: