Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA42.2% 81.6 / 45.2 84intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (38/10);  new base A (0/35);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 1.20e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/5012‑5306
                                                                                                                                                       |                                                                                                                                                    
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                   <  2:897134/151‑1 (MQ=255)
  tcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                                  >  1:380777/1‑150 (MQ=255)
   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                                 <  1:886584/150‑1 (MQ=255)
    aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTGTTATGAAGAAATTATGGAGAAAAATGACAGGGGAAAAGGAGAAATTCTCAATCAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaaca                                                                                                                                               <  2:538206/151‑1 (MQ=255)
      tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                            >  1:811392/1‑151 (MQ=255)
        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                          >  2:342332/1‑151 (MQ=255)
        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                          >  1:726837/1‑151 (MQ=255)
              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtc                                                                                                                                    >  1:877517/1‑151 (MQ=255)
                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                  >  2:362232/1‑150 (MQ=255)
                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtct                                                                                                                                   >  2:595724/1‑149 (MQ=255)
                  aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCAt                                                                                                                                >  1:815097/1‑151 (MQ=255)
                         cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTa                                                                                                                          >  1:1009624/1‑150 (MQ=255)
                              aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                    >  1:53060/1‑151 (MQ=255)
                                 ccGGGTTTCCTTTTAGGCACCCCGTCTTTTTTAAGAAAAAATAAGGAAGAAAATGAAAAGGGAAAAGAAGGAAATCCAAAAAAAGGTGGGTACTGTAGAATTTAGGGTGGCGAAGAAAACCATCCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                  <  2:709815/150‑1 (MQ=255)
                                      tttGTTTTTATGCAGCCCGGTTTTTTTATAAAAAAATTAGGAAAAAAAATAACAGGAAAAAAGAAAAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATGGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                             <  1:396936/150‑1 (MQ=255)
                                       ttGTTTTTAGGAAGCCCGGTTTTTTTAGGAAAAAATTATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAACTTAGAGATTAGAATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCg                                                                                                            <  1:958953/150‑1 (MQ=255)
                                              taggCACCCCGGTTTTTTTATGAAAAAATTATGAAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGATAATTTAGAAATAGGGATTGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                     <  1:641088/147‑1 (MQ=255)
                                               aggCACCCCGGTTTTTTAAGGAAAAAATTATGAAAAAAAATGACAGGAAAAAAGGAAAAATCCTCAAAAAATGCGGTAACTTAGAGATAAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                     <  1:542236/147‑1 (MQ=255)
                                                tGCACCCCGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGAAAAAATTCTCAAAAAAGGCGGTAACTAAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                   <  1:1250132/150‑1 (MQ=255)
                                                 gCAGCCGGGCTTTTTAAAGAAGAAATAAAGGAAAAAAGAGAGAGGAAAAAGGAAGAACTTATCAATAAATGCGTAATCGAAAAAATAGGTAGTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAc                                                                                                   <  2:352707/149‑1 (MQ=255)
                                                  cAGCCCGGTTTTTTTATGAAGAAATTAGGGAGAAAAAGGACAGGAAAAAAGGAGAAATTCTAAAAAAAGGCGTTAATTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACg                                                                                                  <  1:868225/149‑1 (MQ=255)
                                                    gccggttttttttAAAAAAAATTTATGGAAAAAAAAGAGAGGAAAAAAGGAAAATTTCTAAATAAAGGCGGTATTTTAGAGATTGGGGTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                <  1:976360/142‑1 (MQ=255)
                                                     cccGGTTTTTTTATGAAAAAATTAGGAAAAAAAGTGAGAGGAAAAAAGGAAAAAGCCGAAAAAAATGCGATAACTTAGGGATGAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                             <  1:541807/151‑1 (MQ=255)
                                                        gttttttttAGAAAAAAATTAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAGCTGAGAAATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                           <  2:557063/147‑1 (MQ=255)
                                                         gCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         >  1:414011/1‑151 (MQ=255)
                                                         gCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         >  1:337079/1‑151 (MQ=255)
                                                          cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGCTTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAAAGGGCAATGAt                                                                                         >  2:675600/1‑150 (MQ=255)
                                                          cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         >  1:400415/1‑150 (MQ=255)
                                                           ttttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                        >  1:655941/1‑150 (MQ=255)
                                                           ttttttATGAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAATTTCTAAATAAAGGCGTTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         <  1:462160/149‑1 (MQ=255)
                                                           ttttttATAAAGAAATTAAGGAAAAAAATGACAGGAAAAAAGGAAAAATTCAAAAAAAGTGGTGAAACTAAAATATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                         <  1:723111/149‑1 (MQ=255)
                                                            tttttAAAAAAAAATTAGGGAAAAAAATAACAGGAAAAAAGAAAAAATCCAAAAAAAAGGCGTTAATTTAGAAATTAGGATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                       <  2:631938/150‑1 (MQ=255)
                                                             tttaaaaaaaaaaTTGTGAAAAAAAAGAAAAGGAAAAAAGAAAAATTTCACATATAATGCGGTAATTAAAAAATGAGGATTACGAAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                       <  2:1054447/139‑1 (MQ=255)
                                                                taaaaaaaaaTTATGAAAAAAAGAAAAGAGAAAAAGAGAAAAATTCCAAAAAAATGCTGTAAATAAGAAATGGGAATGGGGATAATAACAAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                   <  2:92167/143‑1 (MQ=255)
                                                                    aaaaaaTTATGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAAAAAATGGGGTAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                  <  1:230731/144‑1 (MQ=255)
                                                                       acaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                           >  2:695818/3‑151 (MQ=255)
                                                                       aaaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                           >  2:42762/1‑151 (MQ=255)
                                                                        aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                          >  1:826765/1‑151 (MQ=255)
                                                                        aaTTAGGGAAAAAAAGGAAAGGAAAAAAGAAAAAATCCAAAAAAAATGCGGAAACTTAGAGATAAGGATGGGGGAGAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                          <  1:621690/151‑1 (MQ=255)
                                                                         aTTATGGACACAATTGACAGGTAAAAAGGAGAAATTCTCGATAGCTGAGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACAGACGCCCTCACCGTGTAATCTCACGCTACCGACCCCTAACGGGCAATGACAAAAagagacac                                                                             >  2:858393/1‑143 (MQ=255)
                                                                               gaaaaaaaTAAGAGGAAAAAAGAAAAAATGCAAAATAAGGGGTGTAATTTGGAGAATAGGAGTGCGGAGAAAACAACCCGGCGTTCTCATCGAGTAATCGCCGGATATCGACCCATAACGGGCAATGATAGAAGGAGTAACCTGTGaa                                                                      <  2:148638/145‑1 (MQ=255)
                                                                                aaaaaaaTGACAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                       <  1:518729/144‑1 (MQ=255)
                                                                                aaaaaaaGGACGGGAAAAAAGAAAAATTTCAAAAAAAATGCGTTAATTTAAAGTTAAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                    <  1:690581/147‑1 (MQ=255)
                                                                                 aaaaaaTAAGAGGAAAAAAGAAAAAATCCTAAAAAATTGGGATACATAAAGTTTGAGATTTGGGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                 <  2:1010439/150‑1 (MQ=255)
                                                                                  aaaacaaacagaaaaaaaGAAAAAATTCTAAATAAAGGCGGTAACTTAGAGATAAAGAGGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                <  2:94090/138‑1 (MQ=255)
                                                                                  aaaaaTGACAGAAAAAAAGAAAATATCAAAAATAAGGGGGATAACTAAGAAATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                  <  1:260649/149‑1 (MQ=255)
                                                                                   aaatagagaGGAAAAGAGAAGAAATTTACAAAAACAGCGGTATTATAGTGATTGGGCGGGCGAAGAAAAACTACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  2:1114136/143‑1 (MQ=255)
                                                                                   aaaaTGACAGAAAAAAAGGACAAATTCTAAAGAAATGCTATAATTTAGATATTAGGATGGCGAAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                               <  2:1168185/151‑1 (MQ=255)
                                                                                     aaggACAGGAAAAAAGAAAAAATTCAAAAAAAAGGGGGTAACTAAGAGTTTAGGATTGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                             <  1:822924/148‑1 (MQ=255)
                                                                                     aaagaagggaaaaaaGGAAAAATTCTAAAAAAAGGGGTTATCTAAGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                             <  1:1212791/141‑1 (MQ=255)
                                                                                        gACAGAAAAAAAGGAAAAATCCTCAAAAAATGCGGTAATTTAAAAATTAGGATTGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                           <  1:1008343/150‑1 (MQ=255)
                                                                                          cAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTCCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                         >  1:684761/1‑150 (MQ=255)
                                                                                           aGGGAAAAAGGAGAAATTCTCAATAAATGCGGCAACTTAGAGATTAGGATTGCGGAGAATAACAATCGACGTTCTCATAGAGTAATCTCCGGATATCGACCCATAAAGGGCAATGATAAAAGGAGTACCCTGTGAAAAAGATGCAATCTa                                                        >  2:1031834/1‑150 (MQ=255)
                                                                                            ggaaaaaaGGAAAAATTCTAAATAAATGCGGTATTTTAGAGTTTAGATTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                      <  1:326578/148‑1 (MQ=255)
                                                                                              gAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    >  1:67190/1‑151 (MQ=255)
                                                                                               aaaaaGGAGAAATTCTAAATAAATGGGGTAACTTAGAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTa                                                   <  1:27418/151‑1 (MQ=255)
                                                                                               aaaaaGGAGAAATTCTAAATAAAGGTGTTAACTAAAAGATTGGGATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                          <  2:108896/144‑1 (MQ=255)
                                                                                               aaaaaGGAAAAATTCTCAATAAATGCGTTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    <  1:797141/150‑1 (MQ=255)
                                                                                               aaaaaGAAAAAATTCTCAATAAGGGCGGTAATTTAGAAATTGGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                    <  2:380592/150‑1 (MQ=255)
                                                                                                   aGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCg                                               >  1:1141880/1‑151 (MQ=255)
                                                                                                      gagaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                >  2:384767/1‑148 (MQ=255)
                                                                                                       aaaaaTTCTAAAAAAATGCGGTAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                            <  1:13761/149‑1 (MQ=255)
                                                                                                        gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                           >  1:1093323/1‑151 (MQ=255)
                                                                                                        gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                           >  1:753311/1‑151 (MQ=255)
                                                                                                         aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                          >  2:619181/1‑151 (MQ=255)
                                                                                                          aaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACttt                                         >  2:441976/1‑151 (MQ=255)
                                                                                                            ttCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                       >  1:1228039/1‑151 (MQ=255)
                                                                                                             tctcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                      >  1:20832/1‑151 (MQ=255)
                                                                                                              ctcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCACt                                     >  2:55031/1‑151 (MQ=255)
                                                                                                               tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                    >  1:961236/1‑151 (MQ=255)
                                                                                                                aaaTAAATGATTTGACATAAAGTTTAGGATTGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                    <  2:725484/149‑1 (MQ=255)
                                                                                                                 aaaaaaaGCGGAAATTAAGAGTTAGGGATGCGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                    <  2:1217761/146‑1 (MQ=255)
                                                                                                                  aaaaaTGCGGTAATTAAGAGATTAGGTTGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 <  1:1210975/149‑1 (MQ=255)
                                                                                                                  aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                 >  2:1120604/1‑151 (MQ=255)
                                                                                                                    aattGCGTTAATTTAGAGATTAGGTTTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                    <  1:762644/143‑1 (MQ=255)
                                                                                                                    aaaTGCGGTAATTTAGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                <  1:411138/150‑1 (MQ=255)
                                                                                                                        gCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACAGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTAGCACTTTACCTGGTTCTGGTc                           >  1:486406/1‑151 (MQ=255)
                                                                                                                                ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGATCTGGTCGCTCCCAt                   >  1:147046/1‑151 (MQ=255)
                                                                                                                                       aTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc             >  1:964482/1‑151 (MQ=255)
                                                                                                                                        ttAGGATGGGGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                  <  1:312499/145‑1 (MQ=255)
                                                                                                                                         tgggATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATTATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACGCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca            <  2:53047/148‑1 (MQ=255)
                                                                                                                                            gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg        >  2:261459/1‑151 (MQ=255)
                                                                                                                                             ggttGCGAAGAAAAAATACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCGCTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg         <  2:1222273/148‑1 (MQ=255)
                                                                                                                                                  cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg   >  1:920573/1‑151 (MQ=255)
                                                                                                                                                    gagaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:1058571/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/5012‑5306

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: