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

TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29348‑29642                                                                                                                                                      |                                                                                                                                                   tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACGGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                   <  1:239847/150‑1 (MQ=255) tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                   <  2:537030/150‑1 (MQ=255)  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 <  1:342103/151‑1 (MQ=255)  aaTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaa                                                                                                                                                  >  1:549154/1‑150 (MQ=255)   aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                                 >  2:27146/1‑150 (MQ=255)     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               <  2:177224/150‑1 (MQ=255)     gTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaa                                                                                                                                               >  2:213145/1‑150 (MQ=255)      tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                              <  2:199135/150‑1 (MQ=255)       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTGTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                             <  2:530580/150‑1 (MQ=255)             tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                                       >  2:11868/1‑150 (MQ=255)              gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAt                                                                                                                                     >  2:882660/1‑151 (MQ=255)                ccATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCAtt                                                                                                                                    >  1:12282/1‑150 (MQ=255)                 cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                                   >  2:403087/1‑150 (MQ=255)                     tGGTCCTTTTTTTTCTGCTCTTTTTTATTTCATGAAATCTTTTTGCTGCGCAAGCTTTTTCC‑AGAACAGGTTAAATTATTTTTTTGTCGCTAAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGtt                                                                                                                               <  2:261084/150‑1 (MQ=255)                     tGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                              >  1:508520/1‑151 (MQ=255)                      ggTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                              >  1:576863/1‑150 (MQ=255)                             ttttttCTGCTCTTTTTTATTTCATGAAATTTTCTTGCTGCGCAAGCGTTTACC‑AAAACAGTTAAGATAATTTTTTTGTCGTTAACCCCCTGAAAAACATGCAGGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACttt                                                                                                                         <  2:948436/148‑1 (MQ=255)                             ttttttCTCCTCTTTTTTATTTCATGCATTTTTCTTGCTGCGAAAGCTTTTCCC‑AGAACAGGTTAGATGTTTTTTTTGTCGTTTAATGCCTGAAAAACAGGCATGAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACtttt                                                                                                                        <  1:680921/149‑1 (MQ=255)                               ttttCTGCTCTTTTTTTTTTCATAAAATTTTTTTGCTGCGCAAGCTTTTACA‑AGAAAAGGTTAGATAATTTTTTTGTCGTTTAATGCCTGTAAAAAATGCATGAGCCACAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAgc                                                                                                                     <  2:206970/150‑1 (MQ=255)                                 ttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGcc                                                                                                                  >  1:362061/1‑151 (MQ=255)                                  tCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                 >  1:387389/1‑151 (MQ=255)                                   cTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                 >  2:400500/1‑150 (MQ=255)                                    tGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                               >  2:250099/1‑151 (MQ=255)                                     gCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTATGTTGACTTTTAGCGCCCat                                                                                                               >  2:957144/1‑150 (MQ=255)                                       tCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtc                                                                                                            >  2:500683/1‑151 (MQ=255)                                          ttttttttTCATGCATTCTTCTTGCGCCAAAAGCTTTTTCA‑AAAACAGGTAAAATGTTTTTTTTGTCATTTAATCCCTAAAAAAAATGCATGTCCCAAAAAATAATTAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                                          <  2:540419/150‑1 (MQ=255)                                           ttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                         >  1:664772/1‑150 (MQ=255)                                             ttATTCAATACAATCTTCTTGCTACGCATGCGTTTCCC‑AAAATGGGTTAGATTATTTTTTATCCATTTAACGTTGGAAAAACACTGATTGCCCAAAAAAAAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                       <  2:118115/150‑1 (MQ=255)                                              tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                      >  2:812927/1‑150 (MQ=255)                                              tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                                      >  2:836090/1‑150 (MQ=255)                                                tttCATGCATTCTTTTTGCGGCGCAAGCTTTTTCC‑AGAACAGTTTAAATGTTTTTTTGGTCGCTAAATGCTTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                                    <  1:222096/150‑1 (MQ=255)                                                               ttGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                       >  2:132101/1‑148 (MQ=255)                                                                         atctttttCCAAGAACAGATTAAATTTTCTTTTGGTCTTTTAAGCCCTTAAAAAATTGCATGCCCAAAAAAAAAAAAAAAAAACTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaaa                                                                                  <  2:508515/140‑1 (MQ=255)                                                                           cGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                         >  2:562711/1‑150 (MQ=255)                                                                             ttttCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                                        >  2:322233/1‑149 (MQ=255)                                                                             ttttCC‑AAAACAGTTTAGATGATTTTTTTTTCGTTTAATGCCTATAAAACATGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                          <  2:386877/147‑1 (MQ=255)                                                                             ttttCC‑AAAACAGTTTAGATGATTTTTTTGTCGCTTAATGCCTAAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                          <  2:260682/147‑1 (MQ=255)                                                                                tacaaGAACAGGTAAAATGATTTTTTTTTCGTTTAATCCCTATAAACCATCCAAGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                     <  2:61013/146‑1 (MQ=255)                                                                                          ggTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                          >  2:13542/1‑151 (MQ=255)                                                                                            ttAGATGATTTTTTTGTCGTTTAATGCTTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTg                                                         <  1:851872/150‑1 (MQ=255)                                                                                              aGATTATTTTTTTGTCGTTAAATGCCTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                      <  2:531813/151‑1 (MQ=255)                                                                                               gatgatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     >  2:359055/1‑151 (MQ=255)                                                                                                   atttttttGTCGCTTAATGCCTGTAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     <  1:250101/144‑1 (MQ=255)                                                                                                    tttttttgcccttaaaTGCCTAAAAAACATGCATGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACg                                                  <  2:299454/135‑1 (MQ=255)                                                                                                    tttttttGTCTTTTATTGCTTAAAAAACGTGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                                     <  2:77559/144‑1 (MQ=255)                                                                                                    tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                >  1:441251/1‑151 (MQ=255)                                                                                                    tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                                >  2:840028/1‑151 (MQ=255)                                                                                                       ttttGTCGTTTAATGCCTGAAAAACATGAATGGCCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcat                                              <  1:484358/150‑1 (MQ=255)                                                                                                             ctcTTAATGCCTGTAAAACATGAATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                                          <  2:812123/146‑1 (MQ=255)                                                                                                               cTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                      >  1:113250/1‑150 (MQ=255)                                                                                                               cTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGa                                      >  1:170626/1‑150 (MQ=255)                                                                                                                ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                    >  2:651155/1‑151 (MQ=255)                                                                                                                taaatgctttaaaaaCATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGaa                                     <  1:928152/139‑1 (MQ=255)                                                                                                                       cTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                             >  2:292698/1‑151 (MQ=255)                                                                                                                           aaaaCATACATCACCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCa                           <  2:670934/149‑1 (MQ=255)                                                                                                                                    aggacccAAAAAATAATAAAAAAATCCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAag                    <  1:301782/142‑1 (MQ=255)                                                                                                                                       acccAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa                     <  2:549148/141‑1 (MQ=255)                                                                                                                                         ccACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           >  2:280658/1‑151 (MQ=255)                                                                                                                                         ccAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga             <  1:840041/149‑1 (MQ=255)                                                                                                                                           aaaaaataaaatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat         <  2:5396/149‑1 (MQ=255)                                                                                                                                           aaaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAata          <  2:277803/148‑1 (MQ=255)                                                                                                                                           aaaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAat           <  2:450143/147‑1 (MQ=255)                                                                                                                                           aaaaaaTAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgag                <  1:542251/142‑1 (MQ=255)                                                                                                                                             aaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAattt         <  1:27147/149‑3 (MQ=255)                                                                                                                                              aaaaaaTAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct      <  1:392154/151‑1 (MQ=255)                                                                                                                                               aaaaataAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc     <  2:631407/148‑1 (MQ=255)                                                                                                                                               aaaaaaataAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtctc     <  1:724536/145‑1 (MQ=255)                                                                                                                                                ataataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct    <  1:410694/151‑1 (MQ=255)                                                                                                                                                  aataAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  <  1:53725/151‑1 (MQ=255)                                                                                                                                                   atatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  2:542442/1‑150 (MQ=255)                                                                                                                                                   atatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTgg  >  1:825935/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCC‑AGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTGG  >  NC_000913/29348‑29642

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