Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA47.9% 39.6 / 51.7 72intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (31/6);  new base A (0/34);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 3.19e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAA  >  NC_000913/5012‑5309
                                                                                                                                                       |                                                                                                                                                       
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                      <  1:72837/151‑1 (MQ=255)
     gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                                 >  2:1372669/1‑151 (MQ=255)
     gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                                 >  2:663343/1‑151 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                              <  2:1510716/151‑1 (MQ=255)
         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                             >  1:151401/1‑150 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  2:591744/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  2:560/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  1:1156474/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  2:1242141/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  2:825871/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  2:818307/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGa                                                                                                                                >  1:209903/1‑151 (MQ=255)
                       cacAATCAGGCGGGTTTTGTTTTTATGCACCCGTGTTTTTTTAGAAAAAAATTAGGAAAAAAAATAACAGAAAAAAAGAAAAATTCCTAAATAATGGCGGTAACTTAGAAATTAGGTTGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGt                                                                                                                              <  2:1487474/149‑1 (MQ=255)
                         cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGCGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGaaa                                                                                                                            >  1:649962/1‑151 (MQ=255)
                         cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTaa                                                                                                                            >  2:17244/1‑151 (MQ=255)
                         aaaaggggCGGGTTTTGTTTTTATGCACCCGGGTTTTTTTAAAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCAAAAAAAAGGGGGTAATTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTaa                                                                                                                            <  1:1257629/145‑1 (MQ=255)
                           aTCAGGCGGGGTTTGCTTTTAGCCAGCCGGGTTTTTTTATAAAAAAATTAGGAAAAAAAAGGACAGGAAAAAAGAAAAAATTCTCAAAAAATGCGGTAATTAAGAAATTAGGTTTGCGGAAAAAACCAACCGCCGTTCTCATCGAGTAAtc                                                                                                                          <  2:938370/151‑1 (MQ=255)
                              aaggcGGGTTTTCTTTTAAGGCGCCCGGGTTTTTTAAAAAAAAATTTAGGGAAAAAAATGACAGGAAAAAAGAAAAAATCATAAAAAAAGGCGATAATTTAGAAATTAGGATGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAAtctc                                                                                                                        <  2:125298/147‑1 (MQ=255)
                               ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                      >  1:301325/1‑151 (MQ=255)
                                 cgggttttGTTTTTATGCAGCCCGGTTTTTTTAAGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGa                                                                                                                    <  1:486394/146‑1 (MQ=255)
                                   gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                                  >  2:597571/1‑151 (MQ=255)
                                         gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                             >  2:353991/1‑150 (MQ=255)
                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            >  2:1097406/1‑150 (MQ=255)
                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            >  2:1522810/1‑150 (MQ=255)
                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            >  2:718038/1‑150 (MQ=255)
                                          cTTTTACACACCCCGGTTTTTTTAAAAAAAAATTAGGGAAAAAAAGGAGGGAAAAAAGGAAAAAATTCTAAATAATGCCAGTAATTGAAGGCTGAGAAGCGAGAAAAAAAACCACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            <  1:1014351/150‑1 (MQ=255)
                                            tttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                         >  1:37549/1‑151 (MQ=255)
                                                  gcgCCCAGGTTTTTTAAAAAAAAATTTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCACAAAAAAGGGTTTAATTAAGAAATTGGGATTGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                      <  2:1203738/147‑1 (MQ=255)
                                                  gcacccGGGTTTTTTTAGAAAAAAATAAGGAAAAAAAAGAAGAGGAAAAAAGAAAAAATTCTAAAAAAAGGCGTTAATTTAGAGATTGGGATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                      <  1:325931/146‑1 (MQ=255)
                                                   cgCCCGTTTTTTTTAAAAAAAAATTAGGGAAAAAAAAGACAGGAAAAAAGAAAAAATTCTCAAAAAATGCGGTAATTAAAAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                  <  1:1386832/150‑1 (MQ=255)
                                                         gCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                            >  1:170560/1‑151 (MQ=255)
                                                           tttttttaaaaaaaaTTGAGGAAAAAAATGGCCGGGAAAAAGGGAAAAATTCCAAAAAAGGGCGGTAACTAAGAGATTAGGTTTGGGAAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                            <  1:804948/138‑1 (MQ=255)
                                                               ttaaaaaaaaaTTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGTAATTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:1139316/143‑1 (MQ=255)
                                                               ttaaaaaaaaaTTAGGGAAAAAAAGGAGAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAAATAAGAGATTAGGATGGGGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:811818/143‑1 (MQ=255)
                                                               ttATGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAAATGCGGTAATTTAGAGATTAGTTTTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  2:196346/151‑1 (MQ=255)
                                                               ttATAAAAAAATTAGGAAAAAAAATAACGGGAAAAAAGGAAAAATTCCAAAAAAAGGCGGTAATTAAGAGATTAGGATTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:482807/151‑1 (MQ=255)
                                                               ttATAAAAAAATAAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAATTGGGGAAATTTAGAGTTTGGGGTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:1035038/151‑1 (MQ=255)
                                                               taataaaaaaaTTAGGAAAAAAAAGGACAGGAAAAAAGGAAAAATCCTAAAAAAATGCGGTAATTTAGAAATTAGGATTGGGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  2:22484/143‑1 (MQ=255)
                                                                    aagaaATTAGGAAAAAAAATGACAGGAAAAAAGGAAAATTTTTAAAAAAATGCGGTAACTTAAAGATTAGGTTGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                                 <  1:591671/151‑1 (MQ=255)
                                                                    aaaaaaTTATGGAAAAAAATGACGGGAAAAAAGGAAAAATTCTAAAAAAAGGGGGAAATTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                                 <  1:1115534/148‑1 (MQ=255)
                                                                     agaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                                 >  2:20752/1‑150 (MQ=255)
                                                                      aaaaTTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCACAAAAAAGGCGGTAACTAAAAAATAAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  1:451238/150‑1 (MQ=255)
                                                                      aaaaTTATGGAAAAAAAGGACAGAGAAAAAGGAAAAATTCTAAAAAAGTGGGGTAATTTAGAGATTAGGGTTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  2:258188/150‑1 (MQ=255)
                                                                      aaaaTTATGAAAAAAAATGAAAGGGAAAAAGAAGAAATTCTAAATAACGGCATTACCTTAGAGTTTGGGATTGAGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  2:908401/150‑1 (MQ=255)
                                                                      aaaaTTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAATTGGGGTAATTTAGAAATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  1:881228/150‑1 (MQ=255)
                                                                      aaaaTTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCAAAATAAATGCGTTAACTTAGAAATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  2:372067/150‑1 (MQ=255)
                                                                               ggAAAAAAATGACGGAGAAAAAGGAAAAATTCAAAAAAAATGCGGTAACTTAGAGATAAGGTTGGGGAAAAAAAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                        <  2:167030/150‑1 (MQ=255)
                                                                                     aaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                  >  1:1206320/1‑150 (MQ=255)
                                                                                           aagggaaaaaaGGAAAATTTCAAAAAAAAGGCGGTAATTAAGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                             <  1:284556/144‑1 (MQ=255)
                                                                                           aaaggaaaaaaggaaaaaTTCAAAAAAAATGCGGTAACTTAGAGATTGGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                            <  1:1163236/149‑1 (MQ=255)
                                                                                              gggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                          >  2:644209/1‑150 (MQ=255)
                                                                                                  aaaaGGAAAATTTCTAAAAAAGTGCGTTAATTAAGAGTTTGGGATGGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  1:826331/151‑1 (MQ=255)
                                                                                                  aaaaGGAAAAATTCACAATAAATGGGGTAACTTAAAAATTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  1:927702/151‑1 (MQ=255)
                                                                                                  aaaaGGAAAAATTAAAAAAAAAGGGGGAAATTAAGAAATTGGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  1:1409004/151‑1 (MQ=255)
                                                                                                  aaaaGGAAAAATCCACAATAAATGCGTTAACGAGGAGCTTAGTATGGGGGACAATAGCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCGGTGAAAAAGATGCAATCTATCGTAc                                                     <  2:1518658/151‑1 (MQ=255)
                                                                                                  aaaaGAAAAAATTCTAAAAAAATGGGGTACCTAAAAAATGAGGATTGCGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  2:383505/151‑1 (MQ=255)
                                                                                                          aaaaaTTCACAAAAAGGGCGGAAACGAAAAGATGAGGATGGGGAAGAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                               <  2:1419757/149‑1 (MQ=255)
                                                                                                                tctaAATAAGTGGTGTAACTGAAAAATGAGGGTTCCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                         <  2:573970/151‑1 (MQ=255)
                                                                                                                tcaaaataagggggttAACGTAGATATTAGTAGGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                         <  2:770591/136‑1 (MQ=255)
                                                                                                                   cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                      >  1:1014293/1‑151 (MQ=255)
                                                                                                                          tgcgttacatgagagATAAGGATTGAGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                               <  1:986514/140‑1 (MQ=255)
                                                                                                                                aaCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTc                          >  1:1421319/1‑150 (MQ=255)
                                                                                                                                 aCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                         >  1:310448/1‑150 (MQ=255)
                                                                                                                                 aCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                         >  2:1608469/1‑150 (MQ=255)
                                                                                                                                  cTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                       >  1:710683/1‑151 (MQ=255)
                                                                                                                                  cTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                       >  1:1034741/1‑151 (MQ=255)
                                                                                                                                       agaTTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg                    <  1:494616/149‑1 (MQ=255)
                                                                                                                                       aaattgggattgggAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                   <  1:1305457/148‑1 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag                 >  2:735555/1‑151 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag                 >  2:1262765/1‑151 (MQ=255)
                                                                                                                                           tgggttgggggAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac              <  1:597380/142‑1 (MQ=255)
                                                                                                                                           tAGGATGGCGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac              <  1:415356/151‑1 (MQ=255)
                                                                                                                                                       gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  2:1379227/1‑151 (MQ=255)
                                                                                                                                                       gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  2:1408647/1‑151 (MQ=255)
                                                                                                                                                       gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  1:436156/1‑151 (MQ=255)
                                                                                                                                                       gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  1:355099/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAA  >  NC_000913/5012‑5309

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