Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA43.8% 92.4 / 30.6 64intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (32/4);  new base A (0/28);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 6.43e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29348‑29646
                                                                                                                                                       |                                                                                                                                                      
tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                       >  1:717978/1‑150 (MQ=255)
tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                       <  1:43234/150‑1 (MQ=255)
 aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                     <  2:478223/151‑1 (MQ=255)
 aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaa                                                                                                                                                      >  2:996877/1‑150 (MQ=255)
   tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                                    >  2:52199/1‑150 (MQ=255)
        aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                              >  1:978370/1‑151 (MQ=255)
          tttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                            >  1:637553/1‑151 (MQ=255)
             gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                          <  1:49587/150‑1 (MQ=255)
                 aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                       >  2:747785/1‑149 (MQ=255)
                 aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                       >  2:252258/1‑149 (MQ=255)
                     ggcccATTTTTTTCTGCTCTTTTTTTTTTCATGCAATCTTCTTGCTGCCAAAGCGTTTCCAAAAACAGGTTAAATAATTTTTTTGTCGTTTATTCCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTg                                                                                                                                  <  1:792513/147‑1 (MQ=255)
                      gTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc                                                                                                                                >  2:929153/1‑151 (MQ=255)
                            ttttttCTGCTTTTTTTTTTTTAATGCATTCTTTTTGCTGCGAAGGCTTTTTCAAAAACAGTTAAAATAATTTTTTTGTCTCTAAATCCCTGAAAAACATGAATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                           <  1:920841/150‑1 (MQ=255)
                             tttttCTGCTTTTTTTTATTTCATGCAATTTTCTTGCTGCGAAAGCGTTTTCCAGAAAAGGTTAGATGATCTTTTGGCCCTTAATTGCTTAAAAAACATGCATGAGCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAg                                                                                                                          <  2:332094/150‑1 (MQ=255)
                                ttCTGCTCTTTTTTATTTCATACAATTTTCTTGCTGCGCAAGCTTTTTCCAAAACAGGTTAAATGATTTTTTTTTCTCTTATTCCCTGAAAAACTTGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                                       <  1:459730/150‑1 (MQ=255)
                                     ctctTTTTTTTTTCATGCAATTTTTTTGCTGCGCAAGCTTTTTCCAGAACAGTTTAGATAATTTTTTTGTCGTTTAATCCCTTAAAAACTTGCATGACCCAAAAAAAAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCa                                                                                                                    <  1:220267/148‑1 (MQ=255)
                                     ctctTTTTTATTTCATGCAATCTTCTTGCTGCGCAGGCTTTTTCCAGAACAGTTTAGATGATTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                   <  2:479260/149‑1 (MQ=255)
                                      tcttttttATTTCATGCAATCTTCTTGCTGCGCAAGCTTTTCCCAAAACAGGTTAGATGTTTTTTTTGTCGCTTAATGCCTAAAAAACATGCATGAGCCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                     <  1:249684/143‑1 (MQ=255)
                                         tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                             >  1:836462/1‑151 (MQ=255)
                                             tATTTCATGCATTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGTTTAATGCCTGTAAAACATGCAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                          <  2:340949/150‑1 (MQ=255)
                                               tttCATGCAATCTTCTTGCTGCGCAAGCTTTTCCCAGAACAGTTTAGATGATTTTTTTGTCGCTAAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                        <  2:841048/150‑1 (MQ=255)
                                               tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                       >  1:599960/1‑151 (MQ=255)
                                                   aTGCATTCTTTTTGCTGCGAAAGCTTTTTCCAAAAAGGGTAAGATGATTTTTTTGCCGTTTAATGCTTAAAAAACATGCATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                                   <  2:495515/151‑1 (MQ=255)
                                                         tcttctTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTg                                                                                             >  1:619660/1‑151 (MQ=255)
                                                            tcttGCTGCGCAAGCGTTTTACAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  1:891958/1‑150 (MQ=255)
                                                             cttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                           >  1:713436/1‑149 (MQ=255)
                                                              tggctacaaaaGTGTTTCCCAAAAAAGTTAAGAGTATCTTTTTCTCTTTAAAGCCTTAAAAACCATGCTGGACCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                               <  1:137987/135‑1 (MQ=255)
                                                                 ctgcGAAGGCTTTTCCAAAAACAGGTAAGATGATTTTTTTGTCGCTTAATGCCTTAAAAACATGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                                          <  1:122033/146‑1 (MQ=255)
                                                                  tgcGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                    >  2:446805/1‑151 (MQ=255)
                                                                    cgcAAGCGTTTTCCAGAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                                   <  2:344887/150‑1 (MQ=255)
                                                                      cAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                                >  2:681459/1‑151 (MQ=255)
                                                                         gcttttcCCAAAACGGTTTAGATGTTTTTTTTGTCTTTAAATGCCTAAAAAACATGCATGACCAAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                                  <  2:394052/143‑1 (MQ=255)
                                                                           gTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                            >  2:257218/1‑150 (MQ=255)
                                                                            ttttCCAAAACAGGTTAGATAATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                             <  1:252192/148‑1 (MQ=255)
                                                                            ttttCAAGAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                              <  2:1011374/147‑1 (MQ=255)
                                                                            tttcCCAAAAAAGGTAAGATAATTTTTTTGTCGTTTAATCCCTGAAAAAATTGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                                <  1:4996/145‑1 (MQ=255)
                                                                                ccAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                         >  1:885216/1‑148 (MQ=255)
                                                                                  aGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                                    >  1:663724/1‑151 (MQ=255)
                                                                                      cAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGc                                                                 >  2:598605/1‑150 (MQ=255)
                                                                                                gatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                      >  2:907337/1‑151 (MQ=255)
                                                                                                gatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                      >  2:455288/1‑151 (MQ=255)
                                                                                                  tttttttGTCTTTTAAGCCCTAAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                     <  1:511921/148‑1 (MQ=255)
                                                                                                    tttttGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                                   <  2:458062/150‑1 (MQ=255)
                                                                                                      tttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                                 >  1:50236/1‑150 (MQ=255)
                                                                                                         gTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                             >  2:371126/1‑151 (MQ=255)
                                                                                                          tCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                            >  1:992775/1‑151 (MQ=255)
                                                                                                           cgttaaaTGCTTGAAAAACATGCATGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                              <  2:103003/143‑1 (MQ=255)
                                                                                                           cGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                           >  2:88176/1‑151 (MQ=255)
                                                                                                            gCTTAATGCCTGAAAAACATGCATGACCCAAAAAATATTAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                           <  2:230868/150‑1 (MQ=255)
                                                                                                                taTGCCTTTAAAAACTTGCAACAGCCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                        <  2:718099/149‑1 (MQ=255)
                                                                                                                 aTGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                      >  2:995954/1‑150 (MQ=255)
                                                                                                                      ccttaaaaaCATGAATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                        <  1:996715/140‑1 (MQ=255)
                                                                                                                           aaaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                              >  1:481151/1‑150 (MQ=255)
                                                                                                                           aaaaCAGGCAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                                 <  2:978532/147‑1 (MQ=255)
                                                                                                                           aaaaCAGGAAGGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                              <  1:952295/150‑1 (MQ=255)
                                                                                                                            aaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                            >  1:758675/1‑151 (MQ=255)
                                                                                                                              aCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                           >  2:68310/1‑150 (MQ=255)
                                                                                                                                     atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                      >  2:505115/1‑149 (MQ=255)
                                                                                                                                      tgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt                   >  1:239492/1‑151 (MQ=255)
                                                                                                                                       gAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg                  >  1:265524/1‑151 (MQ=255)
                                                                                                                                       gAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg                  <  2:637663/151‑1 (MQ=255)
                                                                                                                                              aaaaaaaaaaaaaaaaTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc           <  1:123309/151‑1 (MQ=255)
                                                                                                                                              aaaaaaATAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat             <  2:775225/149‑1 (MQ=255)
                                                                                                                                              aaaaTAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc           <  1:257150/151‑1 (MQ=255)
                                                                                                                                               aaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct          >  1:731521/1‑151 (MQ=255)
                                                                                                                                                   aataAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      <  2:151938/151‑1 (MQ=255)
                                                                                                                                                   aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      <  1:521306/151‑1 (MQ=255)
                                                                                                                                                     tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg      >  2:1034644/1‑149 (MQ=255)
                                                                                                                                                       taAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAggg  >  2:779555/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGAGGG  >  NC_000913/29348‑29646

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