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

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAA  >  NC_000913/5012‑5309
                                                                                                                                                       |                                                                                                                                                       
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                      <  1:72837/151‑1 (MQ=255)
     gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                                 >  2:1372669/1‑151 (MQ=255)
     gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                                 >  2:663343/1‑151 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑cc                                                                                                                                              <  2:1510716/151‑1 (MQ=255)
         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑ccg                                                                                                                                             >  1:151401/1‑150 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  2:591744/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  2:560/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  1:1156474/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  2:1242141/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  2:825871/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  2:818307/1‑151 (MQ=255)
                     aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGa                                                                                                                                >  1:209903/1‑151 (MQ=255)
                       cacAATCAGGCGGGTTTTGTTTTTATGCACCCGTGTTTTTTTAGAAAAAAA‑TTAGGAAAAAAAATA‑ACAGAAAAAAA‑GA‑AAAATTCCTAAATAATGGCGGTAACTTAGAAATTAGGTTGGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGt                                                                                                                              <  2:1487474/149‑1 (MQ=255)
                         cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGCGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGaaa                                                                                                                            >  1:649962/1‑151 (MQ=255)
                         cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTaa                                                                                                                            >  2:17244/1‑151 (MQ=255)
                         aaaaggggCGGGTTTTGTTTTTATGCACCCGGGTTTTTTTAAAAAAAAA‑TTAGGGAAAAAAATG‑AAAGGAAAAAA‑GG‑AAAAATTCAAAAAAAAGGGGGTAATTAAGAGATTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTaa                                                                                                                            <  1:1257629/145‑1 (MQ=255)
                           aTCAGGCGGGGTTTGCTTTTAGCCAGCCGGGTTTTTTTATAAAAAAA‑TTAGGAAAAAAAAGG‑ACAGGAAAAAA‑GA‑AAAAATTCTCAAAAAATGCGGTAATTAAGAAATTAGGTTTGCGGAAAAAACCAA‑CCGCCGTTCTCATCGAGTAAtc                                                                                                                          <  2:938370/151‑1 (MQ=255)
                              aaggcGGGTTTTCTTTTAAGGCGCCCGGGTTTTTTAAAAAAAAATTTAGGGAAAAAAATG‑ACAGGAAAAAA‑GA‑AAAAATCATAAAAAAAGGCGATAATTTAGAAATTAGGATGGGGGAAAAAAACAC‑CCGCCGTTCTCATCGAGTAAtctc                                                                                                                        <  2:125298/147‑1 (MQ=255)
                               ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                      >  1:301325/1‑151 (MQ=255)
                                 cgggttttGTTTTTATGCAGCCCGGTTTTTTTAAGAAAAAA‑TTATGGAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAAATTCTAAAAAAAGGGGGTAACTTAGAGATTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGa                                                                                                                    <  1:486394/146‑1 (MQ=255)
                                   gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                                  >  2:597571/1‑151 (MQ=255)
                                         gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                             >  2:353991/1‑150 (MQ=255)
                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            >  2:1097406/1‑150 (MQ=255)
                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            >  2:1522810/1‑150 (MQ=255)
                                          cTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            >  2:718038/1‑150 (MQ=255)
                                          cTTTTACACACCCCGGTTTTTTTAAAAAAAAA‑TTAGGGAAAAAAAGG‑A‑GGGAAAAAA‑GGAAAAAATTCTAAATAATGCCAGTAATTGAAGGCTGAGAAGCGAGAAAAAAAACCA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                            <  1:1014351/150‑1 (MQ=255)
                                            tttATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                         >  1:37549/1‑151 (MQ=255)
                                                  gcgCCCAGGTTTTTTAAAAAAAAATTTAGGGAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAAATTCACAAAAAAGGGTTTAATTAAGAAATTGGGATTGCGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                      <  2:1203738/147‑1 (MQ=255)
                                                  gcacccGGGTTTTTTTAGAAAAAA‑ATAAGGAAAAAAAAGAAGAGGAAAAAA‑GA‑AAAAATTCTAAAAAAAGGCGTTAATTTAGAGATTGGGATGGGGAAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                      <  1:325931/146‑1 (MQ=255)
                                                   cgCCCGTTTTTTTTAAAAAAAAA‑TTAGGGAAAAAAAAG‑ACAGGAAAAAA‑GA‑AAAAATTCTCAAAAAATGCGGTAATTAAAAGATTAGGATTGGGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                                  <  1:1386832/150‑1 (MQ=255)
                                                         gCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                            >  1:170560/1‑151 (MQ=255)
                                                           tttttttaaaaaaaa‑TTGAGGAAAAAAATG‑GCCGGGAAAAAGGG‑AAAAATTCCAAAAAAGGGCGGTAACTAAGAGATTAGGTTTGGGAAGAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                            <  1:804948/138‑1 (MQ=255)
                                                               ttaaaaaaaaa‑TTAGGGAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAAATTCTAAAAAAAGGGGGTAATTAAGAGATTAGGATGGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:1139316/143‑1 (MQ=255)
                                                               ttaaaaaaaaa‑TTAGGGAAAAAAAGG‑AGAGGAAAAAA‑GG‑AAAAATTCTAAAAAAATGCGGTAAATAAGAGATTAGGATGGGGGAGAATAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:811818/143‑1 (MQ=255)
                                                               ttATGAAAAAA‑TTATGGAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAAATTCTCAATAAATGCGGTAATTTAGAGATTAGTTTTGCGGAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  2:196346/151‑1 (MQ=255)
                                                               ttATAAAAAAA‑TTAGGAAAAAAAATA‑ACGGGAAAAAA‑GG‑AAAAATTCCAAAAAAAGGCGGTAATTAAGAGATTAGGATTGCGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:482807/151‑1 (MQ=255)
                                                               ttATAAAAAAA‑TAAGGAAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAAATTCTAAATAATTGGGGAAATTTAGAGTTTGGGGTTGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  1:1035038/151‑1 (MQ=255)
                                                               taataaaaaaa‑TTAGGAAAAAAAAGG‑ACAGGAAAAAA‑GG‑AAAAATCCTAAAAAAATGCGGTAATTTAGAAATTAGGATTGGGAAAAAAAAAAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                      <  2:22484/143‑1 (MQ=255)
                                                                    aagaaA‑TTAGGAAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAATTTTTAAAAAAATGCGGTAACTTAAAGATTAGGTTGGCGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                                 <  1:591671/151‑1 (MQ=255)
                                                                    aaaaaa‑TTATGGAAAAAAATG‑ACGGGAAAAAA‑GG‑AAAAATTCTAAAAAAAGGGGGAAATTTAGAGATTAGGATTGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                                 <  1:1115534/148‑1 (MQ=255)
                                                                     agaaA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                                 >  2:20752/1‑150 (MQ=255)
                                                                      aaaa‑TTATGGAAAAAAATG‑ACAGGAAAAAA‑GG‑AAAAATTCACAAAAAAGGCGGTAACTAAAAAATAAGGATTGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  1:451238/150‑1 (MQ=255)
                                                                      aaaa‑TTATGGAAAAAAAGG‑ACAGAGAAAAA‑GG‑AAAAATTCTAAAAAAGTGGGGTAATTTAGAGATTAGGGTTGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  2:258188/150‑1 (MQ=255)
                                                                      aaaa‑TTATGAAAAAAAATG‑AAAGGGAAAAA‑GA‑AGAAATTCTAAATAACGGCATTACCTTAGAGTTTGGGATTGAGAAAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  2:908401/150‑1 (MQ=255)
                                                                      aaaa‑TTAGGGAAAAAAATG‑AAAGGAAAAAA‑GG‑AAAAATTCTAAAAAATTGGGGTAATTTAGAAATTAGGATGGCGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  1:881228/150‑1 (MQ=255)
                                                                      aaaa‑TTAGGGAAAAAAAGG‑ACAGGAAAAAA‑GG‑AAAAATTCAAAATAAATGCGTTAACTTAGAAATTAGGATGGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                               <  2:372067/150‑1 (MQ=255)
                                                                               ggAAAAAAATG‑ACGGAGAAAAA‑GG‑AAAAATTCAAAAAAAATGCGGTAACTTAGAGATAAGG‑TTGGGGAAAAAAAAAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaa                                                                        <  2:167030/150‑1 (MQ=255)
                                                                                     aaaTG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                  >  1:1206320/1‑150 (MQ=255)
                                                                                           aagggaaaaaa‑GG‑AAAATTTCAAAAAAAAGGCGGTAATTAAGAAATTAGGATTGCGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                             <  1:284556/144‑1 (MQ=255)
                                                                                           aaaggaaaaaa‑gg‑aaaaaTTCAAAAAAAATGCGGTAACTTAGAGATTGGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                            <  1:1163236/149‑1 (MQ=255)
                                                                                              gggAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                          >  2:644209/1‑150 (MQ=255)
                                                                                                  aaaa‑GG‑AAAATTTCTAAAAAAGTGCGTTAATTAAGAGTTTGGGATGGGGGAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  1:826331/151‑1 (MQ=255)
                                                                                                  aaaa‑GG‑AAAAATTCACAATAAATGGGGTAACTTAAAAATTAGGATTGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  1:927702/151‑1 (MQ=255)
                                                                                                  aaaa‑GG‑AAAAATTAAAAAAAAAGGGGGAAATTAAGAAATTGGGATGGGGGAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  1:1409004/151‑1 (MQ=255)
                                                                                                  aaaa‑GG‑AAAAATCCACAATAAATGCGTTAACGAGGAGCTTAGTATGGGGGACAATAGCAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCGGTGAAAAAGATGCAATCTATCGTAc                                                     <  2:1518658/151‑1 (MQ=255)
                                                                                                  aaaa‑GA‑AAAAATTCTAAAAAAATGGGGTACCTAAAAAATGAGGATTGCGGAGAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                     <  2:383505/151‑1 (MQ=255)
                                                                                                          aaaaaTTCACAAAAAGGGCGGAAACGAAAAGATGAGGATGGGGAAGAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                               <  2:1419757/149‑1 (MQ=255)
                                                                                                                tctaAATAAGTGGTGTAACTGAAAAATGAGGGTTCCGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                         <  2:573970/151‑1 (MQ=255)
                                                                                                                tcaaaataagggggttAACGTAGATATTAGTAGGGGGAAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                         <  2:770591/136‑1 (MQ=255)
                                                                                                                   cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                      >  1:1014293/1‑151 (MQ=255)
                                                                                                                          tgcgttacatgagagATAAGGATTGAGGAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGt                               <  1:986514/140‑1 (MQ=255)
                                                                                                                                aaCTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTc                          >  1:1421319/1‑150 (MQ=255)
                                                                                                                                 aCTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                         >  1:310448/1‑150 (MQ=255)
                                                                                                                                 aCTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                         >  2:1608469/1‑150 (MQ=255)
                                                                                                                                  cTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                       >  1:710683/1‑151 (MQ=255)
                                                                                                                                  cTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                       >  1:1034741/1‑151 (MQ=255)
                                                                                                                                       agaTTAGGATGGGGGAAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATgg                    <  1:494616/149‑1 (MQ=255)
                                                                                                                                       aaattgggattgggAAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                   <  1:1305457/148‑1 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag                 >  2:735555/1‑151 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag                 >  2:1262765/1‑151 (MQ=255)
                                                                                                                                           tgggttgggggAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac              <  1:597380/142‑1 (MQ=255)
                                                                                                                                           tAGGATGGCGAAAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac              <  1:415356/151‑1 (MQ=255)
                                                                                                                                                       gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  2:1379227/1‑151 (MQ=255)
                                                                                                                                                       gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  2:1408647/1‑151 (MQ=255)
                                                                                                                                                       gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  1:436156/1‑151 (MQ=255)
                                                                                                                                                       gaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaaa  >  1:355099/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAA‑TTATGGAGAAAAATG‑ACAGGGAAAAA‑GG‑AGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAA  >  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: ‑