Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1620TA38.1% 124.7 / 40.9 85intergenic (+142/‑72)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (40/12);  new base A (0/32);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 2.13e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATT  >  NC_000913/5012‑5311
                                                                                                                                                          |                                                                                                                                                      
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAAt                                                                                                                                                        <  2:897134/151‑1 (MQ=255)
  tcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATa                                                                                                                                                       >  1:380777/1‑150 (MQ=255)
   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATaa                                                                                                                                                      <  1:886584/150‑1 (MQ=255)
    aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTGTTATGAAGAAATTATGGAGAAAAATGACAGGGGAAAAG‑GAGAAATTCTCAATCAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATaaca                                                                                                                                                    <  2:538206/151‑1 (MQ=255)
      tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATaacaa‑c                                                                                                                                                 >  1:811392/1‑151 (MQ=255)
        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑ccg                                                                                                                                               >  2:342332/1‑151 (MQ=255)
        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑ccg                                                                                                                                               >  1:726837/1‑151 (MQ=255)
              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTtc                                                                                                                                         >  1:877517/1‑151 (MQ=255)
                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTtctc                                                                                                                                       >  2:362232/1‑150 (MQ=255)
                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTtct                                                                                                                                        >  2:595724/1‑149 (MQ=255)
                  aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCAt                                                                                                                                     >  1:815097/1‑151 (MQ=255)
                         cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTa                                                                                                                               >  1:1009624/1‑150 (MQ=255)
                              aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTcc                                                                                                                         >  1:53060/1‑151 (MQ=255)
                                 ccGGGTTTCCTTTTAGGCACCCCGTCTTTTTTAAGAA‑AAAATAAGGAAGAAAATGAAAAGGGAAAAG‑AAGGAAATCCAAAAAAAGGTGGGTACTGTAGAATTTAGG‑GTGGCGA‑AGAAAACCAT‑CCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                       <  2:709815/150‑1 (MQ=255)
                                      tttGTTTTTATGCAGCCCGGTTTTTTTATAAAAAAATTAGGAAAAAAAATAACAGGAAAAAAG‑AAAAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATGGGGG‑AGAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                                  <  1:396936/150‑1 (MQ=255)
                                       ttGTTTTTAGGAAGCCCGGTTTTTTTAGGAAAAAATTATGGAAAAAAATGAAAGGAAAAAAG‑GAAAAATTCTAAATAAATGCGGTAACT‑TAGAGATTAGA‑ATGGCGG‑AAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                                 <  1:958953/150‑1 (MQ=255)
                                              taggCACCCCGGTTTTTTTATGAAAAAATTATGAAAAAAAAGGACAGGAAAAAAG‑GAAAAATTCTAAATAAATGCGATAATT‑TAGAAATAGGG‑ATTGGGG‑AGAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                          <  1:641088/147‑1 (MQ=255)
                                               aggCACCCCGGTTTTTTAAGGAAAAAATTATGAAAAAAAATGACAGGAAAAAAG‑GAAAAATCCTCAAAAAATGCGGTAACT‑TAGAGATAAGG‑TTTGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                          <  1:542236/147‑1 (MQ=255)
                                                tGCACCCCGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAG‑AAAAAATTCTCAAAAAAGGCGGTAACT‑AAGAGATTAGG‑ATTGGGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                        <  1:1250132/150‑1 (MQ=255)
                                                 gCAGCCGGGCTTTTTAAAGAAGAAATAAAGGAAAAAAGAGAGAGGAAAAAGG‑AAGAACTTATCAATAAATGC‑GTAATC‑GAAAAAATAGGTAGTGGGG‑AAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                        <  2:352707/149‑1 (MQ=255)
                                                  cAGCCCGGTTTTTTTATGAAGAAATTAGGGAGAAAAAGGACAGGAAAAAAG‑GAGAAATTCTAAAAAAAGGCGTTAATT‑TAGAGATTAGG‑ATTGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                       <  1:868225/149‑1 (MQ=255)
                                                    gccggttttttttAAAAAAAATTTATGGAAAAAAAAGAGAGGAAAAAAG‑GAAAATTTCTAAATAAAGGCGGTATTT‑TAGAGATTGGG‑GTGGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                     <  1:976360/142‑1 (MQ=255)
                                                     cccGGTTTTTTTATGAAAAAATTAGGAAAAAAAGTGAGAGGAAAAAAG‑GAAAAAGCCGAAAAAAATGCGATAACT‑TAGGGATGAGG‑ATGGGGG‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                                  <  1:541807/151‑1 (MQ=255)
                                                        gttttttttAGAAAAAAATTAGGAAAAAAAATGACAGGAAAAAAG‑GAAAAATTCTAAATAAAGGCGGTAGCT‑GAGAAATTAGG‑ATTGGGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                                <  2:557063/147‑1 (MQ=255)
                                                         gCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                              >  1:414011/1‑151 (MQ=255)
                                                         gCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                              >  1:337079/1‑151 (MQ=255)
                                                          cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑CTTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAAGGGCAATGAt                                                                                              >  2:675600/1‑150 (MQ=255)
                                                          cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                              >  1:400415/1‑150 (MQ=255)
                                                           ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                             >  1:655941/1‑150 (MQ=255)
                                                           ttttttATGAAAAAATTAGGGAAAAAAATGACAGGAAAAAAG‑GAAAATTTCTAAATAAAGGCGTTAACT‑TAGAGATTAGG‑ATGGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                              <  1:462160/149‑1 (MQ=255)
                                                           ttttttATAAAGAAATTAAGGAAAAAAATGACAGGAAAAAAG‑GAAAAATTCAAAAAAAGTGGTGAAACT‑AAAATATTAGG‑ATTGGGG‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                              <  1:723111/149‑1 (MQ=255)
                                                            tttttAAAAAAAAATTAGGGAAAAAAATAACAGGAAAAAAG‑AAAAAATCCAAAAAAAAGGCGTTAATT‑TAGAAATTAGG‑ATTGCGA‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                            <  2:631938/150‑1 (MQ=255)
                                                             tttaaaaaaaaaaTTGTGAAAAAAAAGAAAAGGAAAAAAG‑AAAAATTTCACATATAATGCGGTAATT‑AAAAAATGAGG‑ATTACGA‑AAAAAACCAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                            <  2:1054447/139‑1 (MQ=255)
                                                                taaaaaaaaaTTAT‑GAAAAAAAGAAAAGAGAAAAAGAGAAAAATTCCAAAAAAATGCTGTAAAT‑AAGA‑AATGGG‑AATGGGG‑ATAATAACAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                        <  2:92167/143‑1 (MQ=255)
                                                                    aaaaaaTTATGGAAAAAAATGAAAGGAAAAAAG‑GAAAAATTCTCAAAAAATGGGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                       <  1:230731/144‑1 (MQ=255)
                                                                       acaTTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                                >  2:695818/3‑151 (MQ=255)
                                                                       aaaTTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                                >  2:42762/1‑151 (MQ=255)
                                                                        aaTTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                               >  1:826765/1‑151 (MQ=255)
                                                                        aaTTAGGGAAAAAAAGGAAAGGAAAAAAG‑AAAAAATCCAAAAAAAATGCGGAAACT‑TAGAGATAAGG‑ATGGGGG‑AGAAAAACAG‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                               <  1:621690/151‑1 (MQ=255)
                                                                         aTTATGGACACAATTGACAGGTAAAAAG‑GAGAAATTCTCGATAGCTGAGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CAGACGCCCTCACCGTGTAATCTCACGCTACCGACCCCTAACGGGCAATGACAAAAagagacac                                                                                  >  2:858393/1‑143 (MQ=255)
                                                                               gaaaaaaaTAAGAGGAAAAAAG‑AAAAAATGCAAAATAAGGGGTGTAATT‑TGGAGAATAGG‑AGTGCGG‑AGAA‑AACAACCCGGCGTTCTCATCGAGTAATCGCCGGATATCGACCCATAACGGGCAATGATAGAAGGAGTAACCTGTGaa                                                                           <  2:148638/145‑1 (MQ=255)
                                                                                aaaaaaaTGACAGGAAAAAAG‑GAAAAATTCTAAAAAAATGCGTTAATT‑TAGAGATTAGG‑ATGGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                            <  1:518729/144‑1 (MQ=255)
                                                                                aaaaaaaGGACGGGAAAAAAG‑AAAAATTTCAAAAAAAATGCGTTAATT‑TAAAGTTAAGG‑ATGGCGG‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                         <  1:690581/147‑1 (MQ=255)
                                                                                 aaaaaaTAAGAGGAAAAAAG‑AAAAAATCCT‑AAAAAATTGGGATACA‑TAAAGTTTGAG‑ATTTGGG‑AAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                      <  2:1010439/150‑1 (MQ=255)
                                                                                  aaaacaaacagaaaaaaaG‑AAAAAATTCTAAATAAAGGCGGTAACT‑TAGAGATAAAG‑AGGGCGA‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                     <  2:94090/138‑1 (MQ=255)
                                                                                  aaaaaTGACAGAAAAAAAG‑AAAATATCAAAAATAAGGGGGATAACT‑AAGAAATTAGG‑ATGGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                       <  1:260649/149‑1 (MQ=255)
                                                                                   aaatagagaGGAAAAGAG‑AAGAAATTTACAAAAACAGCGGTATTA‑TAGTGATTGGG‑CGGGCGA‑AGAAAAACTA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                    <  2:1114136/143‑1 (MQ=255)
                                                                                   aaaaTGACAGAAAAAAAG‑GACAAATTCTAAAGAAATGCTATAATT‑TAGATATTAGG‑ATGGCGA‑AAAATAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                    <  2:1168185/151‑1 (MQ=255)
                                                                                     aaggACAGGAAAAAAG‑AAAAAATTCAAAAAAAAGGGGGTAACT‑AAGAGTTTAGG‑ATTGGGG‑AAAATAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                                  <  1:822924/148‑1 (MQ=255)
                                                                                     aaagaagggaaaaaaG‑GAAAAATTCTAAAAAAAGGGGTTATCT‑AAGAGATTAGG‑ATTGGGG‑AAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                                  <  1:1212791/141‑1 (MQ=255)
                                                                                        gACAGAAAAAAAG‑GAAAAATCCTCAAAAAATGCGGTAATT‑TAAAAATTAGG‑ATTGGGG‑AAAATAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                                <  1:1008343/150‑1 (MQ=255)
                                                                                          cAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTCCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                              >  1:684761/1‑150 (MQ=255)
                                                                                           aGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGCAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑TCGACGTTCTCATAGAGTAATCTCCGGATATCGACCCATAAAGGGCAATGATAAAAGGAGTACCCTGTGAAAAAGATGCAATCTa                                                             >  2:1031834/1‑150 (MQ=255)
                                                                                            ggaaaaaaG‑GAAAAATTCTAAATAAATGCGGTATTT‑TAGAGTTTAGA‑TTTGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                           <  1:326578/148‑1 (MQ=255)
                                                                                              gAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                         >  1:67190/1‑151 (MQ=255)
                                                                                               aaaaaG‑GAGAAATTCTAAATAAATGGGGTAACT‑TAGAGATTAGG‑ATGGGGA‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTa                                                        <  1:27418/151‑1 (MQ=255)
                                                                                               aaaaaG‑GAGAAATTCTAAATAAAGGTGTTAACT‑AAAAGATTGGG‑ATTGCGA‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                               <  2:108896/144‑1 (MQ=255)
                                                                                               aaaaaG‑GAAAAATTCTCAATAAATGCGTTAACT‑TAGAGATTAGG‑ATTGGGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                         <  1:797141/150‑1 (MQ=255)
                                                                                               aaaaaG‑AAAAAATTCTCAATAAGGGCGGTAATT‑TAGAAATTGGG‑ATGGGGG‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                         <  2:380592/150‑1 (MQ=255)
                                                                                                   aG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                                    >  1:1141880/1‑151 (MQ=255)
                                                                                                      gagaAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                     >  2:384767/1‑148 (MQ=255)
                                                                                                       aaaaaTTCTAAAAAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                                 <  1:13761/149‑1 (MQ=255)
                                                                                                        gaAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                                >  1:1093323/1‑151 (MQ=255)
                                                                                                        gaAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                                >  1:753311/1‑151 (MQ=255)
                                                                                                         aaaTTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                               >  2:619181/1‑151 (MQ=255)
                                                                                                          aaTTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                              >  2:441976/1‑151 (MQ=255)
                                                                                                            ttCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                            >  1:1228039/1‑151 (MQ=255)
                                                                                                             tctcAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                           >  1:20832/1‑151 (MQ=255)
                                                                                                              ctcAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCACt                                          >  2:55031/1‑151 (MQ=255)
                                                                                                               tcAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                         >  1:961236/1‑151 (MQ=255)
                                                                                                                aaaTAAATGATTTGACA‑TAAAGTTTAGG‑ATTGGGG‑AAAAAAACAC‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                         <  2:725484/149‑1 (MQ=255)
                                                                                                                 aaaaaaaGCGGAAATT‑AAGAG‑TTAGG‑GATGCGGAAAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                         <  2:1217761/146‑1 (MQ=255)
                                                                                                                  aaaaaTGCGGTAATT‑AAGAGATTAGG‑TTGGGGG‑AAAAAAAAAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                      <  1:1210975/149‑1 (MQ=255)
                                                                                                                  aTAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                      >  2:1120604/1‑151 (MQ=255)
                                                                                                                    aattGCGTTAATT‑TAGAGATTAGG‑TTTGCGG‑AAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                         <  1:762644/143‑1 (MQ=255)
                                                                                                                    aaaTGCGGTAATT‑TAGAGATTAGG‑ATGGCGG‑AAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                     <  1:411138/150‑1 (MQ=255)
                                                                                                                        gCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CAGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTAGCACTTTACCTGGTTCTGGTc                                >  1:486406/1‑151 (MQ=255)
                                                                                                                                t‑tAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGATCTGGTCGCTCCCAt                        >  1:147046/1‑151 (MQ=255)
                                                                                                                                       aTTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc                  >  1:964482/1‑151 (MQ=255)
                                                                                                                                        ttAGG‑ATGGGGA‑AAAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                       <  1:312499/145‑1 (MQ=255)
                                                                                                                                         tggg‑ATGGCGG‑AAAAAAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATTATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACGCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca                 <  2:53047/148‑1 (MQ=255)
                                                                                                                                            g‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg             >  2:261459/1‑151 (MQ=255)
                                                                                                                                             ggttGCGA‑AGAAAAAATA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCGCTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg              <  2:1222273/148‑1 (MQ=255)
                                                                                                                                                  cGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg        >  1:920573/1‑151 (MQ=255)
                                                                                                                                                    g‑agaATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg       >  1:1058571/1‑150 (MQ=255)
                                                                                                                                                         aTAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAAtt  >  1:650934/1‑151 (MQ=255)
                                                                                                                                                          tAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAAtt  >  1:91228/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAG‑GAGAAATTCTCAATAAATGCGGTAACT‑TAGAGATTAGG‑ATTGCGG‑AGAATAACAA‑CCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATT  >  NC_000913/5012‑5311

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