Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,171,781	0	A	G	20.6%	54.3 / 12.6	34	intergenic (+25/‑276)	rrfB/murB	5S ribosomal RNA/UDP‑N‑acetylenolpyruvoylglucosamine reductase
Reads supporting (aligned to +/- strand):  ref base A (19/8);  new base G (3/4);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.18e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAA  >  NC_000913/4171533‑4172025                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      aGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaa                                                                                                                                                                                                                                                      >  1:500849/1‑250 (MQ=37)               aaCGCCGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCg                                                                                                                                                                                                                                       >  1:449702/1‑251 (MQ=9)                    cGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGCACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGTTCTCCCCATGCGAGGGTAGGGACCTGCCAGGCATCAAATAAAACGCAAGGCTCAGTCGAAAGACAGGGCCTTTCGttttt                                                                                                                                                                                                                                  >  2:116029/1‑250 (MQ=255)                    cGAAGATGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTa                                                                                                                                                                                                                                  <  1:323216/251‑1 (MQ=9)                      aaGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑c                                                                                                                                                                                                                               >  1:586007/1‑251 (MQ=37)                       aGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCCGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑Ct                                                                                                                                                                                                                              >  1:396505/1‑251 (MQ=37)                         cTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑Ctgt                                                                                                                                                                                                                            >  1:97106/1‑251 (MQ=35)                          tGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGAGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGCACTCAGACGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAATACGAAAGGCTCCGTCGGAAAACTGGGCCTTTCGTTTTAT‑Ctttt                                                                                                                                                                                                                           >  2:156970/1‑248 (MQ=11)                            gtttGGCGGATGAGAGAAGATGTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGTCGGCCGGAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTGCGTTTTAT‑Ctgttgt                                                                                                                                                                                                                         <  2:59195/250‑1 (MQ=11)                                gtcGGATGAGAGAAGATTTTCAGCATGATACAGATAAAATCAGAAGGCAGAAGCGTGCTGATGAAACAGAATTTGCCGGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGTGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCTTCAAATAAAACTAAAGGCTCAGTCGGAAGACTTGGCCTTTCGGTTTAT‑CTGttgtttgt                                                                                                                                                                                                                     <  2:418181/249‑1 (MQ=11)                                  cGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGTTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAACCTAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCg                                                                                                                                                                                                                   >  1:689546/1‑251 (MQ=21)                                    gATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATGAAACAGAATTTGCCTGGCGGCCGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGt                                                                                                                                                                                                                 <  2:544788/251‑1 (MQ=11)                                       gagagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGa                                                                                                                                                                                                               >  1:209217/1‑250 (MQ=33)                                              gaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCtctc                                                                                                                                                                                                       <  1:508420/251‑1 (MQ=33)                                                   tCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTg                                                                                                                                                                                                    >  1:603213/1‑249 (MQ=33)                                                        cTGATACAGATTAAATCAGAACGCAGCAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGGGTAgg                                                                                                                                                                                              >  2:679265/1‑250 (MQ=18)                                                                   tAAATCAGAACGCAGACGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGAAAAATccgccg                                                                                                                                                                                  >  2:97198/1‑251 (MQ=18)                                                                            aCGCAGAAGCGGTCTTATAAAACAGAATTTGCCTGGCGGCAGTATCGCGGTGGTCCCACCTGACCCCATGCCGAACTCATAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATTCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAt                                                                                                                                                                         <  2:439975/251‑1 (MQ=16)                                                                            aCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAt                                                                                                                                                                         >  2:326908/1‑251 (MQ=33)                                                                                      ggTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑Ctgttg                                                                                                                                                                                                                          <  1:396747/192‑1 (MQ=34)                                                                                      ggTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑Ctgttg                                                                                                                                                                                                                          >  2:396747/1‑192 (MQ=21)                                                                                           gATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGGGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg                                                                                                                                                          <  2:9455/251‑1 (MQ=18)                                                                                               aaaCAGAATTTGCCTGGCGGCAGTAGCGCGGAGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGACAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAAc                                                                                                                                                      >  1:491721/1‑251 (MQ=18)                                                                                                  cAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGc                                                                                                                                                   <  1:534251/251‑1 (MQ=33)                                                                                                        ttGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGCGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                                                                                                                               >  2:281905/1‑249 (MQ=18)                                                                                                                   cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGCACTCAGAAGTGAAACGCAGTAGCGCCGATGGTCGTGTGGGATCTCCACATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑ATGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGTGGGTGGGGGGcg                                                                                                                                  >  2:406893/1‑250 (MQ=11)                                                                                                                                                          aaaaGTGAACCCCCGAAGCGCGGAGGGTGGTGTGGGGTCCCCCCGTGCGAGAGTAGGGAACTGCCAGGCATCAAAAAAAACGAAAGCCTAAGTCGAAAGTCTGGGCCTTTGTTTTTAT‑TTGTTGTTTGTCGGTGAACGCTCTCCTGAGAAGGACAAACCCGCCGGGTGCGTATTTAAACGTTGAAAAGCAACGGCCGGGAGGGTGGCGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAg                                                                                           <  1:511521/249‑1 (MQ=2)                                                                                                                                                                                    gTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCCTCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCtttt                                                                 >  2:71082/1‑251 (MQ=18)                                                                                                                                                                                         gtgGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAt                                                                                                                                                                         <  1:77113/142‑1 (MQ=35)                                                                                                                                                                                         gtgGGGTCTCCCCATGCGAGAGCAGGGAACTGCCAGGCCTCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGTGCGGAt                                                                                                                                                                         >  2:77113/1‑142 (MQ=255)                                                                                                                                                                                                          gagaGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGCGAAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCCTCACATTAAACAGAAGGCCATCCTGACGGGTGGCCTTTTTGCGTTTCTACAAACTCTTccc                                           >  1:71019/1‑250 (MQ=14)                                                                                                                                                                                                                       tGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGCGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAACGCCATCCTACCGGATGGCATTTTTGCGTTTTTACTAAATCTTCCCCTCGTCatttatc                              >  2:162546/1‑246 (MQ=11)                                                                                                                                                                                                                       ggccaggaaaaaaaTAAAACGAAAGCCTCAGTCAAAAGCCTGGGTCTTT‑GTTTTTTGCTTTTGTTGGGGGGTGACCGCTCCCTTAAGGAGGACAACCCCGCGGGGAGCGGATTTGAACGTTGCGAAGCAAGGGCCGGGGGGGGGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTa                              <  1:326908/240‑1 (MQ=9)                                                                                                                                                                                                                                           gAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGACCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGCTGCGAAGCAACGGCCCGGAGGGGGGCGGGCCGGACGCCCGCCATAAACTGCCAGGCATCAACTTAAGCAGAAGGCCAGCCTGACGGATGGCCTTTTTGCGTTTCTAACAACTCTACCTGGCGTCATATCTACAAGCCATCACCCCACCGAt          >  2:163997/1‑251 (MQ=14)                                                                                                                                                                                                                                                   cAGTCGGAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGATCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCCGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGACATATCTACAAGCCATCCACCCACAGATACGGTaac  >  2:354297/1‑250 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      AGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTAT‑CTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCGTTTCTACAAACTCTTCCTGTCGTCATATCTACAAGCCATCCCCCCACAGATACGGTAAA  >  NC_000913/4171533‑4172025

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