Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	579,671	0	T	C	75.5%	95.3 / 28.5	54	intergenic (‑77/‑209)	ybcV/ybcW	DUF1398 domain‑containing protein YbcV/uncharacterized protein YbcW
Reads supporting (aligned to +/- strand):  ref base T (6/7);  new base C (22/18);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 7.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTG  >  NC_000913/579428‑579915                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     attaACACGGTGTTATCGGTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCGCAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGtt                                                                                                                                                                                                                                               <  1:238475/251‑1 (MQ=255)  ttaACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGttt                                                                                                                                                                                                                                              >  1:133422/1‑251 (MQ=255)         ggTGTTATCGTTTTCTAACACGATGTGAATAGTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTTGTAAATCCTTCCATGTGATACGATGGCGCGTAGTTTGCATTATCGTTTTGAGCGGTTGAATCTGGTCTGAGCTACTTGTGTTTTGTtgat                                                                                                                                                                                                                                       <  2:42506/251‑1 (MQ=255)                    tttCTAACACGATGTTAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTTACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGTTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTTCATTATCTTTTTTTTCGTTTCAATCTGGTCTGACCGCCTTGTGTGTTTTTGAGGATTTATGTc                                                                                                                                                                                                                             <  2:288361/250‑1 (MQ=255)                      tCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCaa                                                                                                                                                                                                                           <  1:19334/250‑1 (MQ=255)                      tCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCaa                                                                                                                                                                                                                           <  1:252108/250‑1 (MQ=255)                         aaCACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAatat                                                                                                                                                                                                                       >  2:147800/1‑251 (MQ=255)                         aaCACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAatct                                                                                                                                                                                                                       >  2:202391/1‑249 (MQ=255)                             cgaTGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAgg                                                                                                                                                                                                                   <  1:298569/249‑1 (MQ=255)                                  tgAATATTATCTATGGCTAGATAGAAAAAATAATGTGAGACGATGTGACGTTTTAGTTCAGAATAAAACAAATCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAACTGATCCATTGGAAAAACCTTCCATGTGATAAGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGTTTTATGTCAAATATTAGGAATGt                                                                                                                                                                                                              >  2:26664/1‑251 (MQ=255)                                   gAATATTATCTGTGGCTAGATAGTAAATATAATGGTAGACATTGTGACGTTTTAGTTTAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCTAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAACCCTTCCATGTGATACGAGTGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGTCCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGtt                                                                                                                                                                                                             <  1:378487/251‑1 (MQ=255)                                   gAATATTATCTGTGGCTAGATAGTAAATAGAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCTCACTTGATCGAATATTTCTTTAAAAATTGCATCCTGATCTATTGGTAAAGCCTTCCATGTGATACGAGGGTGCGTTGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCGTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGtt                                                                                                                                                                                                             <  1:428901/251‑1 (MQ=255)                                                      aTAGTACATATAAAGTGAGACGTTGTCACGTTTTCGTTCCGGATAAAACAATTACCAGTCTCAATCTTTTCGCACTTGATAGAATATATCTTTACAAATGGAAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGAATTGTCGTTTTTATCGTTTCAATCTGGTCTGACATCCTTGTGTTTTGTTGTTGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGt                                                                                                                                                                                          >  2:61153/1‑251 (MQ=255)                                                       tAGTAAATATAATGTGAGACATTGTGACGATTTAGTTCAGAATAACACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAACAATGGCAACCTGAGCCATTGGTAAAAACTTCCATGTGATACAAGGGAGCGTAGTTTGCATTATCATTTTTATCGCTTCACTCTGGTCTGACCTCTTTGTGTTTTGTTTATGATTTATGTAAAATATTAGGATAGTTTTCATTTATTAGTATTGGtt                                                                                                                                                                                         >  2:411786/1‑251 (MQ=255)                                                        aGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCGGAGCCATTGGTAAAACCTTCCATGTGATACGTGGTCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGTTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGtt                                                                                                                                                                                         <  1:69206/250‑1 (MQ=255)                                                           aaaTATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCg                                                                                                                                                                                      <  2:71358/250‑1 (MQ=255)                                                              tataATGTGAGACGTTGTGACGGTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTAGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTATTAGTATTGGTTGCGTAAc                                                                                                                                                                                  >  1:17301/1‑251 (MQ=255)                                                                 aaagtgAGACGTTGGGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGAAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACaa                                                                                                                                                                                >  2:76200/4‑250 (MQ=255)                                                                  aTGTGATACATGGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAg                                                                                                                                                                              <  2:109997/251‑1 (MQ=255)                                                                      gagacgctgtgACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCg                                                                                                                                                                          >  1:117584/8‑251 (MQ=255)                                                                            ttGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCctg                                                                                                                                                                    >  1:230155/1‑251 (MQ=255)                                                                              gtgACGTTTTAGTTCAGAATAAAACAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCctgct                                                                                                                                                                  >  1:274739/1‑251 (MQ=255)                                                                                  cGTTTTCGTCCAGAATAAAACAATTCACAGTCTAACTCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAAACTGAGCCATTGGTAAACCCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGc                                                                                                                                                               >  2:304009/1‑250 (MQ=255)                                                                                                  aaaaCCATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATAc                                                                                                                                              >  1:348444/1‑251 (MQ=255)                                                                                                  aaaaCAATTCACAGTCTAAATCTTTGCGCACTTGATCGAAGATTTCTTTAAATATGGCAACCTGAGCCATTGGTAAATCCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTTTCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATAGTAGGAAGGTTTTCACTTAATAGTATTGGTTGCGTATCAAAGTTCGGTCGTGCTGGCATTCTGGTGGGAAATa                                                                                                                                               <  2:260894/250‑1 (MQ=255)                                                                                                     aCAATTCACAGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAAc                                                                                                                                           <  1:297143/251‑1 (MQ=255)                                                                                                            acaGTCTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGAc                                                                                                                                       >  1:362960/1‑248 (MQ=255)                                                                                                                   aaaTCGTTTCGAACTTGATCGAATATTTCTTGAAAAATGGCAACCTGAGCCATTGGTATAATCTTCCATGTGATACGAGGGCGCGTAGTTTTCATTATCGTTTTTATCGTTTCAATCTGGTCTGATCTCCTTGTTTTTTGTTGTTGATTTATGTCAAATATTAGGAATTTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGTAAATTCAACCGCCAGatgta                                                                                                                                <  2:244985/248‑1 (MQ=255)                                                                                                                         tttCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAg                                                                                                                          <  2:18358/248‑1 (MQ=255)                                                                                                                         tttCGCACTTGAAAGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATAAGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGACTGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGcc                                                                                                                       >  2:386831/1‑251 (MQ=255)                                                                                                                          ttCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGcc                                                                                                                       <  1:269120/250‑1 (MQ=255)                                                                                                                              cACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGt                                                                                                                  <  1:158714/251‑1 (MQ=255)                                                                                                                                      cGAATATTTCTTTAAAAATGGCACACTGAGCCATTGGTAAAACCTTCTATGTGATAAGAGGGCGCGTAGTTTGAATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGCTGATTTATGTCACATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCaaat                                                                                                          >  2:454345/1‑250 (MQ=255)                                                                                                                                      cGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTTATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGCTGTATGTAAGGACAACGTGCTCAAAc                                                                                                          >  2:106928/1‑251 (MQ=255)                                                                                                                                      cGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTAGTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCaaa                                                                                                           <  2:126713/250‑1 (MQ=255)                                                                                                                                             ttCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCAt                                                                                                    >  1:6835/1‑250 (MQ=255)                                                                                                                                                tttAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGAAAGATGTATGAAAGGCCAACGTGCTCAAATCTTCATAc                                                                                                  >  1:1781/1‑249 (MQ=255)                                                                                                                                                  tAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGa                                                                                               <  1:144877/250‑1 (MQ=255)                                                                                                                                                    aaaaTGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGAGTAGTATACATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAAAATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGaa                                                                                              >  1:133487/1‑249 (MQ=255)                                                                                                                                                      aaTGGCAACCTGAGACATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGaa                                                                                              >  2:78235/1‑247 (MQ=255)                                                                                                                                                             aCCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGt                                                                                   <  1:195083/251‑1 (MQ=255)                                                                                                                                                               cTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTa                                                                                  >  1:246250/1‑250 (MQ=255)                                                                                                                                                                      aTTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTa                                                                          <  1:63408/251‑1 (MQ=255)                                                                                                                                                                      aTTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTa                                                                          <  2:362960/251‑1 (MQ=255)                                                                                                                                                                       ttGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTACTATTTTa                                                                          >  2:232167/1‑250 (MQ=255)                                                                                                                                                                                                   gcgcGTAGGTTGCATTATCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGATTGTTTTCACTTAATAGTATTGGTGGCGTAACAAAGGGCGGCCCTGCTGTCATTCTAGAGGGAAATACAGCCTACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATATAGAAATATTTGAATTGATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAg                                             <  1:215413/251‑1 (MQ=255)                                                                                                                                                                                                       gTAGTTTGCATTATCGTTTTTATCGTTTAACTCTAGTCTGACCTCCTTGTGTTTTGTTGTTAATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGTTTGCGTAACAAAGTGCGGTACTACTGGCATTCTGGAGGGAAATACAACCGACAGATGAAAGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGTTTTGAAGTAATATTTTAACCGCTAGATGCAGAGCAAGCGCATGGAGCgc                                          >  2:268610/1‑249 (MQ=255)                                                                                                                                                                                                             ttccttttCGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAAt                                    <  1:195444/243‑1 (MQ=255)                                                                                                                                                                                                             tGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACAGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGCTTTGAAGTAATATTTTAACCGATAGATGAAGAGCAAGCGCATGGAGCGACAAAATg                                   >  1:245650/1‑251 (MQ=255)                                                                                                                                                                                                                   aacgttgttaacgtttCAATCTGGTCGGACCTACTTGTGTTTTGTTGATGATTGATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGTAAGCGCATGGAGCGACAAAATGAATaa                              <  2:246250/236‑1 (MQ=255)                                                                                                                                                                                                                     cGTTTTTATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATaa                              >  1:108599/1‑248 (MQ=255)                                                                                                                                                                                                                           tATCGTTTCAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAAt                      <  1:67232/250‑1 (MQ=255)                                                                                                                                                                                                                                   cAATCTGGTCTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCtgatg             >  1:376311/1‑251 (MQ=255)                                                                                                                                                                                                                                              gACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTg  >  1:283614/1‑251 (MQ=255)                                                                                                                                                                                                                                                ccgctTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGaa                                                                                     <  2:229281/163‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     ATTAACACGGTGTTATCGTTTTCTAACACAATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACATTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTTTAAATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAACCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGCTTCAATCTGGTCTGACCTCTTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCAATTAATAGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGACAGATGTATGTAAGGCCAACGTGCTCAAACCTTCATACAGAAAGATTTGAAGTAATATTTTAACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGATGATCCCTCCGTG  >  NC_000913/579428‑579915

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