Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,496	0	T	A	44.3%	93.6 / 38.5	70	intergenic (+301/‑155)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (32/7);  new base A (0/31);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 2.26e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGA  >  NC_000913/29354‑29643                                                                                                                                                |                                                                                                                                                    aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                             >  2:957299/1‑151 (MQ=255)   aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                                           <  2:254711/151‑1 (MQ=255)      ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                         <  1:1051545/150‑1 (MQ=255)      ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                         <  1:751218/150‑1 (MQ=255)         aCCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                       <  2:79369/149‑1 (MQ=255)          ccttttGGTCCATTTTTTTCTGCTCGTTTTTATTTCATCCAATCTTTTTCCTGCGCAAGCTTTTTCA‑AGAAAGGTTTAGATAATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCAt                                                                                                                                      <  2:327239/146‑1 (MQ=255)            aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                    >  1:878845/1‑149 (MQ=255)              ttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGtt                                                                                                                                >  1:90830/1‑151 (MQ=255)                    cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                          >  2:842331/1‑151 (MQ=255)                       ttttttCTGCTCTTTTTTATTTCATCCAATCTTCTTGCGGCGCAAGCTTTTCCC‑AGAAAAGGTTAGATGATTTTTTTGTCTTTAATTGCCTGAAAAACTTGCATGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                          <  1:21/148‑1 (MQ=255)                       ttttttCTCCTTGTTTTTATTAAACAAATTTTTTTTCCGGCGCAGGCGTTTTCA‑AAAAGAGTTTAGATAATTTTTTTGTCATTTAATCCCTAAAAAACTGGTATGAGCCAAAAAAAAATATAAAAAACCCCGCCATTAAGTTGACTTTTa                                                                                                                        <  2:557095/150‑1 (MQ=255)                                cTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                              >  1:755751/1‑151 (MQ=255)                                 tctttttttttTCATGCATTCTTCTTGCTGCGCAAGCTTTTCCA‑AGAACAGTTTAGATGTTTTTTTTGTCGTTTATTGCTTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                                <  1:1039340/145‑1 (MQ=255)                                   tttttttttTCATCCATTCTTCTTGCTGCGCAAGCGTTTCCC‑AAAACGGGTTAGATGTTTTTTTTGTCCTTTAATGCCTGAAAAACTTGCATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                              <  1:7013/147‑1 (MQ=255)                                    ttttttttTCATGCATTCTTCTTGCTGCGAAAGCTTTTCCC‑AGAACAGGTTAGATGATTTTTTTGTCTTTAAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                             <  1:808675/148‑1 (MQ=255)                                    tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                          >  1:570909/1‑151 (MQ=255)                                     ttttttttCATGCAATCTTCTTGCTGCGAAAGCTTTTTCC‑AAAACAGTTTAGATGATTTTTTTGTCGCTTAATGCCTAAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                          <  1:524736/150‑1 (MQ=255)                                      ttttttTCATGCATTTTTTTTGCTGCGCAAGCTTTTCCC‑AAAACAGTTTAAATGATTTTTTTGTCTTTTAATCCCTATAAAACTTGCATGCGCCAAAAAATAATATAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                         <  2:595847/150‑1 (MQ=255)                                      ttttttTCATGCAATCTTCTTGCTGCGCAAGCTTTTCCC‑AGAACGGTTTAGATGTTTTTTTTGCCGCTTAATGCTTGAAAAACATGCATGAGCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                                         <  1:363523/150‑1 (MQ=255)                                        tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                       >  1:725993/1‑150 (MQ=255)                                          tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                    >  2:759102/1‑151 (MQ=255)                                          tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                    >  1:881563/1‑151 (MQ=255)                                           ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                   >  1:774680/1‑151 (MQ=255)                                             cATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                                 >  2:671271/1‑151 (MQ=255)                                                  attcttctTGCCGCGCAAGCTTTTCCA‑AAAACAGTTTAGATTTTTTTTTTGTCGTTTAATGCTTAAAAAACAGTCTTGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                    <  2:1012183/141‑1 (MQ=255)                                                      ttcttGCTG‑GCAAAGTTTTTCCAAAAACAGGTTAGTTGTTTTTTTTGTCGTTTAATGCTTAAAAAACATGCATGACCCAAAAAAAAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                         <  1:969199/150‑1 (MQ=255)                                                          tgctgcGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                     >  1:858786/1‑150 (MQ=255)                                                          tgctgcGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGa                                                                                     >  1:147534/1‑150 (MQ=255)                                                            ctgcGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                  >  1:1069214/1‑151 (MQ=255)                                                                gcAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                              >  1:720491/1‑151 (MQ=255)                                                                  aaGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                            >  1:3471/1‑151 (MQ=255)                                                                       ttttCC‑AGAACAGGTTAGATGATTTTTTTGTCGTTAAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                           <  1:546454/147‑1 (MQ=255)                                                                        ttcac‑AGAATAGGTTAGATGATTTTTTTGTCGTTTAAGGCCTAAAAAACATGCATGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                       <  2:701610/146‑1 (MQ=255)                                                                          tCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                      >  1:832028/1‑149 (MQ=255)                                                                            caaaaaaggtttAGATTTTTTTTTTTTCTTTTAACCCCTAAAAAACATGAAGGCCCAAAAAAAAAAAAAAAAAACCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                   <  2:423659/141‑1 (MQ=255)                                                                             aaaaatggtttATATTATCTTTTTGTCGCTTAATCCCTAAAAAACATTCAGGCCCCAAAAAAAAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                                     <  2:302074/139‑1 (MQ=255)                                                                                aaCAGGTTAGATGTTTTTTTTTCCGTTTAATTCCTGTAAAACTTGCAGGCCCCACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                                <  1:265045/150‑1 (MQ=255)                                                                                    ggTTAGATGATTTTTTTGTCGTTAAATGCTTGAAAAACATGAAGGCCCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCa                                                             <  1:1018542/149‑1 (MQ=255)                                                                                        aGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                       >  1:353581/1‑151 (MQ=255)                                                                                           ttttttttttATCGCTTAATCCCTAAAAAACATGCTTGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAt                                                        <  2:566011/142‑1 (MQ=255)                                                                                                 ttttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                              >  1:814824/1‑151 (MQ=255)                                                                                                  tttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                              >  2:344120/1‑150 (MQ=255)                                                                                                     gTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                          >  1:118142/1‑151 (MQ=255)                                                                                                        gCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                        >  2:532322/1‑150 (MQ=255)                                                                                                          ttaatccataaaaaaCATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                                            <  2:856838/133‑1 (MQ=255)                                                                                                          ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                     >  1:746383/1‑151 (MQ=255)                                                                                                          ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                     >  2:824471/1‑151 (MQ=255)                                                                                                          ttAATCCCTAAAAAACATTAATGAGCCAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                        <  2:497634/148‑1 (MQ=255)                                                                                                             agcctttaaaaaCATGTATG‑GCCA‑AAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                                  <  2:5980/141‑1 (MQ=255)                                                                                                              tGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                                 >  1:289262/1‑151 (MQ=255)                                                                                                                    aaaaaCATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                            <  1:1057210/149‑1 (MQ=255)                                                                                                                    aaaaaCATCAATGACCAAAAAAATAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                                <  2:708608/145‑1 (MQ=255)                                                                                                                     aaaaCAGGCAGGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                           <  1:260887/150‑1 (MQ=255)                                                                                                                      aaaggttggcgcgCCAAAAAAAAAAATAAAAACGCCCGCCATTAAGTTTACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCATAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                         <  2:90828/139‑1 (MQ=255)                                                                                                                       aacgggcttgACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                        <  1:842356/143‑1 (MQ=255)                                                                                                                       aacagacatgACCAAAAAAAAAATCAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                             <  1:449166/140‑1 (MQ=255)                                                                                                                             catgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAATGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                   >  2:712927/1‑150 (MQ=255)                                                                                                                                gccccAAAAAAAAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg               <  2:878819/148‑1 (MQ=255)                                                                                                                                 acccAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga              <  1:543041/149‑1 (MQ=255)                                                                                                                                       aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc        >  1:119983/1‑151 (MQ=255)                                                                                                                                       aaaaTAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc        <  1:1020796/151‑1 (MQ=255)                                                                                                                                       aaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc        <  1:1066808/151‑1 (MQ=255)                                                                                                                                        aaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct       <  1:671294/151‑1 (MQ=255)                                                                                                                                        aaaTAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct       <  1:491410/151‑1 (MQ=255)                                                                                                                                          ataataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct     >  2:452293/1‑151 (MQ=255)                                                                                                                                           taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg    >  1:409867/1‑151 (MQ=255)                                                                                                                                           aaataAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg    <  2:722052/150‑1 (MQ=255)                                                                                                                                            aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg   >  1:842774/1‑151 (MQ=255)                                                                                                                                            aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg   >  1:765065/1‑151 (MQ=255)                                                                                                                                            aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg   >  1:442485/1‑151 (MQ=255)                                                                                                                                             atatAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGa  <  2:687839/151‑1 (MQ=255)                                                                                                                                              tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg   >  2:231173/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    AAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGGA  >  NC_000913/29354‑29643

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