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

ATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5010‑5303
                                                                                                                                                       |                                                                                                                                                 
atcatcAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                 <  1:945554/151‑1 (MQ=255)
  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGaa                                                                                                                                                <  2:488408/150‑1 (MQ=255)
  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                               >  2:582413/1‑151 (MQ=255)
  catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                               >  2:91786/1‑151 (MQ=255)
     caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                            >  2:665383/1‑151 (MQ=255)
     caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                            >  1:381416/1‑151 (MQ=255)
                aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAATTTCTCAATAAATGCGGTCACTTAGAGATTAGGATTGCGGAGAATAACAAacgccg                                                                                                                                    >  2:284228/1‑148 (MQ=255)
                 ttCATTATCTAAATCAGGCGGGGTTGGGTTTTAGGCAGCCCGGTTTTTTTATGAAGAAATGATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTAAATAAATGCGGTAATTTAGAGATGAGGATTGCGGAGAATAACAAccgccg                                                                                                                                    <  2:100570/147‑1 (MQ=255)
                     ttATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc                                                                                                                            >  2:898890/1‑151 (MQ=255)
                              tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                   >  2:351667/1‑151 (MQ=255)
                              tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                   >  1:656529/1‑151 (MQ=255)
                               cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                  >  2:929176/1‑151 (MQ=255)
                                 ggCGGGTTTTGTTTTTATGCGCCCCGGTTTTTTTATGAAAAAATTAGGGAAAAAAATGACGGGAAAAAAGGAAAATTTCCAAATAAATGCGGAAATTAAGAGATTAGGATTGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                  <  1:619872/149‑1 (MQ=255)
                                 ggCGGGGTTTGTTTTTATGCACCCCGGTTTTTTTATGAAAAAATTAGGAAAAAAAATGACAGGAAAAAAGGAAAATTTCCAAAAAAATGGGGTAACTTAGAGATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAAtct                                                                                                                   <  1:224810/148‑1 (MQ=255)
                                  gCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                               >  1:228283/1‑151 (MQ=255)
                                        tttttttttACGCACCCCGGTTTTTTTATAAAAAAATTATGAAAAAAAATAAAAGGAAAAAAGGAAAAATTCTCAATAAATGCGTTACCTAAGAGATAAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                          <  2:904873/145‑1 (MQ=255)
                                        ttgttttttCTGCAGCCCGGTTTTTTAAAAAAAAATTTAGGAAAAAAAATAACAGGAAAAAAGGAGAAATTCTAAAAAAAGGCGGTAACTGAGAGATTAGTATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                          <  2:437378/145‑1 (MQ=255)
                                            ttttaaagCAACCCGGTTTTTTAATGAAGAAATTAAGAAGAAAAAGGGAGGGAAAAAAGGAAAAATTTCAAAAAAGGTTGGTATCATAGAGTTTGGGGTTGGGAGAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGa                                                                                                        <  2:687465/141‑1 (MQ=255)
                                              tttATGCACCCCGGTTTTTTTATAAAAAAATAAGGAAAAAAAAGGACGGGAAAAAAGAAAAAATTCCCAAAAAATGGGGAAATTGAGAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                    <  2:753980/150‑1 (MQ=255)
                                               ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCa                                                                                                    >  2:677944/1‑149 (MQ=255)
                                                       ccgggtttttttATGAAAAAATAAGGAAAAAAAAGGACGGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATTGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                           <  2:208873/144‑1 (MQ=255)
                                                        ccGGTTTTTTAAAAAAAAAATTAGGGAAAAAAATGAAGGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAATTAAGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCa                                                                                           <  1:1012761/149‑1 (MQ=255)
                                                        ccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                         >  2:783360/1‑151 (MQ=255)
                                                         cGGCTTTTTAATAAAAAAATTATGGAAAAAAATGACAGGGAAAAAGGAAAAATTCAAAATAACCGCGATAAATTAAATCTTAGGATGGGGCAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                        <  1:438342/151‑1 (MQ=255)
                                                          ggtttttttAGAAAAAAATAAGGGAAAAAAATGACGGGAAAAAGGGAAAAATTCTAAATAAATGGGGTAAATTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                        <  2:461037/147‑1 (MQ=255)
                                                             ttttttATGAAAAAATTATGAAAAAAAAGGACAGGAAAAAAGAAAAAATCCTAAAAAAAGGCGTTAATTTAGAAATTAGAATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                      <  2:868436/149‑1 (MQ=255)
                                                                tttATGCAGAAATTTTGGAGAAAAATGAAAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCACCGATCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                        >  2:540103/1‑144 (MQ=255)
                                                                  tATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                >  2:974806/1‑150 (MQ=255)
                                                                   agaaaaaaaTTAGGGAGAAAAATGAAAGGAAAAAAGAAAAAATTACAAAAAAAGGCGGTAATTAAGAGATTAGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                <  2:303398/143‑1 (MQ=255)
                                                                     gaagaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                             >  1:137985/1‑150 (MQ=255)
                                                                      aagaaATTATGAAAAAAAAGGACAGGAAAAAAGAAGAATTTCAAAATAAATGCGTTACCTTAGAGATTAGGATTGCGGAAAATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                           <  1:1088627/151‑1 (MQ=255)
                                                                       aaaaaTTAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAACTAAGAGATTAGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                          <  1:506640/149‑1 (MQ=255)
                                                                       aaaaaTTAAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCAAAATAAAGGGGATAAATAAAAGATAAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                             <  1:521327/147‑1 (MQ=255)
                                                                        aaaaTAATGAAAAAAAATGACAGGAAAAAGGGAAAAATTCTAAATAAAGGCTGTAATTTAGAGATTAGGATTGAGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                         <  2:350181/150‑1 (MQ=255)
                                                                         aaaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                        >  2:238847/1‑151 (MQ=255)
                                                                          aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                       >  1:728806/1‑151 (MQ=255)
                                                                            ttATGAAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                      <  1:91806/150‑1 (MQ=255)
                                                                                 gagaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                 >  1:1081365/1‑150 (MQ=255)
                                                                                     aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                 >  2:306989/1‑147 (MQ=255)
                                                                                     aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                 >  1:717520/1‑147 (MQ=255)
                                                                                      aaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                            >  2:987584/1‑151 (MQ=255)
                                                                                       aaaTGAAAGGAAAAAAGGAGAATTTCTAAATAAATGCGTTACCTAAGAGATTAGGTTGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                           <  1:314271/151‑1 (MQ=255)
                                                                                         aTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                          >  2:330859/1‑150 (MQ=255)
                                                                                            aCAGGAAAAAAGGAGAAATTCCAAAAAAAGGGGGAAATTTAGAGATTGGGATTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                        <  1:255562/149‑1 (MQ=255)
                                                                                            aCAGGAAAAAAGGAAAATTTCTAAAAAAATGCGTTAACTAAGAGTTTGGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCGGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                       <  1:438951/150‑1 (MQ=255)
                                                                                            aCAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                        <  1:869757/149‑1 (MQ=255)
                                                                                               gggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                   >  1:686955/1‑151 (MQ=255)
                                                                                               gggAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGCGAAAAAGATGCAATCTATc                                                   >  1:392231/1‑151 (MQ=255)
                                                                                               ggaaaaaaGAAAAATTTCTAAATAAATGGGTTATCTTAGAGGTAAGGTTGGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                       <  1:629423/144‑1 (MQ=255)
                                                                                                 aaaaaaGGAGAAATTCACAAAAAAGGCGTTAACTTAGAGATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGt                                                 <  2:409977/150‑1 (MQ=255)
                                                                                                  aaaaaGGAGAATTTTAAAAAAAGTGCGGTAATTAAGAGATTAGGATGGGGGAAAAAAACAGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                  <  2:895575/149‑1 (MQ=255)
                                                                                                  aaaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACTTGTGAAAAAGATGCAATCTATCGt                                                 >  2:738693/1‑150 (MQ=255)
                                                                                                   aaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                               >  1:797598/1‑151 (MQ=255)
                                                                                                      agaaaaaaTTCTCAAAAAAGGGGGTAATTAAGAGATTAGGATTTGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                             <  2:228220/145‑1 (MQ=255)
                                                                                                         agaAATTCTAAATAAAAGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACAGCAGTTCTCATCGAGTAATCTCCAGATATCGACCCATAACGAGCAGTGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                         >  2:158629/1‑151 (MQ=255)
                                                                                                          gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                        >  2:1065446/1‑151 (MQ=255)
                                                                                                           aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                       >  1:525943/1‑151 (MQ=255)
                                                                                                                 tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                 >  1:824748/1‑151 (MQ=255)
                                                                                                                  cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                >  1:484902/1‑151 (MQ=255)
                                                                                                                       aaggCGGTAACTTAGAAATTGGGATTGCGAAAAATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                           <  2:121400/148‑1 (MQ=255)
                                                                                                                            cGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCg                       >  2:110144/1‑151 (MQ=255)
                                                                                                                                aaCTTAGAAATTAGAATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                   <  1:331949/151‑1 (MQ=255)
                                                                                                                                aaCTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCt                     <  1:306965/149‑1 (MQ=255)
                                                                                                                                  cTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                  >  2:421317/1‑150 (MQ=255)
                                                                                                                                    tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg               >  2:200054/1‑151 (MQ=255)
                                                                                                                                    tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATATCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTACCCTGTGAAAAAGATCCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg               >  2:27127/1‑151 (MQ=255)
                                                                                                                                       agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc             >  1:51032/1‑150 (MQ=255)
                                                                                                                                       agaTTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc             <  1:823152/150‑1 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc             >  1:648323/1‑149 (MQ=255)
                                                                                                                                          ttAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc          >  2:286034/1‑150 (MQ=255)
                                                                                                                                            aGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac        >  1:762537/1‑150 (MQ=255)
                                                                                                                                              gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg     >  1:6272/1‑151 (MQ=255)
                                                                                                                                                  ggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGccg  <  1:899182/148‑3 (MQ=255)
                                                                                                                                                  ggggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc    <  1:138085/146‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
ATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTG  >  NC_000913/5010‑5303

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: