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

CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAA  >  NC_000913/5012‑5308
                                                                                                                                                      |                                                                                                                                                     
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                    <  2:1114512/151‑1 (MQ=255)
catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                    >  2:369277/1‑151 (MQ=255)
   caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                                  <  2:1172898/150‑1 (MQ=255)
      tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                              >  1:931192/1‑151 (MQ=255)
      tAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaa                                                                                                                                               >  1:839441/1‑150 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                             >  1:1124824/1‑151 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                             <  1:660349/151‑1 (MQ=255)
       aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                             <  1:873043/151‑1 (MQ=255)
        aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                            >  1:513566/1‑151 (MQ=255)
         aaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgc                                                                                                                                           >  1:358155/1‑151 (MQ=255)
           tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgccg                                                                                                                                         >  1:395659/1‑151 (MQ=255)
              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                          <  2:749409/147‑1 (MQ=255)
              aTTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                          <  1:1208231/147‑1 (MQ=255)
                tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                    >  1:602010/1‑151 (MQ=255)
                 cATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                    >  2:9450/1‑150 (MQ=255)
                  aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCAt                                                                                                                                  >  1:583292/1‑151 (MQ=255)
                      tctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                                >  1:963231/1‑149 (MQ=255)
                          aaTCGGGCGGGTTTTGTTTTTATGCAGCCCGTTTTTTTTATAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTCAATAAAGGGGTTAATTTAGAGATTAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCg                                                                                                                                <  1:177808/145‑1 (MQ=255)
                           aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                         >  1:475502/1‑151 (MQ=255)
                            tCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                        >  2:484902/1‑151 (MQ=255)
                               ggCGGGGTTTGTTTTTATGCACCCCGGTTTTTTTAAAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGAAAAAATTCTAAAAAAAGGCGGTAACTTAGAAATTAGGATGGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                      <  2:884691/150‑1 (MQ=255)
                                   gggTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                                  >  1:1251786/1‑150 (MQ=255)
                                      tttttttttATGCAGCCCGGTTTTTTTATGAAGAAATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATCCTAAAAAAATGCGGTAACTTAGAGTTTAGGATGGGGGAGAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGatat                                                                                                                <  1:945853/144‑1 (MQ=255)
                                             ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                        >  2:447396/1‑150 (MQ=255)
                                             ttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTCCTCATCGAGTAATCTCCGGATATCGACCCAt                                                                                                        <  1:1202595/150‑1 (MQ=255)
                                              tATGCAGCCCGGTTTTTTTAAAAAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCTCAATAAAGGCGTAAATTTAGAAATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATa                                                                                                       <  1:442460/150‑1 (MQ=255)
                                                     ccgggtttttttATAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAATTTAGAGATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCaa                                                                                               <  1:445899/145‑1 (MQ=255)
                                                        gttttttttATAAAAAAATTAGGGAAAAAAATGACAGGAAAAAAGGAAAATTTCCAAAAAAATGCGTTAATTAAGATATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                              <  2:813006/146‑1 (MQ=255)
                                                        gtttttttaATAAAAAAATTAGGAAAAAAAAAGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAACTAAGAGATTAGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                             <  2:1258533/147‑1 (MQ=255)
                                                          tttttttATAAAAAAATTATGAAAAAAAAGGAAAGGAAAAAAGGAAAAATTCACAAAAAATGCGATAATTTAGAAATTGGGATTGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                            <  2:281580/148‑1 (MQ=255)
                                                           ttttttATAAAAAAATAAGGAAAAAAAAGGACAGGAAAAAAGAAAAATTCCAAAATAAAGGCGGTAACTAAGAAATAAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGa                                                                                            <  2:1127629/148‑1 (MQ=255)
                                                            tttttATAAAAAAATTATGGAAAAAAAGGACAGGAAAAAAGGAGAAATTCTAAAAAAAGGCGTTAACTAAAAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                         <  1:503634/150‑1 (MQ=255)
                                                             ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        >  1:287309/1‑150 (MQ=255)
                                                             ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        >  1:996502/1‑150 (MQ=255)
                                                             tttaaaaaaaaaaTAAGGGAAAAAAATGAAAGGAAAAAAGGAAAATTCTAAAAAAAATGGGTAAATTAAGATTTTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                        <  2:501039/140‑1 (MQ=255)
                                                               ttaaggaagaaTTTTGGAAAAAAAATGCAGGGAAAAAAAGAAAAATTCTTAATAATGGCGTTAATTTAGAGATTAGGATTGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                       <  1:777273/141‑1 (MQ=255)
                                                                 agaaagaaATTATGGAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGGTAATTAAGAAATTAGGATGGGGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                    <  1:550386/147‑1 (MQ=255)
                                                                 aTGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                     >  1:1188611/1‑149 (MQ=255)
                                                                    aggaaATTATGAAAAAAAATGACAGGAAAAAAGAAAAAATCCTAAAAAAATGGGGTAATTTAGAGATTAGGATTGCGAAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                     <  1:1182848/144‑1 (MQ=255)
                                                                    aaaaaaTTAGGAAAAAAAATGACAGGAAAAAAGAAAAAATTCTAAAAAAAGGCGGTAATTAAGAAATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                  <  2:669642/146‑1 (MQ=255)
                                                                      gaaATTATGGCGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                              >  2:462967/1‑151 (MQ=255)
                                                                        aaTTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            >  2:1225397/1‑151 (MQ=255)
                                                                          ttATGAAAAAAAAGGAAAGGAAAAAAGAAAAAATTCTAAAAAAGGGCGTTAATTTAGAGATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgt                                                                           <  1:1047402/150‑1 (MQ=255)
                                                                          taagagaaaaaaaaGACAGAAAAAAAAGAAAAATTCACAAAAAATGCGAAAATTAAGAAATGAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                            <  2:296424/141‑1 (MQ=255)
                                                                            aTGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGa                                                                         >  1:26393/1‑150 (MQ=255)
                                                                                aaaaaaaTGAAAGAGAAAAAGGAAAAATTCTCAAAAATTGCGGTAACTTAGAGATTAGGATTGGGAAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaaa                                                                     <  2:1059263/148‑1 (MQ=255)
                                                                                 aaaaaggaaCAGGGAAAAAGAAAAAACTGCAAATAAGTGCTTTAACTGATAAATTGGGATGGAAGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                   <  2:666011/143‑1 (MQ=255)
                                                                                  aaaaaTGAAAGGAAAAAAGAAAAAATCCAAAAAAAAGGCGTTAATTAAGAGATTGGGTTGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGAt                                                                  <  1:1202237/151‑1 (MQ=255)
                                                                                        gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                            >  1:20249/1‑151 (MQ=255)
                                                                                        gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                            >  1:218612/1‑151 (MQ=255)
                                                                                           agggaaaaagggaaaaTTTCACAAAAAGTGCGATAACTTAGAATTTGGGATGGCAGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                          <  2:439082/151‑1 (MQ=255)
                                                                                             ggaaaaaaGGAAAAATTCTCAAAAAATGGGGAAATTTAGAGATTGGGATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                         <  1:272640/147‑1 (MQ=255)
                                                                                             ggaaaaaaGGAAAAATTCTAAAAAAAGGGGGTAATTTAGAGATTAGGATTGCGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCt                                                           <  1:861882/145‑1 (MQ=255)
                                                                                             ggaaaaaaGGAAAAATTCAAAATAAATGCGTTAACTGAGATATTAGGTTTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTa                                                          <  2:279736/146‑1 (MQ=255)
                                                                                               aaaaaaGGAAAAATTCTAAAAAAAGGGGGTAACTTAGAAATTAGGTTTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCg                                                       <  1:1235385/149‑1 (MQ=255)
                                                                                                aaaaaGGAAAAATTCTCAAAAAGGGCTTTAACTAGGATATTGGGATGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATAAATCg                                                       <  1:132352/149‑1 (MQ=255)
                                                                                                aaaaaGAAAAAATTCAAAAAAAATGCGTAAACTTAGAGATGAGGATTGAGAAAAAAAACACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATc                                                        <  2:1124810/148‑1 (MQ=255)
                                                                                                  aaaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                   >  1:6970/1‑151 (MQ=255)
                                                                                                  aaaGGAAAAATTCAAAATAAATGGGGAAACTAAGAATTGGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACt                                                   <  1:933272/151‑1 (MQ=255)
                                                                                                        aaaaaTTCTAAAAAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                              <  1:532156/149‑1 (MQ=255)
                                                                                                         gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             >  1:167648/1‑151 (MQ=255)
                                                                                                         gaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             >  2:1120757/1‑151 (MQ=255)
                                                                                                          aaaTTCTAAAAAAATGCGTTAACTTAAAGATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                             <  1:708731/150‑1 (MQ=255)
                                                                                                                 cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                     >  1:1234778/1‑151 (MQ=255)
                                                                                                                  aaaaaagggggTAATTAAGAGTTTAGGATTGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                     <  1:899924/141‑1 (MQ=255)
                                                                                                                   aaaaaggggtaaatttAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                   <  1:507967/137‑1 (MQ=255)
                                                                                                                   aTAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGtt                                   >  1:329075/1‑151 (MQ=255)
                                                                                                                    tAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTc                                  >  1:1182700/1‑151 (MQ=255)
                                                                                                                              aaCTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTcc                        >  1:860912/1‑151 (MQ=255)
                                                                                                                                 ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                     >  1:66116/1‑151 (MQ=255)
                                                                                                                                      agaTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                  >  1:224891/1‑150 (MQ=255)
                                                                                                                                        aTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcag                >  1:386651/1‑150 (MQ=255)
                                                                                                                                         ttAGGTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               <  1:129822/150‑1 (MQ=255)
                                                                                                                                          tAGAATTGAGGAAAATACAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca              <  1:1031244/150‑1 (MQ=255)
                                                                                                                                            ggATGGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcaca            <  2:1062823/150‑1 (MQ=255)
                                                                                                                                             gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg          >  2:1021965/1‑151 (MQ=255)
                                                                                                                                             gATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAgg          >  1:753655/1‑151 (MQ=255)
                                                                                                                                               tgggggAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagc               <  1:786581/140‑1 (MQ=255)
                                                                                                                                                 gggaagaAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGATGc      <  2:358084/147‑1 (MQ=255)
                                                                                                                                                 gCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc      <  1:344479/151‑1 (MQ=255)
                                                                                                                                                   ggAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg    >  1:791935/1‑151 (MQ=255)
                                                                                                                                                   ggAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg       <  1:1168713/148‑1 (MQ=255)
                                                                                                                                                    gagaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg    >  1:118035/1‑150 (MQ=255)
                                                                                                                                                    gagaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg     >  1:908933/1‑149 (MQ=255)
                                                                                                                                                    gagaAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGc      <  1:43592/148‑1 (MQ=255)
                                                                                                                                                      gaATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGaa  >  2:1235683/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAA  >  NC_000913/5012‑5308

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