Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,072	0	A	C	55.4%	12.3 / 61.2	56	V8V (GTA→GTC)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (11/14);  new base C (19/12);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 2.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTT  >  NC_000913/1428830‑1429314                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     atgaCGGTCGCCACAGCAACCCGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGAAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGGTTTCAGGTGTACTGaa                                                                                                                                                                                                                                                 >  1:481989/1‑248 (MQ=255)          gCCACAGAAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAgg                                                                                                                                                                                                                                     <  2:137466/251‑1 (MQ=255)           ccACAGCAACACGTCTCTGGCGGGGGGAGCTGGCTTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGa                                                                                                                                                                                                                                    <  2:190533/251‑1 (MQ=255)                              gCGTGGGGAGCTGGAGTGACGTTTAACCCAAACGAATCCGTGGACATTGATATTGCTTATGCAGGCAACGGCAGGGGCGACTGGCGAACTGACGGTTTAATCGTGGGTGTCGGTTATAAGTTCTGATTAGTCAGGTACCACCGTGTTATGAAAGCCCGCCGGTAGAGGCGGGCTTTTTTGTGGGGTGAATCTGGCAGGACAGCTTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAAAACTGc                                                                                                                                                                                                                  >  2:244021/1‑250 (MQ=255)                                        cTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGTGTACTGAAAGACGGCAAAGGAAAAACGGGACAGAAATGCACAATCCAGCt                                                                                                                                                                                                       >  2:436614/1‑251 (MQ=255)                                                        tccggcagcgtccgtcgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAAc                                                                                                                                                                                       >  2:110908/21‑251 (MQ=255)                                                         ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGCCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACAAAGTGTTATGACAGCCCGACGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCCGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCCAAACGTAAc                                                                                                                                                                                       >  1:564329/1‑250 (MQ=255)                                                                 aTCCGGGGCCAGTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCGCAATCCAGCTGAAAGCAAAACGTAACAGCACCACg                                                                                                                                                                              <  2:465789/251‑1 (MQ=255)                                                                     gtcgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATTCAGCTGAAAGCCAGACGTCACAGCACCACggagg                                                                                                                                                                          >  2:447235/8‑248 (MQ=255)                                                                       cgttgttGATCTTGATTATGACGGCTCCGGCAGTGGCAACTGGCGTACTGTAGGATTACTCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCACGACGGACCCGGAGGGCTTTTTTGTGGGGGGAATATGGAAGTAAAGATTTCAGGCGACCTGACAGACGGCACAGGAAAACCGGCACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACCGCAACACggtg‑g‑tg                                                                                                                                                                      >  2:381028/6‑251 (MQ=255)                                                                       cgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACggtg‑g‑tg                                                                                                                                                                      >  1:187014/6‑251 (MQ=255)                                                                                       tATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTggg                                                                                                                                                        <  1:334552/249‑2 (MQ=255)                                                                                        aTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTc                                                                                                                                                     <  2:234528/251‑1 (MQ=255)                                                                                           aaGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCTCCACGGTG‑G‑TGGGGAACACGGTGGGCTCAGa                                                                                                                                                  <  1:160798/251‑1 (MQ=255)                                                                                                tccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTgggctacgagaatc                                                                                                                                             >  2:177645/2‑239 (MQ=255)                                                                                                  cGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCg                                                                                                                                           >  1:317236/1‑251 (MQ=255)                                                                                                   ggCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCgg                                                                                                                                          >  1:346426/1‑251 (MQ=255)                                                                                                      agttgCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAGAAATTCTGATTAGCCAGGTAACCCAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATAGGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATg                                                                                                                                       <  2:185104/247‑1 (MQ=255)                                                                                                       gTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCAACACGGTG‑GATGG‑GAACACGGTGGGCTCAGAGAATCCGGATgc                                                                                                                                      >  1:495757/1‑250 (MQ=255)                                                                                                                  cggaCTGACGGATGCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAATCCGGGCGtt                                                                                                                           <  2:460696/248‑1 (MQ=255)                                                                                                                   gtaCTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGtt                                                                                                                           >  1:175382/3‑250 (MQ=255)                                                                                                                        gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTGCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTTAAAGCAAAACGTAACAGCACCACTGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATTAATCCGGTCGTTACAGc                                                                                                                      <  2:378134/250‑1 (MQ=255)                                                                                                                           ggattCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATgg                                                                                                                  >  1:64105/4‑251 (MQ=255)                                                                                                                                   cGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGAGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGa                                                                                                           <  1:139373/250‑1 (MQ=255)                                                                                                                                            cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAg                                                                                                 <  2:171667/251‑1 (MQ=255)                                                                                                                                                tataAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAc                                                                                              <  1:110908/250‑1 (MQ=255)                                                                                                                                                      ttCTGATTAGCCAGGGAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGCCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGCACGGCCGGTACAgtgtcc                                                                                       >  2:315031/1‑246 (MQ=255)                                                                                                                                                          gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTACCAGCACCACGGAG‑G‑TGGTGAACACGGTGGGCTAAGAGAATCCGGATGAAGACGGGCGTTACAGCATGGATGGGGAGGACGGTCAGTACAgtggcaccct                                                                                   >  2:349461/1‑242 (MQ=255)                                                                                                                                                            ttAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGt                                                                                 >  2:125431/1‑251 (MQ=255)                                                                                                                                                                ccgggTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCGTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAGCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTTAGTACGGTCAGTACAGCGCTATTCTGTtggt                                                                             <  2:481989/248‑1 (MQ=255)                                                                                                                                                                 cAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggtt                                                                            >  2:163313/1‑235 (MQ=255)                                                                                                                                                                  aGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttg                                                                           <  1:163313/251‑18 (MQ=255)                                                                                                                                                                      aaCACAGTGTGATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGGACAgtgtcatcctgcaggttgacg                                                                        <  2:62126/250‑21 (MQ=255)                                                                                                                                                                         acaGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggttt                                                                    >  1:448740/1‑227 (MQ=255)                                                                                                                                                                                           gccgGTTCAGGCGGGTTTTTTTGGGGGGGGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCAAAGGAAAACCGGGACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGGGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGATTACGGTCATTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCgg                                                  <  2:447809/251‑1 (MQ=255)                                                                                                                                                                                               gtccAGGCGGGCTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCGGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAcc                                              <  1:537475/248‑1 (MQ=255)                                                                                                                                                                                                 gcgggtgggCTTTTGTGTGGGGGGAATAGGTCAGGAAAGATTTCAGTTGTACTTAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCGGCTGAAAGAAAAACGAAACAGCCCCA‑GGTGTG‑GGGTGAACGCGCTGGCTTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAGGGTCATTTCAGCGGGATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAt                                            <  1:34094/245‑1 (MQ=255)                                                                                                                                                                                                    ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                         >  1:320049/1‑251 (MQ=255)                                                                                                                                                                                                       gggCTTTTTTGTGGGGTGAAGATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGAAAAACGTAACAACACCACGGTG‑G‑TGGTTAACACGCTGGCCTCGGAAAATCCGGATGAACCCGGGCGTGACAGCATGGACGTTGAGTAGGGTAAGTACAGCGTTATTCTGTTGGTGGAAGTATTCCCGCCGTCACATGCCGGGACCATTACCgt                                      <  1:369497/251‑1 (MQ=255)                                                                                                                                                                                                        ggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGCCCATTACCgtg                                     >  2:744/1‑251 (MQ=255)                                                                                                                                                                                                         gCTTTTTTGTGGGGTGAATATGGCAGTAAAGCTTTCAGGTGTACTGAAAGACGGCCCAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑AGGCGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTACCATGCCGGGACCATTCCCgtgt                                    >  2:284960/1‑251 (MQ=255)                                                                                                                                                                                                              tttGTGGGGTGAATATGGCAGTAAATATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCTCCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGTGCTCATAGAATCCGGATGAAGCCGGGCGGTACAGCATGGATGTGTTGTACGGTCAGTACAGTGTCATCCTGCGGGGTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGa                                <  1:43405/250‑1 (MQ=255)                                                                                                                                                                                                               ttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGaa                               <  1:366256/250‑1 (MQ=255)                                                                                                                                                                                                                  tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTCCTGAAAGCCGGAACAGGAACACAGGCACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCAAGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACGGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCAAATGCCGGGACCATTACCGTGTATGAAGAtt                           >  2:561195/1‑251 (MQ=255)                                                                                                                                                                                                                  gggggTGAATAGGGCAGTAAAGATTTCAGGTGGACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGATGAAAGCAAAACGTAACAGCACCGCGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAAACCGGTTGAAGCCGGGCGTTACAGCATGGAGGTGGAGTACGTTCAGTACAGCGTTATTTTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGGATGAAGAtg                           <  2:182601/250‑2 (MQ=255)                                                                                                                                                                                                                    gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACCGGAAAACCGGTAAAGAACTGCACCATTCAGCTGCAAGCCAGACGTAACAGCCCCACGGTG‑A‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGCGTACGGCAAGTACAGTGTCATCCTGCCGGCTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTtc                          >  2:493375/1‑250 (MQ=255)                                                                                                                                                                                                                    gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAACATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGCTGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCAACGTGTATGAAGATtcc                         >  2:502073/1‑250 (MQ=255)                                                                                                                                                                                                                      gTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAACACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGCAGATtctc                        >  1:26086/1‑250 (MQ=255)                                                                                                                                                                                                                       tGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGACAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcaca                       >  1:502116/1‑247 (MQ=255)                                                                                                                                                                                                                           aatGGCTGTAAAGATTGCAGGTTAACTGAAAGCCGGCACAGGAAAACCGGTACAGAACTGCAATACCCAGCGGAAAGCAAAACTGAAAAGAACCACGGTG‑G‑GGGTGAACACGCGGGCCTCAGAAAATCCGGAGGAAGCCGGGCGTTACAGCATTGACGTGGAGTCGGGTCAGTCCAGCGTTATTCTTTTGGTGAAATTATTCCCGCCGGCACATGCCGGGACCATTACCGTGTATGAAGATGCTCAACCCg                  <  2:510030/250‑1 (MQ=255)                                                                                                                                                                                                                                cAGTAAAGATTTCAGGTGTCCTGAAAGACGGCACAGGAAATCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTG‑G‑TGGTGACCCCGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg             <  1:137197/251‑1 (MQ=255)                                                                                                                                                                                                                                cAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggggacg             <  2:198451/251‑8 (MQ=255)                                                                                                                                                                                                                                   tAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTg          >  2:156236/1‑251 (MQ=255)                                                                                                                                                                                                                                    aaaGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCAGCGCACGCCGGGACCATCACCGGGTATGAAGATTCACAAccgtggacgctga         <  2:88738/251‑12 (MQ=255)                                                                                                                                                                                                                                       gATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATg      >  1:180217/1‑251 (MQ=255)                                                                                                                                                                                                                                           tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt  >  1:517741/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     ATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTT  >  NC_000913/1428830‑1429314

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