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

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642                                                                                                                                                        |                                                                                                                                                   ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                    <  2:338517/150‑1 (MQ=255)    aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 >  2:21900/1‑150 (MQ=255)      gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                              <  1:1033353/151‑1 (MQ=255)             ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:517600/1‑151 (MQ=255)             ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGc                                                                                                                                        <  2:1055093/150‑1 (MQ=255)               gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATAT‑AAAAATCCCGCCAtt                                                                                                                                    >  2:376034/1‑151 (MQ=255)                  cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  2:910009/1‑150 (MQ=255)                     ttGGTCCACTTTTTTCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGAAAGCTTTTTCA‑AAAAAAGTTT‑AGATAATTTTTTTGTCTCTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGt                                                                                                                                <  1:490784/150‑1 (MQ=255)                        gTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAc                                                                                                                            >  2:1060614/1‑151 (MQ=255)                            aCTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        >  2:572305/1‑151 (MQ=255)                             tttttttttCCTCTTTTTTTTTTCATACATTTTTTTTCTTGCGAAAGCTTTTTCA‑AAAACAGTTT‑AAATAATTTTTTTTTCGTTTATCCCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACttt                                                                                                                         <  1:440428/148‑1 (MQ=255)                              ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                        >  2:633545/1‑149 (MQ=255)                              ttttttCTCCTCTTTTTTATTTCATGCATTCTTCTTGCTCCACAAG‑GTTTTCCAAAAAAAGGTA‑AAATATTTTTTTTGCCGTTAAATGCCTATAAAACATGAATGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACttt                                                                                                                         <  1:390104/148‑1 (MQ=255)                               tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                     >  1:922935/1‑151 (MQ=255)                               tttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgc                                                                                                                     >  2:362744/1‑151 (MQ=255)                                          gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                          >  2:11638/1‑151 (MQ=255)                                             tttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  2:723363/1‑149 (MQ=255)                                              ttATTTCATGCATTCTTCTTGCTGCGCAAGCGTTTTCC‑AAAACAGGTT‑AGATTTTTTTTTTGTCTTTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                       <  2:838981/150‑1 (MQ=255)                                                 tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                   >  1:796466/1‑151 (MQ=255)                                                  ttCATGCATTTTTCTGGCGGCGAAAGCTTTTTCC‑AAAAAGGTTT‑AGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATGCAGGACCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                                   <  2:1022460/150‑1 (MQ=255)                                                  ttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgc                                                                                                  >  1:839450/1‑151 (MQ=255)                                                         attcttctGGCTGCGCAAGCTTTTCCA‑AAAACAGTTT‑AGATGTTTTTTTTGTCGTTTAATCCCTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                              <  2:910343/146‑1 (MQ=255)                                                         aaTCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                           >  2:766358/1‑151 (MQ=255)                                                           tcttctTGCTGCGCAAGCTTTTCCA‑AAAACGGTTT‑AGATGATTTTTTTGTCTTTAAACGCCTAAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                              <  1:429766/146‑1 (MQ=255)                                                              tcttGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  1:244064/1‑150 (MQ=255)                                                                tttctgcGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  1:419355/4‑148 (MQ=255)                                                                   ctgccaaaGCTTTTCCA‑AAAACAGTTT‑AAATTTTTTTTTTGCCTTTTAATCCCTGAAAAACATTCTTGGGCCAGAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaaa                                                                                  <  2:800153/144‑1 (MQ=255)                                                                      cgcaattttttcCC‑AGAACGGTTT‑AAATATTTTTTTTTTCCTTTAATGCTTGAAAAACATGCATGACCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                 <  1:643798/140‑1 (MQ=255)                                                                       gaaaGCTTTTTCA‑AAAACTGTTT‑AGATTTTTTTTTTGTCGTTTAATGCCTTAAAACCATGCTTGCGCCAAAAAATAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                                   <  1:434416/143‑1 (MQ=255)                                                                       ccacgtttttcCC‑ATAACAGATA‑AGATCATTTTTTTGTCTTTTATGCCCTGTAAAACATTCAGGCGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcg                                                                              <  2:820398/143‑1 (MQ=255)                                                                         aagttttttCC‑AAAACAGTTT‑AGATTTTTTTTTTTTCGTTTATGGCTTAAAAAACATGCAGGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgt                                                                             <  1:1192249/144‑1 (MQ=255)                                                                         aaGCTTTTTCA‑AAAACAGTTA‑AGATGATCTTTTTGTCGCTTAATGCCTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                            <  2:1144276/150‑1 (MQ=255)                                                                         aaGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           >  1:87627/1‑151 (MQ=255)                                                                         aaGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                           >  2:640616/1‑151 (MQ=255)                                                                          aagcttttCC‑AAAACAAGTTAAAATTATTTTTTTGCTGTTAAATCCCTGTAAAACATGCTGGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                                 <  1:393955/141‑1 (MQ=255)                                                                            cGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                         >  1:778977/1‑150 (MQ=255)                                                                              ttttCA‑AAAACGGTTT‑AGATATTTTTTTTTTCGTTTAATCCCTTAAAAACATGCATGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                        <  1:1025125/149‑1 (MQ=255)                                                                              tttcCC‑AAAACAGGTT‑AGATGTTTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                          <  2:318361/147‑1 (MQ=255)                                                                                 tccaaGAACGGGTT‑AGATTTTTTTTTTGTCTTTTAATGCTGAAAAAACATGAAGGACCCACAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     <  2:1099122/146‑1 (MQ=255)                                                                                 tCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     >  2:293041/1‑149 (MQ=255)                                                                                     aGAACGGTTT‑AGATGATTTTTTTGTCTCTTAAGGCCTGAAAAACATGCAGGAGCAAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAt                                                                 <  1:485036/150‑1 (MQ=255)                                                                                       aaCAGGTT‑AGATTTTTTTTTTTTCTTTTAATGCTTGAAAAAAATGCATGACCAAAAAAATAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCGGAttt                                                               <  2:761368/150‑1 (MQ=255)                                                                                             tt‑AAATGTTTTTTTTGTCGTTTAATGCCTGAAAAACATACATGACCAAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                         <  1:922515/150‑1 (MQ=255)                                                                                                     atCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                                 >  2:345387/1‑151 (MQ=255)                                                                                                      tttttttGTCTTTTAATCCTTTAAAAACATGCATGCCCCAAAAAAAAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                   <  1:751423/146‑1 (MQ=255)                                                                                                       cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                >  1:143755/1‑150 (MQ=255)                                                                                                           ttGTCGGTTATTGCCTAAAAACCATGCAAGAACCAAAAAAGAAAATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                            <  2:168737/150‑1 (MQ=255)                                                                                                               cgtttAATGCCTAAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                        <  2:67023/147‑1 (MQ=255)                                                                                                               cgttaaatccctaaaaaaCATGCATGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                                        <  1:21899/136‑1 (MQ=255)                                                                                                               cctttAATGCTTGAAAAACATGCAGGACCCAAAAAAAAATAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                          <  2:786869/145‑1 (MQ=255)                                                                                                               cGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                       >  1:943533/1‑151 (MQ=255)                                                                                                                  ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                    >  1:931940/1‑151 (MQ=255)                                                                                                                        cctaaaaaaCATGAATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAatat                              <  2:796507/146‑1 (MQ=255)                                                                                                                          tGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAACTTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             >  1:344255/1‑150 (MQ=255)                                                                                                                          tGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGc                            >  2:45208/1‑151 (MQ=255)                                                                                                                             aaaaCATGCTAAAACCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa                         <  1:529856/151‑1 (MQ=255)                                                                                                                             aaaaCATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaa                          <  2:312159/150‑1 (MQ=255)                                                                                                                               aaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        >  1:678307/1‑150 (MQ=255)                                                                                                                               aaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                        >  2:775223/1‑150 (MQ=255)                                                                                                                               aaCATGCATGAGCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                       <  2:971180/151‑1 (MQ=255)                                                                                                                                aCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa                       >  2:349326/1‑150 (MQ=255)                                                                                                                                  atgcatgACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                     <  2:764614/150‑1 (MQ=255)                                                                                                                                     catgAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg                  <  2:71475/150‑1 (MQ=255)                                                                                                                                        gAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt               <  2:623904/150‑1 (MQ=255)                                                                                                                                         agcaaaaaaaaaaaaaaaaaAGTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagt               <  1:882289/143‑1 (MQ=255)                                                                                                                                           caaaaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa            <  1:600347/146‑1 (MQ=255)                                                                                                                                             acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         >  2:1122750/1‑151 (MQ=255)                                                                                                                                               aaaaaaaaaaaaaaaaTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       <  1:1039812/151‑1 (MQ=255)                                                                                                                                               aaaaaaATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATAtt        <  2:689788/150‑1 (MQ=255)                                                                                                                                               aaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       <  2:1038081/151‑1 (MQ=255)                                                                                                                                               aaaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtc       >  2:850506/1‑151 (MQ=255)                                                                                                                                                  ataataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    >  1:906858/1‑151 (MQ=255)                                                                                                                                                    aataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  2:101194/1‑151 (MQ=255)                                                                                                                                                    aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg   <  1:1098243/150‑1 (MQ=255)                                                                                                                                                     atatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  1:34765/1‑150 (MQ=255)                                                                                                                                                      tataAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  2:1098966/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                   TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTT‑AGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29347‑29642

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