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

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642                                                                                                                                                        |                                                                                                                                                   ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                    >  1:64269/1‑150 (MQ=255)      gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCCCAAAATAATATaaaaa                                                                                                                                               >  2:828547/1‑150 (MQ=255)      gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               >  2:213217/1‑150 (MQ=255)         aaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATcc                                                                                                                                           >  1:363868/1‑151 (MQ=255)               gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                      <  1:1002325/150‑1 (MQ=255)                  cATTTGGTTCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  2:852717/1‑150 (MQ=255)                  cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  1:803893/1‑150 (MQ=255)                  cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  1:1156628/1‑150 (MQ=255)                    tttGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGt                                                                                                                                >  1:263878/1‑151 (MQ=255)                           cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                         >  1:142293/1‑151 (MQ=255)                             cTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        >  2:785833/1‑150 (MQ=255)                              ttttttCTGCTCTTTTTTTTTTCATGAATTTTTTTTGCGGCGCAAGCTTTTCCA‑AGAACAGGTTAGATGATTTTTTTGTCGTTAAATGCTT‑GAAAAACATGAATGGCCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACttt                                                                                                                         <  1:672324/148‑1 (MQ=255)                               tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                     >  2:809175/1‑151 (MQ=255)                                 tttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcgc                                                                                                                   >  1:435799/1‑151 (MQ=255)                                       ttCGTTTTTATGTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATCATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCGCata                                                                                                              >  2:1152340/2‑150 (MQ=255)                                         cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                            >  2:360158/1‑150 (MQ=255)                                           tttttATTTCATGAAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:532862/1‑151 (MQ=255)                                            ttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:901726/1‑150 (MQ=255)                                              ttttttCATGCATTTTTCTTGCTGCGCAAGTTTTTTCA‑AAAACAGTTTAGATAATTTTTTTGTCGCTTAATCCTT‑AAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                       <  2:1029917/147‑1 (MQ=255)                                                 ttccaagaactcttgtTGCTCCGAACGTTTTTCCC‑AGAAAAGTTTAGATATTCTTTTTTCCGATTAGTGCCT‑GAAAACCTTGTATGAACCAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                    <  1:21182/140‑1 (MQ=255)                                                     aTGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                               >  2:1146962/1‑151 (MQ=255)                                                       gCAATCTTTTTGCTGCACAACCGTTTTCC‑AGAACAGGTTCGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCCTGAGCCACAAAATAATATAAACAATCCCGACATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                              >  2:747342/1‑150 (MQ=255)                                                          attcttcttaccgacAAGTGTTTTCCAAAAACGGTTTAGATTATTTTTTTGACGTATAATCCTT‑GTAAAAAATGATTGACCCAAAAAAAAATATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgc                                                                                               <  2:131810/133‑1 (MQ=255)                                                          aTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                          >  2:477648/1‑151 (MQ=255)                                                          aTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                          >  1:390581/1‑151 (MQ=255)                                                            cttcttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGc                                                                                        >  2:656035/1‑151 (MQ=255)                                                                    tgcACAAGCTTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGCTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAtt                                                                                >  2:823243/1‑151 (MQ=255)                                                                         agggtttttCA‑AGAACAGTTTAGATGATTTTTTTGTCGTTAAATGCTT‑GAAAACCATGCAGGGGCCAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtc                                                                               <  1:1052160/142‑1 (MQ=255)                                                                          aggcttttCCAAGAACAGGTTAGATTTTTTTTTTGTCGTTTAATGCCT‑GAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                                     <  1:1143823/137‑1 (MQ=255)                                                                            ctttttCA‑AGAAAAGTATAGTTAATTTTTTTGCCTC‑TAATTCCTAGTAAAAAAAGCATGACCCAGAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           <  1:466838/146‑1 (MQ=255)                                                                             ttttcCC‑AAAACAGTTTATATGTTTTTTTTGTCTTTTACTGCCT‑AAAAAACATGAATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                          <  2:1114669/147‑1 (MQ=255)                                                                              ttttCCAAGAACAGGTTAGATATTTTTTTTTCCGTTAAATCCTT‑GAAAAACATGCAGGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                            <  1:981453/146‑1 (MQ=255)                                                                               tttccaaaaaaaGGTTAGATGTTTTTTTTGTCGTTTAATCCCT‑GAAAAACATGCAGGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                          <  1:146033/136‑1 (MQ=255)                                                                               tttCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  1:471824/1‑151 (MQ=255)                                                                                     agactaGGTTAGTTTTTTTTTTTGTGGCTAAATGCAA‑ATAAAACATGCAGAAACAAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                   <  2:843094/143‑1 (MQ=255)                                                                                       aaCAGGTTAGATGTTTTTTTTGCCATTTAATGCTT‑GAAAAACTTGCAGGCGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAg                                                                   <  1:815653/146‑1 (MQ=255)                                                                                       aaCAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                              >  1:1199016/1‑151 (MQ=255)                                                                                       aaCAGGTAAGATGTTTTTTTTGTCGTTTAATCCTT‑GAAAAACATGAATGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                               <  1:971245/150‑1 (MQ=255)                                                                                             ttaattgatTTTTTTGTCGTTAAAGCCCT‑GAAAAACTTCAATGACCCAAAAAAAAAAAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                          <  1:654300/144‑1 (MQ=255)                                                                                                     tttttttGTCTTTTAATGCCT‑AAAAAACATGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                   <  2:690311/146‑1 (MQ=255)                                                                                                     tttttttGTCGTTTAATGCCT‑GAAAAACAGGAATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                   <  1:109437/146‑1 (MQ=255)                                                                                                        ttttGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                             >  2:1133170/1‑151 (MQ=255)                                                                                                         tttGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatc                                             >  1:134814/1‑150 (MQ=255)                                                                                                              cctttAATCCCT‑AAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                          <  1:760459/145‑1 (MQ=255)                                                                                                               gCTTATTGCCT‑GTAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                      <  2:171351/151‑1 (MQ=255)                                                                                                                   aaTGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                                   >  1:896591/1‑150 (MQ=255)                                                                                                                             aaaaCATGAATCACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGA‑TTGCATTGATTTACGTCATCATTGTGAAt                                    <  2:1093292/139‑1 (MQ=255)                                                                                                                             aaaaCAGGCATGACCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                          <  1:136852/150‑1 (MQ=255)                                                                                                                              agacttctttgAGCCAAAAAAAAATATAAAACGTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        <  2:93434/142‑1 (MQ=255)                                                                                                                                  atgcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                      >  1:484216/1‑149 (MQ=255)                                                                                                                                   tgcttaacccAAAAAATAAAAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                     <  2:195302/141‑1 (MQ=255)                                                                                                                                   tgcatgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                   >  2:987967/1‑151 (MQ=255)                                                                                                                                    gcaggAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                    <  1:1045186/149‑1 (MQ=255)                                                                                                                                      atgAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt                   <  1:61647/148‑1 (MQ=255)                                                                                                                                         aGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg              >  2:19723/1‑150 (MQ=255)                                                                                                                                            cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           >  2:1073559/1‑150 (MQ=255)                                                                                                                                            cacaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           >  1:1142303/1‑150 (MQ=255)                                                                                                                                             acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         >  2:1141963/1‑151 (MQ=255)                                                                                                                                               aaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa            <  1:934290/146‑1 (MQ=255)                                                                                                                                                    aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  <  1:882091/151‑1 (MQ=255)                                                                                                                                                     atatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  1:726037/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                   TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCT‑GTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642

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