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

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATT  >  NC_000913/5012‑5311
                                                                                                                                                          |                                                                                                                                                      
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)
                                                                                                                                                         aTAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAAtt  >  1:650934/1‑151 (MQ=255)
                                                                                                                                                          tAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAAtt  >  1:91228/1‑150 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATT  >  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: