Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA39.5% 92.0 / 50.0 77intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (41/5);  new base A (0/30);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 2.51e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5009‑5307
                                                                                                                                                          |                                                                                                                                                      
aaTCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGg‑ag                                                                                                                                                        <  2:416890/151‑1 (MQ=255)
    atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAAt                                                                                                                                                     >  1:950295/1‑150 (MQ=255)
     tcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATa                                                                                                                                                    >  2:600316/1‑150 (MQ=255)
       aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTACCTTAGAGATTAGGATTGCGG‑AGAATaa                                                                                                                                                   <  2:709201/149‑1 (MQ=255)
        gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATaac‑a                                                                                                                                                >  1:1103269/1‑150 (MQ=255)
           aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑ccg                                                                                                                                           <  2:435664/151‑1 (MQ=255)
           aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑cc                                                                                                                                            >  2:1219026/1‑150 (MQ=255)
            aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑ccgc                                                                                                                                          >  1:143988/1‑151 (MQ=255)
            aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑ccgc                                                                                                                                          >  1:936209/1‑151 (MQ=255)
                 aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑ccgcc                                                                                                                                         <  1:530734/147‑1 (MQ=255)
                    cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTtctc                                                                                                                                   >  2:750092/1‑150 (MQ=255)
                    cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTtctc                                                                                                                                   >  2:562460/1‑150 (MQ=255)
                       tATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATc                                                                                                                                >  1:135554/1‑150 (MQ=255)
                           tcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGt                                                                                                                            >  2:1095974/1‑150 (MQ=255)
                              aTCAGGCCGGGTTTGTTTTTATGCACCCGGTTTTTTTTAAAAAAAAATTATGGAA‑AAAAATAACAGGAAAAAA‑GGAAAAATTCTC‑AAAAAATGCGGTAACTTAGATATTAGGTTTGCGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAAt                                                                                                                         <  2:781090/150‑1 (MQ=255)
                               tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGAGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAAtct                                                                                                                       >  1:923392/1‑151 (MQ=255)
                                cAGGCCGTTTTTTTTTTTACCCGGCCCGTTTTTTTAAGAAAAAAATTGGGGAA‑AAAAATGACAGGAAAAAA‑GAAAAATTTCAA‑AATAAAGGCTTTAACTTAGAAATTAGGATGGGAG‑AAAATAAC‑AA‑CCGCCGTTCTCATCGAGTAAtctc                                                                                                                      <  1:473426/151‑1 (MQ=255)
                                cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAAtctc                                                                                                                      >  2:519731/1‑151 (MQ=255)
                                 aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTa                                                                                                                           >  1:592917/1‑145 (MQ=255)
                                  ggCGGGTTTTGTTTTTAGCCACCCCGGTTTTTTTAGAAAAAAATTAGGGAA‑AAAAATGA‑AGGAAAAAAGGGAAAAATCCAA‑AATAAATGGGTTAATTTAGAGATAGGGATTGCGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTa                                                                                                                           <  2:592934/144‑1 (MQ=255)
                                  ggCCGGTTTTGCTTTTATGCAGCCCGTTTTTTTTAAAAAAAAATTAGGGAA‑AAAAAGGAAAGGAAAAAA‑GGAAAAATTCTC‑AAAAAATGCGGTAACTAAGAGATTAGGATTGGGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAAtct                                                                                                                       <  2:407603/148‑1 (MQ=255)
                                     gggttttGTTTTTATGCAGCCCGGTTTTTTTATGAAAAAATTAGGAAA‑AAAAATGACGGGAAAAAA‑GAAAAAATTCTA‑AAAAAGGGGGTTATCTTAAGAATGGGGATGGGGA‑AGAATAAA‑AA‑CCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  1:426974/145‑1 (MQ=255)
                                     gggttttGCTTTTAGGCACCCCGTTTTTTTTATAAAAAATTTAGGGAA‑AAAAATGAAAGGAAAAAA‑GGAAAAATTCTC‑AATAAATGCGGTAATTTAGAAATTAGGATGGCGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 <  1:1087142/147‑1 (MQ=255)
                                     cGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 >  2:82040/1‑151 (MQ=255)
                                       ggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               >  2:195721/1‑151 (MQ=255)
                                         tttttttttATGCACCCCGTTTTTTTTAGAAAAAAATTAGGAAA‑AAAAATGACAGAAAAAAA‑GGAAAAATTCAA‑AAAAAATGGGATAATTTAGAAATTGGGTTGGGGG‑AAAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                               <  2:526032/144‑1 (MQ=255)
                                         tttGCTTTTATGCAGCCCGGCTTTTTTATGAAGCAATTATGGAG‑AAAACTGACAGGGAAAAC‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CTGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             >  1:159802/1‑151 (MQ=255)
                                          ttGTTTTTTTGCAGCCGGATTTTTTTAGGAAAAAATATTGGAA‑AAAAAGGATGGGAAAAAA‑GGA‑AAATTTTCAAATAAAGCCGTAAACTAAGAGGTTAGGAGGGGGG‑AAAATAAC‑CG‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             <  2:344631/150‑1 (MQ=255)
                                          ttGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                             >  2:84216/1‑150 (MQ=255)
                                            gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                          >  2:236152/1‑151 (MQ=255)
                                               tttATGCAGCCCGTTTTTTTTAAAAAAAAATTATGGAA‑AAAAATGAAAGGAAAAAA‑GAAAAAATCCAA‑AAAAAATGCGTTACTTAAGATATAAGGTTTGGGG‑AAAAAAAC‑AG‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                           <  1:300369/147‑1 (MQ=255)
                                               tttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                        >  1:1170052/1‑150 (MQ=255)
                                                 tATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      >  2:939796/1‑150 (MQ=255)
                                                  aggccgCCGGGTTTTTTAAAAAAAAAATGATGGAA‑AAAAAGAAAAGGAAAAAA‑GAAAAATTTCCA‑AATAAATGCGGTAATTTAGAAATTAGGTTTGCGG‑AAAAAACC‑AA‑CCGCCGTTCTCATCTAGTAATCTCCGGATATCGACCCATAAc                                                                                                    <  1:74861/146‑1 (MQ=255)
                                                   gccgccccGGTTTTTTTATGAAAAAATTATGGAA‑AAAAATAACGGGAAAAAA‑GGAGAATTCCTA‑AAAAACTGGGGAAATTAAGGGTGGAGATGGGGGG‑AGAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                   <  2:592356/146‑1 (MQ=255)
                                                        cccgttttttttATAAAAAATTTAGGGAA‑AAAAATGACAGGAAAAAA‑GGAGAAATTCTA‑AATAAATGGGTAAATTTAAAAATAAGGTTGGGGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                <  2:1120901/143‑1 (MQ=255)
                                                        cccGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                               >  2:74219/1‑150 (MQ=255)
                                                         ccctgttttttAAAAAAAAAATAAGGAAAAAAAAGGAAAGAAAAAAA‑GAAAAAATTCTA‑AAAAAATGCGATAATTTAGAAATTAGGATTGGGA‑AAAAAAAA‑AG‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                               <  1:658213/145‑1 (MQ=255)
                                                           ggtttttttAAAAAAAAATTAGGAAA‑AAAAAGGACAGGAAAAAG‑GAAAAAATTCAA‑AAAAAATGCGGTAATTTAGAGATTGGGTTGGCGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                           <  1:1199939/148‑1 (MQ=255)
                                                                           aattaatgaaaaaaaaGGACAGAAAAAAA‑GGAAAAATTCTA‑AAAAAAGGGGTTAATTAAGAGATGAGGATGGGGG‑A‑AAAAAC‑AACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                            <  1:398745/140‑1 (MQ=255)
                                                                           aaTTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  1:176675/1‑151 (MQ=255)
                                                                           aaTTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  2:68572/1‑151 (MQ=255)
                                                                           aaTTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                           >  2:494773/1‑151 (MQ=255)
                                                                             ttATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         >  1:1023394/1‑151 (MQ=255)
                                                                                  gag‑aAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                     >  1:850862/1‑150 (MQ=255)
                                                                                  gaaaaaaaaTGAAAGGAAAAAG‑GAAAAAATTCAA‑AAAAAAGGCGGAATTTAAGAAATTA‑GATTGGGGAAAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                         <  1:452393/143‑1 (MQ=255)
                                                                                  gaaaaaaaaGGAAAGGAAAAAA‑GGAAAAATTCTA‑AAAAAAGGGGGTAACTTAGATTTTAGGATTGCGA‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        <  1:303507/144‑1 (MQ=255)
                                                                                        aaaagagagaaaaaaa‑GGAAAAATTCTA‑AAAAAAGGGGGTAATTTAGAGATTAGGTTTGGGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               <  2:1061715/140‑1 (MQ=255)
                                                                                        aaaTGAGGGGAAAAAG‑GGAAAATTTCAA‑AAAAAATGCGTAATCTAAGAAATTAGTATGGGGA‑AAAAAAAC‑AG‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               <  2:234417/151‑1 (MQ=255)
                                                                                         aaggaaaggaaaaaa‑ggaaaaaTTCAA‑AAAAAATGCGTTAACTAAGAGATTAGGATTGGGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                              <  2:686359/145‑1 (MQ=255)
                                                                                         aaaaaCAGGAAAAAA‑GGAAAAATTCAA‑AAAAATGGCGGAAACTGAGAGTTTAGGATGGGGA‑AAAAAAAA‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATgaa                                                              <  2:1103335/147‑3 (MQ=255)
                                                                                          ataaaaGGAAAAAA‑GAAAAAATTCTA‑AAAAAAGGCGGTAATTTAGAGATTAGGATGGGGA‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                             <  2:1041196/146‑1 (MQ=255)
                                                                                                ggaaaaaa‑GGAAAAATTCAA‑AAAAAAGGCGATAACTAAGAGATTAGGATGGCGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                         <  2:809433/146‑1 (MQ=255)
                                                                                                  gAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     >  1:664692/1‑151 (MQ=255)
                                                                                                  aaaaaa‑GGAAAAATTCTA‑AAAAAAGGCGGTAACTTAGAGATAAGGATGGGGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                     <  2:143994/150‑1 (MQ=255)
                                                                                                      aa‑GGAAAAATTCTA‑AAAAAATGCGGTAATTAAGAGATTGGGATTGGGG‑AAAAAACC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                 <  2:472912/151‑1 (MQ=255)
                                                                                                           aaaaaTTCTA‑AAAAAAGGGGTTAACTTAGAGATTAGGTTTGGGG‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                               <  1:609469/147‑1 (MQ=255)
                                                                                                            gaAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                            >  1:818467/1‑151 (MQ=255)
                                                                                                                    taaaaaaaGGCGGAAATTAAGAAATTGGGATGGGGA‑AAAAAAAA‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                        <  1:689757/143‑1 (MQ=255)
                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                  >  2:405747/1‑151 (MQ=255)
                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                  >  2:608897/1‑151 (MQ=255)
                                                                                                                       aTAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   >  1:1171886/1‑150 (MQ=255)
                                                                                                                         aaaTGCGGAAAGTAAGA‑AATAGGAATGGGG‑AAAAAAACAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                <  1:936058/151‑1 (MQ=255)
                                                                                                                          aaTGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               >  1:640778/1‑151 (MQ=255)
                                                                                                                             gggaaaattaaaggaTGAGTATGGCGG‑AAAAAAAC‑AG‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                               <  1:737209/135‑1 (MQ=255)
                                                                                                                             gcttaaagtaagagATT‑GGAATGGGG‑AAAAAAAA‑AAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            <  2:1208732/141‑1 (MQ=255)
                                                                                                                             gCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            >  2:1218336/1‑151 (MQ=255)
                                                                                                                               ggTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGc                          >  2:20646/1‑151 (MQ=255)
                                                                                                                                  aaCTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                       >  2:585722/1‑151 (MQ=255)
                                                                                                                                   aCTTAGTGATTAGGATTGCGG‑AGAATAAC‑AA‑CTGCCGTTCTCATCGAGTAATTTCCGGATATCGACCCATAACGGGCAATGATAAAACGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCGTGGTTCTGGTCGCTc                        >  1:14245/1‑149 (MQ=255)
                                                                                                                                     ttAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    >  1:343785/1‑151 (MQ=255)
                                                                                                                                     taagatATTGGGATGGGGA‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                           <  2:1063203/142‑1 (MQ=255)
                                                                                                                                             tAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac            >  2:461839/1‑151 (MQ=255)
                                                                                                                                                   tGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg      >  2:633790/1‑151 (MQ=255)
                                                                                                                                                    gggg‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTGTCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc     <  2:762209/149‑1 (MQ=255)
                                                                                                                                                      gg‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg   >  1:31546/1‑151 (MQ=255)
                                                                                                                                                      gg‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTATGGTCGCTCCCATGGCAGCACAGGCTGCgg   >  2:690120/1‑151 (MQ=255)
                                                                                                                                                      gg‑AAAAAAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg   <  2:224328/151‑1 (MQ=255)
                                                                                                                                                         agaATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  1:1127758/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
AATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAG‑AAAAATGACAGGGAAAAA‑GGAGAAATTCTC‑AATAAATGCGGTAACTTAGAGATTAGGATTGCGG‑AGAATAAC‑AA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5009‑5307

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