Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,160	0	C	G	67.2%	104.6 / 51.5	64	L38V (CTG→GTG)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (14/7);  new base G (24/19);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGT  >  NC_000913/1428914‑1429400                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 cTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGGGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGAGGTGGTGAACACGGTggg                                                                                                                                                                                                                                              >  1:51030/1‑250 (MQ=255)  ttATGAGGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGCAATGCCCAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACAAGCTGGcc                                                                                                                                                                                                                                             >  2:167829/1‑251 (MQ=255)   tATGAAGGCTCCGGCAGTGGCTACTGGCGTACTGACGGATTCAGCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTggg                                                                                                                                                                                                                                              <  2:369488/249‑2 (MQ=255)    agtatgtcccgggCAGTGGCGACTGGCGCATTGACGGTTTCTGCGTGGGTTTCGGTTATAAGTTCTGATTTGCCTGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCt                                                                                                                                                                                                                                            <  2:354781/239‑1 (MQ=255)       aaGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTc                                                                                                                                                                                                                                           >  2:124966/1‑248 (MQ=255)           ctccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGAACCACGGTGGTGGTGAACACGGTGGGCTCAGagact                                                                                                                                                                                                                                    >  2:126145/3‑247 (MQ=255)             ccGGCAGTGGCGACTGGCGGACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTTAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATcc                                                                                                                                                                                                                                  <  1:292426/251‑1 (MQ=255)                 cAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGAGTAGCCAGGTAACACAGTTTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCGCCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGAt                                                                                                                                                                                                                              <  1:286365/251‑1 (MQ=255)                    tGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGa                                                                                                                                                                                                                            >  1:114762/1‑250 (MQ=255)                           tGGCGCACTGACGGTTTAATCGTGGGTGCCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGCCAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGGGTACTGAAAGACGGCCCAGGAAAACCGGTACCGAACTGCACAATCCAGCTGAAAGCAAAACGTAAAAACACCACGGTGGGGGTGAACACGCAGGCCTCAGAAAATCCGGATGAAGCAgggg                                                                                                                                                                                                                    >  2:366390/1‑250 (MQ=255)                               gtaCTGACGGCTTCATCGTTGGCGTCGGATATAACTTCTGATTAGCCAGGTAACACAGCGTTATGACAACACGCCGGAACCGGTAGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCATGAAAGACGGAACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACAAGGTGGGCTCAGAGCATCAGGATCAAGCCGGGCGtt                                                                                                                                                                                                                 >  2:48353/3‑250 (MQ=255)                                  ttGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGGAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGACAGACGGCCCAGGAAAACCGGTACAGAACTGCCCCATTCAGCTGAAAGCAAGACGTAACAGCACCACGGTGGTGGTGAGCACCGGGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACa                                                                                                                                                                                                              >  2:80997/2‑250 (MQ=255)                                   tGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGc                                                                                                                                                                                                            >  1:227407/1‑251 (MQ=255)                                                gttggggtCGGTTATAAATTATGATTCGCCAGGCAACACAGTGTTATGACAGCCCGCCGGACCCGGTGGGCTTTTTTGTGGGGTGCATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCACAGAGAATCCGGATGAAGCCGGGCGTTACAACATggctgtggagg                                                                                                                                                                                               >  2:298073/7‑242 (MQ=255)                                                   ggggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTAc                                                                                                                                                                                             <  1:426892/247‑1 (MQ=255)                                                    tggtctgTTATAAATTCTGATTAGCCAGGTACCATAGTGTTATGACAGCCCGTCGGACCCGGTGGGCTTTGTTGTGGGGTGAATATGGTAGTAAAGATTTCAGGCGTCCTGAAAGATGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACGCGGTGGGCTCAGAGAAGCCGGTTGAAGCTGGGCGTTACAGCATGGATGTGGAGTACg                                                                                                                                                                                            <  2:223864/244‑1 (MQ=255)                                                     ggtCGGTTATAAGTTCTGAGTAGCCAGGTAATACAGGGTGATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTTAAAGCAAAACGTAACAGGACCACGGTGGTGGTGAGCACGCTGGCGTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGt                                                                                                                                                                                          <  2:191792/250‑1 (MQ=255)                                                     ggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGt                                                                                                                                                                                          >  1:72875/2‑251 (MQ=255)                                                               aaattCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtg                                                                                                                                                                                <  1:350811/248‑3 (MQ=255)                                                                 attCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGTTTTTTTGTGGGGTGAATATGGCAGTATAGATTTCAGGAGGCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGTATGTGGAGTACGGTCAGTACAgtgtc                                                                                                                                                                              <  2:192808/250‑5 (MQ=255)                                                                      gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcc                                                                                                                                                                          <  1:14988/250‑9 (MQ=255)                                                                      gATTAGCCAGGTAACACAGTGTTATGACAGCCCACCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAAAGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCt                                                                                                                                                                         >  1:344460/1‑251 (MQ=255)                                                                      gATTAGCAAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGAGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcct                                                                                                                                                                         >  1:324221/1‑242 (MQ=255)                                                                                gTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgac                                                                                                                                                               <  2:9245/251‑20 (MQ=255)                                                                                   acacAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCACGACGTAACGGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTCCGGTCCGTACAgtgtcatcctgcaggatgacggt                                                                                                                                                            >  1:356569/1‑229 (MQ=255)                                                                                   acacAGTGTTATGACAGCCAGCCGGAACCGGTGGGCTTTTTTGCGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGCGAATCCGGATGCAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAccgtgtcatcctgcaggttgacggt                                                                                                                                                            >  2:54950/1‑227 (MQ=255)                                                                                    cacaGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggtt                                                                                                                                                           >  1:391761/1‑228 (MQ=255)                                                                                      caGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggtttt                                                                                                                                                         <  2:314857/251‑26 (MQ=255)                                                                                                                  cggacTTTTTTGTGGGGTGAATATGGCAATAAAGATTTCAGATGTACTGAACGACGGGACAGGAAAACCGGTACAGAACTGCACAATACAGCTGAACGCAAAAAGTAACAGCACCACGGTGGGGGTGAACACGCTGGCCTCAGAAACTCCGGCTGAAGCAGGGCGTTACAGCATGGACGTTGAGTACGGACAGTACAGCGTTATTCGGTTGGTGTAAGGATTCCCGCCGTCACATGCCGGGACCATTCCCg                                                                                                                             >  2:456643/5‑251 (MQ=255)                                                                                                                   gggtttttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGAGGGCATAGGAAAACCGTGACAGAACTGAACAACCCAGCTGAAAGCAAAACTTAGCAGCTCCGGGGTGGTGGTGAGCACGCTGGCCTCAGAAAATCCGGATGAAGCGGGGTGTTACATCTTGGACTTTGGGTACGGTAAGGGTAGGGTTATTCTGTTGTGGGAAGTGTTTCCGCCGTCACATGCCGGGAGCATTACCgt                                                                                                                            <  2:251722/247‑1 (MQ=255)                                                                                                                       ttttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTAt                                                                                                                        >  1:12564/1‑251 (MQ=255)                                                                                                                        ttgtgggggggTGAATAGGGCAGTAAAAATTTCAGGTGTACTAAAAGACGGCCCAGGAAAACCGGTACAAAACTGCACAATCCAGCTAAAAGCAAAACGTAACGGCACCACGGGGGTGGTGAACACGCGGGCCTCAGAAAATCCGGAGGAAGCCGGGTGTTACAGCATGGACGTTGAGTCCGGTCAGTACAGCGTTATTCTTTTGGTGTAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTa                                                                                                                         <  2:377348/242‑1 (MQ=255)                                                                                                                               ggggAGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGAAAAACCGGTACAGACCTGCACACTCCAACTGAAAGCAAAACGTAACAGCACCACGGTGGAGGTGAACACGCTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGCCGTTGCGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATtc                                                                                                                >  2:297954/1‑251 (MQ=255)                                                                                                                                      atatGGCAGTAATGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGAGTATGAAGATTCACAAccg                                                                                                         >  1:228729/1‑250 (MQ=255)                                                                                                                                      atatGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAAccc                                                                                                         >  1:200756/1‑251 (MQ=255)                                                                                                                                      atatGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCAGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccg                                                                                                         >  1:216607/1‑250 (MQ=255)                                                                                                                                        atGGCAGTAAAGATTTCAGGAGTTCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccgg                                                                                                        <  2:336621/250‑3 (MQ=255)                                                                                                                                          ggCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccgggg                                                                                                      <  2:107834/250‑5 (MQ=255)                                                                                                                                           gCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggggac                                                                                                    <  2:271516/251‑7 (MQ=255)                                                                                                                                              gTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCt                                                                                                 >  1:50414/1‑251 (MQ=255)                                                                                                                                                  aGATTTCAGGAGTTCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCTGGATTAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGGACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCTCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAAt                                                                                             <  1:444898/251‑1 (MQ=255)                                                                                                                                                  aGATTTCAGGAGTCCTGAAAGACGGCACAGTAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGGGGAGAACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAAt                                                                                             >  2:107083/1‑251 (MQ=255)                                                                                                                                                     tttCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGACCTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGAGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGAt                                                                                          >  1:194147/1‑251 (MQ=255)                                                                                                                                                       tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                                                                        >  1:204163/1‑251 (MQ=255)                                                                                                                                                               aCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTg                                                                                >  1:251064/1‑251 (MQ=255)                                                                                                                                                                  gAAAGAAGGCACAGGAAAACCGGTACAGAACAGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCGGGCCACAGAAAATCCGGATGAAGCCGGGCGTTACAGACTGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATACCCGCCGTCACATGCCGGGACCATTCCCGTGTATGAAGATTCTCAACCCGGGACACTGAATGATTTTCTCggagca                                                                              >  2:356267/1‑246 (MQ=255)                                                                                                                                                                       aCGGCACACGAAAACCGGTACAGAACTGCACAATCCAGCTGCAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACCCGCTGGCCTCAGAAAATCCGGATAAAGCCGGGCGTTACAGCCTGGACGTTGAGTACGGTCCGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGCTTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCCTGACg                                                                        >  2:380228/1‑251 (MQ=255)                                                                                                                                                                        cGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                       >  1:102459/1‑251 (MQ=255)                                                                                                                                                                         ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgga                                                                      <  2:72875/251‑1 (MQ=255)                                                                                                                                                                         ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                       >  1:73145/1‑250 (MQ=255)                                                                                                                                                                         ggCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACgg                                                                       >  1:399198/1‑250 (MQ=255)                                                                                                                                                                                    aCCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGccc                                                           <  1:90255/251‑1 (MQ=255)                                                                                                                                                                                        gTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggcc                                                       <  1:348194/251‑4 (MQ=255)                                                                                                                                                                                             gAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACAGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCTGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTGCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAgg                                                  <  1:54486/251‑1 (MQ=255)                                                                                                                                                                                             gAACTGCACAATCCAGCTGAAAGCAAAACGTCACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAgg                                                  >  1:223011/1‑251 (MQ=255)                                                                                                                                                                                                       acccAGTTGAAACCAAAACGTACAAGCCCCAGGGTGGTGGTGAACACGCGGGCCTCAGAAAATCCGGATGAGCCCGGGCGTTACAGCATGGACGTTGAGTTCGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCg                                        <  2:317456/249‑1 (MQ=255)                                                                                                                                                                                                         ccAGCTGAAAGCAAAACGTAACAGAACCAGGGGGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGTGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtt                                      <  1:16027/251‑1 (MQ=255)                                                                                                                                                                                                                aaagccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCACGGCCGGAGGTGCTGCGTCGTCTTGAACt                               >  1:86258/9‑251 (MQ=255)                                                                                                                                                                                                                 aagccagaCGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTg                              >  1:24147/8‑251 (MQ=255)                                                                                                                                                                                                                            aaCAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAag                    <  1:141190/250‑1 (MQ=255)                                                                                                                                                                                                                              cAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACTGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAgg                 >  1:48075/1‑251 (MQ=255)                                                                                                                                                                                                                                  accacGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGgc             >  1:63478/1‑251 (MQ=255)                                                                                                                                                                                                                                           gtggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAAcgc    <  2:12564/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggtgAACACGCTGGACTCAGCAAATCCGGATGAAGCCGGGAGATACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCCTGACGGAGGATGATGCCCGTCCGGAGGCACTGAGCCGCTTTGAACTGAGGGGGGAAGAGGTGGCGCGTAACGCGt  >  2:274331/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 CTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGT  >  NC_000913/1428914‑1429400

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: ‑