Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,1590GA41.1% 127.0 / 62.6 108intergenic (+139/‑75)thrC/yaaXL‑threonine synthase/DUF2502 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (55/8);  new base A (0/44);  total (55/53)
Fisher's exact test for biased strand distribution p-value = 3.23e-22
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 109 total aligned reads displayed.

TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5011‑5305
                                                                                                                                                     |                                                                                                                                                  
tcatcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGgaga                                                                                                                                                   >  2:1181463/1‑150 (MQ=255)
 catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                 <  2:644256/151‑1 (MQ=255)
 catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAt                                                                                                                                                 <  2:486785/151‑1 (MQ=255)
    caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                              <  3:322687/151‑1 (MQ=255)
    caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaa                                                                                                                                               <  2:1202671/150‑1 (MQ=255)
     aGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaac                                                                                                                                              >  2:303615/1‑150 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaaca                                                                                                                                          <  2:403251/151‑2 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                          >  2:932567/1‑151 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                          >  2:272888/1‑151 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                          <  2:259824/151‑1 (MQ=255)
        aaaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAcc                                                                                                                                          >  2:1212478/1‑151 (MQ=255)
         aaaTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccg                                                                                                                                         >  1:382060/1‑151 (MQ=255)
           aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgcc                                                                                                                                       >  4:760381/1‑151 (MQ=255)
            tCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAAccgccg                                                                                                                                      >  4:1392690/1‑151 (MQ=255)
                 tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTtctc                                                                                                                                 >  2:110412/1‑151 (MQ=255)
                   aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCAt                                                                                                                               >  3:145343/1‑151 (MQ=255)
                      aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                             >  1:265761/1‑150 (MQ=255)
                      aTCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                             >  3:86154/1‑150 (MQ=255)
                       tctcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCg                                                                                                                             >  3:138542/1‑149 (MQ=255)
                         tcAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGt                                                                                                                          >  4:386378/1‑150 (MQ=255)
                           aaTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAt                                                                                                                       >  2:1438954/1‑151 (MQ=255)
                            aTCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtc                                                                                                                      >  1:419682/1‑151 (MQ=255)
                              cAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAAtctc                                                                                                                    >  4:303516/1‑151 (MQ=255)
                               aGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                   >  4:272645/1‑151 (MQ=255)
                                ggCGGGTTTTGTTTTTATGCAGCCGGTTTTTTTTATAAAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAACTAAGAGATTAGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTcc                                                                                                                   <  4:1539954/150‑1 (MQ=255)
                                ggCCGTTTTTGTTTTTTTGCACCCCGGTTTTTTTAAGAAAAAATTATGAAAAAAAATCAGAGGAAAAAAGGAGAAATCTTAAAAAAGTGCGGTAACTTAGAGATTGGGGTTGGGGAAAAAAACACCCGCCGTTCTCATCGAGTAAtctc                                                                                                                    <  4:806317/149‑1 (MQ=255)
                                ggCCGGGTTTGTTTTTAGGCAGCCCGGTTTTTTTATAAAAAAATTATGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGTTAACTTAGAAATTAGGTTTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAAtct                                                                                                                     <  3:350028/148‑1 (MQ=255)
                                ggCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCg                                                                                                                  >  2:92446/1‑151 (MQ=255)
                                    ggttttttttttACGCAGCCCTGTTTTTTTAAGAAAAAATTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTCAATAAATGCGGTAACTTAGAGATTAGTATTGCGGAAAATACCAACCGCCGTTCTCATCGAGTAATCTCCGGat                                                                                                               <  1:460158/142‑1 (MQ=255)
                                    ggttttGTTTTTATGCAGCCGGGTTTTTTTATAAAAAAATTAGGGAAAAAAATGAAAGGAAAAAGGGAAAAATTCTCAATAAATGCGTTAATTTAGAGATTGGGAGTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGata                                                                                                              <  1:203914/148‑1 (MQ=255)
                                    ggtttgttttttATGCACCCCGGTTTTTTAAGAAAAAAATTAGGAAAAAAATTGACAGGAAAAAAGGAAAATTTCTAAAAAAATGCGGTAACTAAGAAATGAGTTTGGCGAAAAAAACCAACCGCCGTTCTCATCGAGTAATCTCCGGa                                                                                                                <  4:1388587/141‑1 (MQ=255)
                                       tttGTTTTTATGCACCCCGGTTTTTTTATAAAAAAATAAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAATAAAGGCGTTAACTTAGAAATTGGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                            <  3:241904/150‑1 (MQ=255)
                                       tttGTTTTTAGGCGCCCGGTCTTTTTTATAAAAAAATTATGAAAAAAAATAACAGAAAAAAAGGAAAATTTCTAAAAAAATGCGGTAATTTAGAGATTGGGTTTGAGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATc                                                                                                            <  2:1209052/150‑1 (MQ=255)
                                         tGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAc                                                                                                         >  4:546438/1‑151 (MQ=255)
                                          gCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAcc                                                                                                        >  4:540600/1‑151 (MQ=255)
                                            ttttATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGAccc                                                                                                       >  2:200250/1‑150 (MQ=255)
                                                   cAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACggg                                                                                               >  2:259811/1‑151 (MQ=255)
                                                    aGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                              >  2:1647808/1‑151 (MQ=255)
                                                    aGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGc                                                                                              >  3:164381/1‑151 (MQ=255)
                                                       ccGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAAt                                                                                           >  4:411906/1‑151 (MQ=255)
                                                        cGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATg                                                                                          >  4:930165/1‑151 (MQ=255)
                                                            tttttttAAGAAAAAATTAGGGAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAACTTAGAGATTAGGATTGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACgg                                                                                                <  1:108542/141‑1 (MQ=255)
                                                            ttttttaATAAAAAAATTATGGAAAAAAAGGAAAGGAAAAAAGAAAAAATCCAAAATAAGGGCGGTAATTAAGAGATGAGTTTGGGGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATgt                                                                                         <  4:1569612/148‑2 (MQ=255)
                                                            cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                        >  4:1458735/1‑150 (MQ=255)
                                                            cTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                        >  3:363463/1‑150 (MQ=255)
                                                             ttttttATAAAAAAATAAGGAAAAAAAAGAACAGGAAAAAAGGAAAAATTCTAAAAAAAGGCGTTAATTAAGAGATTAGTATGGGGGAAAAAACCACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATa                                                                                       <  2:1661729/150‑1 (MQ=255)
                                                              tttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaa                                                                                      >  4:11486/1‑150 (MQ=255)
                                                               ttttATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaa                                                                                     >  4:1412561/1‑150 (MQ=255)
                                                                tttATAAAAAAATTATGAAAAAAAAGGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGTTAATTAAGAGATTAGGTTGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                    <  2:569897/150‑1 (MQ=255)
                                                                   aTGAAAAATTTATGAAAAAAAATGACAGGAAAAAAGGAAAAATTCTAAATAAATGCGGTAATTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                  <  4:434556/149‑1 (MQ=255)
                                                                   aTGAAAAAATTAGGGAAAAAAAGGACAGGAAAAAAGAAAAATTTCTAAAAAAATGGGGTAACTAAGAGATAAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGa                                                                                 <  4:837196/150‑1 (MQ=255)
                                                                    tGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                               >  4:1634693/1‑151 (MQ=255)
                                                                      aaaataTTAGGGAAAAAAAAGAGGGGAAAAAAGGAAAAATTTAAAAAGAATGTGGAAACTTAGATAAGAGATTGGCGAAGAATAACATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                               <  4:325800/144‑1 (MQ=255)
                                                                      aaaaaaTTAGGAAAAAAAATGACAGGAAAAAAGAAAAAATTCTCAAAAAAGGCGGTAATTTAGAGATTAGGATTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                              <  4:975698/147‑1 (MQ=255)
                                                                      aaaaaaTGATGGAAAAAAATGAGGGGAAAAAGGGAAAAATTCAAAAAAAGTGGGTTAACTAAGAGATTGGGATGGGGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                              <  2:193194/147‑1 (MQ=255)
                                                                       agaaATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                            >  4:618794/1‑151 (MQ=255)
                                                                       aaaaaTTAGGGAAAAAAATGAAAGGAAAAAAGGAAAAATTCTAAAAAAATGCGGTAACTTAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAc                                                                            <  3:208403/149‑1 (MQ=255)
                                                                         aaaTTATGAAAAAAAAGAACAGGAAAAAAGGAAAATTTCAAAAAAAAGGGGTTAACTTAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                          <  4:1316075/151‑1 (MQ=255)
                                                                         aaaTTAGGGAAAAAAATGAAAGGAAAAAAGAAAAAATCCTAAATAAAGGGGTTAATTAAGAGATTGGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAAcc                                                                           <  4:325580/150‑1 (MQ=255)
                                                                                ggAAAAAAAGGAAAGGAAAAAAGGAAAAATTCACAAAAAACGGGGAAATTAAGAGATGGGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtg                                                                         <  2:1302992/145‑1 (MQ=255)
                                                                                  agaAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaaaa                                                                   >  4:1224996/1‑149 (MQ=255)
                                                                                  aaaaaaaaGAGGGGAAAAAAGGAAAAATCCACAACATAGGGGGATATTAAGAGATTGGGGTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGaa                                                                     <  2:697014/145‑1 (MQ=255)
                                                                                  aaaaaaaGGAAAGGAAAAAAGGAAAAATTCAAAAAAAATGCGGTAATTTAGAGTTTGGGTTTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtgtg                                                                       <  2:385817/143‑1 (MQ=255)
                                                                                   aaaaaaaGAAAGGAAAAAAGAAAAAATTTAAAAAAAATGCGAAAAGTAAGAAATAGGGATGGCGAAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGa                                                                <  4:1457029/150‑1 (MQ=255)
                                                                                       aaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                             >  2:289299/1‑150 (MQ=255)
                                                                                       aaTGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGc                                                             >  2:909624/1‑150 (MQ=255)
                                                                                       aaTGAAAGGAAAAAGGGAAAAATTCTAAAAAAATGCGGTAATTTAGAGATTAGGTTTGCGGAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCa                                                            <  3:85205/151‑1 (MQ=255)
                                                                                          gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                         >  4:1474503/1‑151 (MQ=255)
                                                                                          gACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATc                                                         >  4:1437489/1‑151 (MQ=255)
                                                                                          gACAGGAAAAAAGGAAAAATTCTCAAAAAAGGGGGTAACTTAGAGATTAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCaa                                                           <  4:1053426/149‑1 (MQ=255)
                                                                                             aGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                      >  4:445519/1‑151 (MQ=255)
                                                                                              ggaaaaaaGAAAAAATTCTAAAAAAAGGCGGAAATTAAGAGATTAGGAGGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                      <  1:410702/147‑1 (MQ=255)
                                                                                              gaaaaaaaGAAAAAATTCTAAAAAAGGTGGGAAACTAAGAAATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAAt                                                          <  4:1214520/143‑1 (MQ=255)
                                                                                                aaaaaaGGAAAAATTCTAAAAAAGTGCGTTAATTTAGAAATAAGGATGGGGAAAAAAAAAACCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTAt                                                      <  2:21338/147‑1 (MQ=255)
                                                                                                 aaaaaGGAAAAATTCTCAATAAGTGGGGTAATTTAGAGATTAGGTTGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTa                                                  <  4:1182950/151‑1 (MQ=255)
                                                                                                    aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                               >  1:107791/1‑151 (MQ=255)
                                                                                                    aaGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTc                                               >  4:328246/1‑151 (MQ=255)
                                                                                                       gagaAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGc                                             >  2:966350/1‑150 (MQ=255)
                                                                                                        aaaatttcaaaaaaaaTGCGGTAACTAAGAGATTAGGTTGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCAc                                           <  4:333887/138‑1 (MQ=255)
                                                                                                         aaaaTTCTAAAAAAAGGGGGTAACTAAGAGATTAGGATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACt                                          <  4:110592/150‑1 (MQ=255)
                                                                                                              tctcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTccc                                     >  1:87493/1‑151 (MQ=255)
                                                                                                               ctcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                    >  1:109174/1‑151 (MQ=255)
                                                                                                                tcAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                   >  3:206224/1‑151 (MQ=255)
                                                                                                                 cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                  >  3:433998/1‑151 (MQ=255)
                                                                                                                 cAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                   >  4:635937/1‑150 (MQ=255)
                                                                                                                 aaaaaaaggggttAACTAAGAAATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTg                                   <  1:167/138‑1 (MQ=255)
                                                                                                                  ataaagtgtggTAAGTTAGAATTGAGGATGGGGGACAAAAACATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCt                                    <  4:273327/139‑1 (MQ=255)
                                                                                                                     aaaTGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                              >  4:505401/1‑151 (MQ=255)
                                                                                                                      aaTGCGGTAATTAAGAGATTAGGATTGCGGAAAATAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTg                             <  2:199981/151‑1 (MQ=255)
                                                                                                                                cTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTccc                    >  1:255209/1‑150 (MQ=255)
                                                                                                                                 ttAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                  >  4:767390/1‑151 (MQ=255)
                                                                                                                                   aaaaaTTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCAt                  <  4:1183638/145‑1 (MQ=255)
                                                                                                                                     agaTTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc               <  3:383630/150‑1 (MQ=255)
                                                                                                                                         tAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGcac          >  2:1487733/1‑151 (MQ=255)
                                                                                                                                           ggATTGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCa                   <  1:316091/140‑1 (MQ=255)
                                                                                                                                           ggATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAg        <  2:792146/151‑1 (MQ=255)
                                                                                                                                           ggAGTGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgcagca           <  2:759372/148‑1 (MQ=255)
                                                                                                                                                gaggAAAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTg    <  2:545681/148‑1 (MQ=255)
                                                                                                                                                 cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCg  >  4:1660917/1‑151 (MQ=255)
                                                                                                                                                  ggAGAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGc      <  4:1038298/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCG  >  NC_000913/5011‑5305

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: