Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,135,026 A→G V449A (GTC→GCC)  pitB ← phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,0260AG96.3% 366.3 / 1.3 107V449A (GTC→GCC) pitBphosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (67/36);  total (70/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-02

Only 100 of 107 total aligned reads displayed.

TCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACGCCAGCC  >  NC_000913/3134877‑3135171
                                                                                                                                                     |                                                                                                                                                 
tCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGAc                                                                                                                                                  >  1:649790/1‑151 (MQ=255)
      gCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCgtgtg                                                                                                                                             >  2:1278854/1‑150 (MQ=255)
        aaTGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGtgttg                                                                                                                                          >  2:194468/1‑151 (MQ=255)
        aaTGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGtgttg                                                                                                                                          <  1:1157610/151‑1 (MQ=255)
         aTGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGtgttgt                                                                                                                                         <  2:1011739/151‑1 (MQ=255)
                        gCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATccc                                                                                                                          >  1:1124899/1‑151 (MQ=255)
                         cccACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATccc                                                                                                                          >  1:1178849/1‑150 (MQ=255)
                         cccACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATccc                                                                                                                          >  1:1178861/1‑150 (MQ=255)
                         cccACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCa                                                                                                                         <  2:1009799/151‑1 (MQ=255)
                           caccaGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAat                                                                                                                       >  2:1573785/1‑151 (MQ=255)
                           caccaGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAat                                                                                                                       >  1:453515/1‑151 (MQ=255)
                            accaGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAat                                                                                                                       >  1:393787/1‑150 (MQ=255)
                                  aagaaaAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTg                                                                                                                <  2:1051071/151‑1 (MQ=255)
                                   agaaaAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTgg                                                                                                               >  2:1577068/1‑151 (MQ=255)
                                      aaaaTTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCa                                                                                                             >  1:750660/1‑150 (MQ=255)
                                       aaaTTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCa                                                                                                             >  1:1062833/1‑149 (MQ=255)
                                           tGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGAcc                                                                                                        <  1:1043206/150‑1 (MQ=255)
                                           tGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCg                                                                                                       >  2:178963/1‑151 (MQ=255)
                                           tGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCg                                                                                                       >  1:1574417/1‑151 (MQ=255)
                                            gccgccGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCg                                                                                                       >  1:468371/1‑150 (MQ=255)
                                            gccgccGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCg                                                                                                       <  1:935133/150‑1 (MQ=255)
                                             ccgccgGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCg                                                                                                       <  1:1083365/149‑1 (MQ=255)
                                             ccgccgGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGAt                                                                                                     >  1:16785/1‑151 (MQ=255)
                                                ccgGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGAt                                                                                                     >  2:706056/1‑148 (MQ=255)
                                                ccgGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGAt                                                                                                     >  2:818660/1‑148 (MQ=255)
                                                ccgGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCATCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGAt                                                                                                     >  2:315635/1‑148 (MQ=255)
                                                 cgGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGac                                                                                                 >  2:1688748/1‑151 (MQ=255)
                                                  ggTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGaca                                                                                                <  1:539769/151‑1 (MQ=255)
                                                  ggTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGac                                                                                                 >  2:1080078/1‑150 (MQ=255)
                                                     aaagtaaaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACt                                                                                              <  1:854845/150‑1 (MQ=255)
                                                     aaagtaaaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTg                                                                                             >  1:1291411/1‑151 (MQ=255)
                                                     aaagtaaaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTg                                                                                             <  2:571303/151‑1 (MQ=255)
                                                      aagtaaaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTgc                                                                                            >  1:1111363/1‑151 (MQ=255)
                                                      aagtaaaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTgc                                                                                            >  2:1057145/1‑151 (MQ=255)
                                                         taaaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTgccgc                                                                                         >  2:171761/1‑151 (MQ=255)
                                                           aaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGccgccg                                                                                       <  1:1357517/151‑1 (MQ=255)
                                                           aaTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGccgccg                                                                                       >  2:489308/1‑151 (MQ=255)
                                                            aTACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTgccgcc                                                                                        >  2:232810/1‑149 (MQ=255)
                                                               cccACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCAt                                                                                   <  1:1417191/151‑1 (MQ=255)
                                                                 cACGCCATCAGGATGCTGGTTACCGTTTTGCGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGCCGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCAttt                                                                                 <  2:87394/151‑1 (MQ=255)
                                                                 cACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCAttt                                                                                 >  2:1546306/1‑151 (MQ=255)
                                                                  aCGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTtg                                                                                >  2:1088155/1‑151 (MQ=255)
                                                                  aCGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTtg                                                                                >  2:653209/1‑151 (MQ=255)
                                                                     ccATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTgc                                                                             <  2:347850/151‑1 (MQ=255)
                                                                     ccATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTgc                                                                             >  1:116331/1‑151 (MQ=255)
                                                                         cAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTgccgcc                                                                         >  1:764785/1‑151 (MQ=255)
                                                                         cAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTgccgc                                                                          <  2:161433/150‑1 (MQ=255)
                                                                            gATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATg                                                                      <  1:216087/151‑1 (MQ=255)
                                                                               gCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGcc                                                                    <  2:1428542/150‑1 (MQ=255)
                                                                               gCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACCGCCGCCGTCATTTGTGCCGCCATGCCt                                                                   <  1:811211/151‑1 (MQ=255)
                                                                                cTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCtt                                                                  >  2:304142/1‑151 (MQ=255)
                                                                                 tGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTg                                                                 >  2:769058/1‑151 (MQ=255)
                                                                                  ggTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTgc                                                                >  2:1373718/1‑151 (MQ=255)
                                                                                    ttACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTgcgc                                                              >  2:1112040/1‑151 (MQ=255)
                                                                                     tACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTgcgc                                                              <  2:1206829/150‑1 (MQ=255)
                                                                                       ccGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATa                                                           <  1:112264/151‑1 (MQ=255)
                                                                                          ttttACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGt                                                        >  2:1627661/1‑151 (MQ=255)
                                                                                          ttttACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGt                                                        >  1:1423614/1‑151 (MQ=255)
                                                                                           tttACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCTTACGTc                                                       >  2:776055/1‑151 (MQ=255)
                                                                                            ttACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTc                                                       >  2:850221/1‑150 (MQ=255)
                                                                                                 cTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCcg                                                 >  1:1692112/1‑151 (MQ=255)
                                                                                                 cTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCcg                                                 >  1:597415/1‑151 (MQ=255)
                                                                                                  tGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCcgc                                                >  2:1178475/1‑151 (MQ=255)
                                                                                                  tGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCcgc                                                >  2:17394/1‑151 (MQ=255)
                                                                                                   gTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCcgcg                                               >  2:641981/1‑151 (MQ=255)
                                                                                                    tAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCgcgc                                              >  1:1637850/1‑151 (MQ=255)
                                                                                                    tAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCgcgc                                              >  1:329353/1‑151 (MQ=255)
                                                                                                     aaCCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGTATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCAGCCATGCCTTGCGCATACGTCATGCCGCGCt                                             >  2:356636/1‑151 (MQ=255)
                                                                                                         cACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTAcc                                         >  1:1656831/1‑151 (MQ=255)
                                                                                                          aCCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTAcc                                         >  2:457667/1‑150 (MQ=255)
                                                                                                           ccgccgTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCa                                        <  2:721131/150‑1 (MQ=255)
                                                                                                           ccgccgTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCa                                        <  1:1279541/150‑1 (MQ=255)
                                                                                                               cgTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAAtctt                                   <  2:649848/151‑1 (MQ=255)
                                                                                                                gTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAAtcttc                                  >  1:571015/1‑151 (MQ=255)
                                                                                                                 tCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAAtcttct                                 >  1:257924/1‑151 (MQ=255)
                                                                                                                        ccaTCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGa                           <  1:354158/150‑1 (MQ=255)
                                                                                                                         caTCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATg                         >  1:263136/1‑151 (MQ=255)
                                                                                                                         caTCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATg                         >  2:133334/1‑151 (MQ=255)
                                                                                                                           tCGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATgg                        <  1:896398/150‑1 (MQ=255)
                                                                                                                           tCGTCCCTGCAACTGCAGACGAGAGGGAGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGt                       >  1:1089135/1‑151 (MQ=255)
                                                                                                                            cGTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGt                       <  2:285744/150‑1 (MQ=255)
                                                                                                                             gTCCCTGCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTc                      >  1:489880/1‑150 (MQ=255)
                                                                                                                                   gCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCt                <  2:198629/150‑1 (MQ=255)
                                                                                                                                   gCAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTa               <  1:1576971/151‑1 (MQ=255)
                                                                                                                                    cAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTac              >  1:1451695/1‑151 (MQ=255)
                                                                                                                                    cAACTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTac              >  2:1423637/1‑151 (MQ=255)
                                                                                                                                     aaCTGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTaca             <  2:782508/151‑1 (MQ=255)
                                                                                                                                        tGCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACacg           >  1:130767/1‑150 (MQ=255)
                                                                                                                                         gCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACacgac         >  2:194004/1‑151 (MQ=255)
                                                                                                                                         gCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACacgac         <  1:1088070/151‑1 (MQ=255)
                                                                                                                                         gCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACacga          <  2:667496/150‑1 (MQ=255)
                                                                                                                                         gCAGACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACacg           <  1:1241872/149‑1 (MQ=255)
                                                                                                                                            gACGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgcc      >  2:410751/1‑151 (MQ=255)
                                                                                                                                             aCGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgcc      >  1:1063155/1‑150 (MQ=255)
                                                                                                                                              cGAGAGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgcc      >  1:964831/1‑149 (MQ=255)
                                                                                                                                               gagaGGGCGTGTGTTGTGGAGACGGGCATTCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgccagc   >  2:1563601/1‑151 (MQ=255)
                                                                                                                                               gagaGGGCGTGTGTTGTGGAGACGGGCATCCCAATATACCTGGCAAGACCGATAGACACTGCCGCCGTTATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgccag    <  1:1690638/150‑1 (MQ=255)
                                                                                                                                               gagaGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgccagc   >  1:191737/1‑151 (MQ=255)
                                                                                                                                               gagaGGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgccagc   >  1:1536809/1‑151 (MQ=255)
                                                                                                                                                agagGGCGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACgccagcc  >  1:400792/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAAATACCCACGCCATCAGGATGCTGGTTACCGTTTTACGCTGTAACCCACCGCCGTCCACCACCATCGTCCCTGCAACTGCAGACGAGAGGACGTGTGTTGTGGAGACGGGCATCCCAATATAACTGGCAAGACCGATAGACACTGCCGCCGTCATTTGTGCCGCCATGCCTTGCGCATACGTCATGCCGCGCTTACCAATCTTCTCACCGATGGTCATCGCTACACGACGCCAGCC  >  NC_000913/3134877‑3135171

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