Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA47.9% 78.8 / 75.3 96intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (46/4);  new base A (0/46);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 3.93e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5013‑5307
                                                                                                                                                        |                                                                                                                                                    
atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                   >  1:1372267/1‑150 (MQ=255)
atcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATa                                                                                                                                                  <  2:865363/151‑1 (MQ=255)
  caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                                <  1:1359778/151‑1 (MQ=255)
    gTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaaca                                                                                                                                               >  2:97409/1‑150 (MQ=255)
     tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                             >  1:86342/1‑151 (MQ=255)
       aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                           >  2:550388/1‑151 (MQ=255)
          tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgccg                                                                                                                                        >  2:871731/1‑151 (MQ=255)
           cTATTCATTATCTCAATCAGGCAGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGAGGAGAATAACAAccgccg                                                                                                                                        >  1:552401/1‑150 (MQ=255)
              ttCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtc                                                                                                                                     >  2:478106/1‑150 (MQ=255)
                cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                   >  1:1634525/1‑150 (MQ=255)
                  ttATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc                                                                                                                                >  1:690097/1‑151 (MQ=255)
                     tctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATc                                                                                                                                >  1:141793/1‑148 (MQ=255)
                      ctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                            >  2:1597623/1‑151 (MQ=255)
                      ctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                            >  2:425331/1‑151 (MQ=255)
                        cAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTa                                                                                                                           >  1:1606995/1‑150 (MQ=255)
                         aaTCAGGCGGGTTTTTTTTTTAGGCACCCCGTTTTTTTTATAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGAAAAAATTCTAAAAAAATGCGTTAATTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTa                                                                                                                           <  1:1515596/149‑1 (MQ=255)
                           tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                       >  1:210209/1‑151 (MQ=255)
                              ggCGGGTTTTGTTTTTATGCAGCCCGTTTTTTTTAAAAAAAAATTAGGAAAAAAAATGAGAGGAAAAAAGGAAAAATTCTAAAAAAAGGGGTTATTTTAGAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                    <  1:1163378/151‑1 (MQ=255)
                              ggCGGGTTTTGTTTTTATGCAGCCCGTTTTTTTAATAAAAAAATTAGGAAAAAAAAGAAGAGGAAAAAAGGAAAAATTCAAAAAAAAGGGGGTAACTAAGAGATTAGGAGGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                       <  1:642155/148‑1 (MQ=255)
                              ggCGGGTTTTGTTTTTAGGCACCCCGGTTTTTTAATAAAAAAATTAGGAAAAAAAAGTAAGGGAAAAAAGGAAAATTTCTAAAAAAGGGTGTTAATTTAGAAATGGGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                    <  1:177190/151‑1 (MQ=255)
                              ggCGGGGTTTGTTTTTATGCAGCCCGGTTTTTTTATGAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAATTTCTCAAAAAAGGCGTTAATTAAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                       <  2:1245112/148‑1 (MQ=255)
                               gcgggttttGTTTTTATGCAGCCCGTTTTTTTTAGGAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAGTGCGGTAACTTGGAGATTAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCgg                                                                                                                   <  2:1500294/145‑1 (MQ=255)
                                  ggttttGTTTTTATGCACCCCGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAACTTAGAGATTAGGTTGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                                <  2:896820/148‑1 (MQ=255)
                                  gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                 >  2:870296/1‑150 (MQ=255)
                                     ttttCTTTTATGCCGCCGGTTTTTTTTATAAAAAATTTATGAAAAAAAATAACGGGAAAAAAGGAAAAATTCAAAAAAAGTGCGTTACCTTAGAAATTGGTTTTGCGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              <  1:704897/150‑1 (MQ=255)
                                     tttGTTTTTAGGCAGCCGGTTTTTTTTATAAAAAATTTATGAAAAAAAAGGAAAGGAAAAAAGAAAAAATCCTAAAAAAATGGGGTAACTAAGAGATAAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                              <  2:1402660/150‑1 (MQ=255)
                                       tGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                           >  2:1331786/1‑151 (MQ=255)
                                          ttttATGCAGCCCGGTTTTTTTATAAAAAAATTAGGAAAAAAAAGGACAGAAAAAAAGGAAAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                           <  1:637527/148‑1 (MQ=255)
                                              aggCAGCCGGGTTTTTTAATAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGGAAAATTTCTAAAAAAATGCGGTAACTTAGAGATTGGGATTGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      <  2:1411125/147‑1 (MQ=255)
                                                 cAGCCCTGTTTTTTTAGGAAAAATTGTTGGAAAAAAAAGACAGGAAAAAGGGAGAAATCCTAAAAAAGGGGGTTATGTCGGATATTAGTGTGGGGAAAAAAAACATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                      <  2:522845/146‑1 (MQ=255)
                                                    cccGGTTTTTTTAAAAAAAAATTAGGAAAAAAAATGAAAGGAAAAAAGGAAAATTTCTAAAAAAAGGGGGTACTTTAGAGATAAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                              <  2:585568/151‑1 (MQ=255)
                                                     ccGGTTTTTTTTAAAAAAAAATAAGGAAAAAAAAGGAAGGGAAAAAAGAAAAATTTCTAAAAAATTGGGTTAATTAAAAGTTTGGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                              <  1:1633768/151‑1 (MQ=255)
                                                     ccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                            >  1:904713/1‑151 (MQ=255)
                                                      cGGTTTTTTTAAAAAGAAATTAAGAAAAAAAAAAAAAGGAAAAAGGGATAAATTCCGAAAAAGGCGGATAAGTTATGGATTGGGATGGCGAAAAACAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                                 <  2:170581/146‑1 (MQ=255)
                                                       ggCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                            >  2:16794/1‑150 (MQ=255)
                                                       ggCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                            >  1:871000/1‑150 (MQ=255)
                                                         tttttttATGAAGAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAAAGGCGGTAACTGAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                            <  1:1682584/147‑1 (MQ=255)
                                                         tttttttAAAAAAAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAAGGCGGTAACTTAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                            <  2:810274/147‑1 (MQ=255)
                                                          ttttttATAAAAAAATTAGGGAAAAAAATGACGGGAAAAAAGGAAAAATTCTAAATAAAGGCGTTAATTTAGAGATTAGGTTGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                         <  2:241402/150‑1 (MQ=255)
                                                             tttaaaaaaaaaTTAGGAAAAAAAAGGAAGGGAAAAAAGGAAAATTTCAAAAAAAGTGGGTTAATTTAGAGATTAGAATGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                      <  2:627313/141‑1 (MQ=255)
                                                                agaaaaaaaTTAGGAAAAAAAAGGAAGGGAAAAAAGGAAAATTTCTAAAAAAAGGGGTTACTTAAAAAATAAGGTTGGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                   <  2:671349/144‑1 (MQ=255)
                                                                aagaaaaaATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAAAGGCGGTAACTTAGAGATTAGGATGGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                   <  1:1117689/148‑1 (MQ=255)
                                                                aaaaagaaATTAGAGAAAAAAAAAAAAGGAAAAAAGGAAAAATTCAAAAAAATTGCGATACTTTAAAAATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                                  <  1:170815/148‑1 (MQ=255)
                                                                aaaaaaaaaTTAGGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAATGGCGGTAACTTAGAAATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    <  1:1206494/143‑1 (MQ=255)
                                                                aTGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    >  1:1414507/1‑149 (MQ=255)
                                                                 atgaaaaaTTAAGGAAAAAAAACGAGAGGGAAAAAGGGAAATTTCTCAATAAAGGGGGAAAGTGAGATATCGGGTTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAg                                                                                     <  1:442178/143‑1 (MQ=255)
                                                                   aaaaaaTTATGAAAAAAAATGACGGGAAAAAAGGAGAAATTCAAAAAAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAg                                                                                  <  2:788088/145‑1 (MQ=255)
                                                                   aaaaaaTTAGGAAAAAAAAAGAAAGGAAAAAAGGAAAAATTCAAAAAAAATGCGGTAATTTAGAGATTAGGATGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTaa                                                                               <  1:1437497/148‑1 (MQ=255)
                                                                   aaaaaaTTAAGAAAAAAAAGGAAAGAAAAAAAGAAAAAATTCCAAAAAATTGCTATATTATAGAAATTAGATTGGAGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    <  1:1147778/143‑1 (MQ=255)
                                                                          ataaggaagaaaaAGACAAGGGAAAAGAGAAAATTCCTCAATGAACGGGGAAATTTAGAGATAAGACTGGCGAAAAAAACCCGCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                          <  1:216600/142‑1 (MQ=255)
                                                                           aTGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                        >  2:525535/1‑150 (MQ=255)
                                                                                aaaaaaaaGAAAGAAAAAAAGAAAAAATTCAAAAAAAATGCGGTAATTTAGAAATAAGGTTGGAGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAg                                                                   <  1:1506012/149‑1 (MQ=255)
                                                                                  aaaaaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                     >  1:287269/1‑147 (MQ=255)
                                                                                  aaaaaTGACAGAAAAAAAGGAGAATTTCTAAATAATTGGGGTAACTTAGAGATGAGGATGGCGAAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                 <  2:1313849/151‑1 (MQ=255)
                                                                                   aaaaTGAAGGGAAAAAAGGAAAAATTCAAAAAAAAGGCGGAAATTAAGAGATTGGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                <  2:808872/151‑1 (MQ=255)
                                                                                   aaaaTGAAAGGAAAAAAGGAAAATTTCAAAAAAAAGGCGGTAATTAAGAGATTAGGATGGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATg                                                                <  2:1398500/151‑1 (MQ=255)
                                                                                     aaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                               >  1:1477355/1‑150 (MQ=255)
                                                                                          cAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                          >  2:1344314/1‑150 (MQ=255)
                                                                                            aGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                        >  2:440771/1‑151 (MQ=255)
                                                                                                  aaaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTAc                                                   >  2:882610/1‑151 (MQ=255)
                                                                                                    aaGGAAAAATTCAAAAAAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                                 <  2:1163830/151‑1 (MQ=255)
                                                                                                        agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             >  1:57184/1‑151 (MQ=255)
                                                                                                        agaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             >  1:215779/1‑151 (MQ=255)
                                                                                                        aaaaaTTCTAAAAAAATGCGGTAATTTAAAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             <  1:51352/149‑1 (MQ=255)
                                                                                                          aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                           >  1:450311/1‑151 (MQ=255)
                                                                                                          aaaTTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACtt                                           >  2:1458721/1‑151 (MQ=255)
                                                                                                           aatttCTAAATAAATGCGTTAATTTAGAGATTGGGTTTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                            <  2:354123/147‑1 (MQ=255)
                                                                                                              ttCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTcc                                        >  2:1004684/1‑151 (MQ=255)
                                                                                                               tctcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                       >  1:1448685/1‑151 (MQ=255)
                                                                                                                  cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                    >  2:1031734/1‑151 (MQ=255)
                                                                                                                   aaaaaaGGCGTTAACTAAAAGAATGGGATGGCGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTc                                         <  1:1087297/142‑1 (MQ=255)
                                                                                                                   aaaaaaGGCGGTAACTAAAAGATTAGGTTGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                     <  2:1313733/146‑1 (MQ=255)
                                                                                                                   aaTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   >  2:468695/1‑151 (MQ=255)
                                                                                                                   aaTAAAGGCGTTAACTTAGATTTTAGGATTGCGAAAAATACCAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                             <  1:1438888/141‑1 (MQ=255)
                                                                                                                   aaTAAAGGCGGAACCTAAGAGATTGGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGt                                   <  2:207317/151‑1 (MQ=255)
                                                                                                                      aaaGGCGGTAACTAAGAGATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                <  1:1029515/151‑1 (MQ=255)
                                                                                                                        aTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              >  1:1170201/1‑151 (MQ=255)
                                                                                                                          gCGGTAATTAAGAGATTAGGATGGGGGAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                            <  2:763882/151‑1 (MQ=255)
                                                                                                                          gCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              <  2:1313296/149‑1 (MQ=255)
                                                                                                                          gCGGAAATTTAGAGTTTGGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                              <  1:822337/149‑1 (MQ=255)
                                                                                                                                aaCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                       >  2:1091266/1‑151 (MQ=255)
                                                                                                                                   ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    >  1:1227750/1‑151 (MQ=255)
                                                                                                                                   ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    >  2:267111/1‑151 (MQ=255)
                                                                                                                                    tAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                   >  2:413840/1‑151 (MQ=255)
                                                                                                                                     agagATTAGGATGGGGAAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                    <  2:1035158/149‑1 (MQ=255)
                                                                                                                                       agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  2:1167645/1‑150 (MQ=255)
                                                                                                                                       agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  1:871553/1‑150 (MQ=255)
                                                                                                                                        gaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                 >  1:626704/1‑149 (MQ=255)
                                                                                                                                         aTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc              >  2:311060/1‑151 (MQ=255)
                                                                                                                                            gggATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag               <  1:51213/146‑1 (MQ=255)
                                                                                                                                             ggtttGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg          <  1:65453/148‑1 (MQ=255)
                                                                                                                                             ggATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca           <  2:1518704/150‑1 (MQ=255)
                                                                                                                                              gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg         >  2:559315/1‑151 (MQ=255)
                                                                                                                                                  gCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg      <  2:1634589/150‑1 (MQ=255)
                                                                                                                                                    ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg   >  2:384649/1‑151 (MQ=255)
                                                                                                                                                    ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg   >  1:853428/1‑151 (MQ=255)
                                                                                                                                                       agaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  2:1478471/1‑150 (MQ=255)
                                                                                                                                                       agaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGa  >  1:1037129/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                    
ATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGA  >  NC_000913/5013‑5307

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