Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	29,486	0	C	A	41.9%	108.7 / 40.8	63	intergenic (+291/‑165)	dapB/carA	dihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (28/8);  new base A (0/26);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.42e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGA  >  NC_000913/29336‑29631                                                                                                                                                        |                                                                                                                                                  tttttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCcac                                                                                                                                                   >  1:654924/1‑151 (MQ=255)   tttGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCcaca                                                                                                                                                  >  1:708630/1‑150 (MQ=255)    ttgttgTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaa                                                                                                                                                 >  1:234198/1‑150 (MQ=255)    ttgttgTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaa                                                                                                                                                 <  2:502645/150‑1 (MQ=255)      gttgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaaaa                                                                                                                                               <  2:1093650/150‑1 (MQ=255)        tgttTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAaata                                                                                                                                             >  1:284633/1‑150 (MQ=255)                   aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                 >  2:957299/1‑151 (MQ=255)                     aTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATccc                                                                                                                               <  2:254711/151‑1 (MQ=255)                        ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                             <  1:751218/150‑1 (MQ=255)                        ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                             <  1:1051545/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)                                                                                                                                                        |                                                                                                                                                  TTTTTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGA  >  NC_000913/29336‑29631

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