Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,164	0	C	G	56.7%	34.0 / 78.7	60	A39G (GCC→GGC)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (18/8);  new base G (18/16);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCAC  >  NC_000913/1428923‑1429412                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         gctccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGagaa                                                                                                                                                                                                                                                 <  1:236472/248‑4 (MQ=255)    ccGGCAGTGGCGGCTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATc                                                                                                                                                                                                                                               <  2:94541/250‑1 (MQ=255)    ccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATc                                                                                                                                                                                                                                               <  2:219515/250‑1 (MQ=255)              cGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCAAGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGcc                                                                                                                                                                                                                                    >  1:62573/1‑251 (MQ=255)                   ggCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCAACACGGTGGTGGTGAACAAGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCg                                                                                                                                                                                                                               >  2:4148/1‑251 (MQ=255)                          tGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGAACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTCCAGc                                                                                                                                                                                                                        >  1:221935/1‑251 (MQ=255)                           gACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAg                                                                                                                                                                                                                         <  1:151649/249‑1 (MQ=255)                             cggattCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATg                                                                                                                                                                                                                     <  1:265776/247‑1 (MQ=255)                             cGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATg                                                                                                                                                                                                                     >  1:19160/1‑251 (MQ=255)                               gTTTCATCGTCGGTATCGGTTATCAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGATTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGGGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GAAAAACGTAACAGGACCACCTTGGTGGAGAACACGCTGGCCTAAGAAAATCCGGATACAGCCGGTCGTGACAGCATgtg                                                                                                                                                                                                                   >  2:27811/1‑249 (MQ=255)                                actCATCGTTGGGGCCGGTTATAAATTCTGATTAGACAGGTACCACAGTGTTATGACAGCCCGCCGGAACCGGTGGGTTTTTTCGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTGCTGAAAGACGGCACAGGGAAACCGGGAAAGAAATGCACCATTAAGCTGAAA‑GCCAGACGGAAAAGCACCACAGTGGAGGTGAAGAAGGTGGGCTAAGAGAATCAGAATGAAGACGGGCGATAAAGCaanga                                                                                                                                                                                                                   >  2:252723/3‑246 (MQ=255)                                      cgttggggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGatgtggag                                                                                                                                                                                                            <  2:230321/244‑8 (MQ=255)                                         tggggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTa                                                                                                                                                                                                          <  1:64067/246‑1 (MQ=255)                                           gggtCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACgg                                                                                                                                                                                                       <  2:95662/249‑1 (MQ=255)                                                   tataAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAc                                                                                                                                                                                                >  1:178478/1‑250 (MQ=255)                                                    ataAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCTGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGGAACCGCACCACGGTGGTGGTGAACACGCTGACCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTAc                                                                                                                                                                                                >  1:110224/1‑249 (MQ=255)                                                     taAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGc                                                                                                                                                                                             >  2:30462/1‑251 (MQ=255)                                                           cTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTAtt                                                                                                                                                                                       >  1:61033/1‑251 (MQ=255)                                                                tAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTTCAGCGTTATTCTGt                                                                                                                                                                                   <  1:192861/250‑1 (MQ=255)                                                                 aGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGGGGAGTACGGTCAGTACAgtgtcatcctnc                                                                                                                                                                                   >  2:259729/1‑238 (MQ=255)                                                                 aGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGGGGAGTACGGTCAGTACAgtgtcatcctgc                                                                                                                                                                                   >  2:141663/1‑238 (MQ=255)                                                                    cAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggtt                                                                                                                                                                              >  1:45618/1‑235 (MQ=255)                                                                       gTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGt                                                                                                                                                                                                  <  2:83077/228‑1 (MQ=255)                                                                       gTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGt                                                                                                                                                                                                  >  1:83077/1‑228 (MQ=255)                                                                              aGTGTTATGACAGCACGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAAAGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTcc                                                                                                                                                                    >  2:36495/1‑251 (MQ=255)                                                                                               ccgGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GAAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTGGGTGGAAGGATTCCCGCCGTCACATGCCggg                                                                                                                                                   <  1:229607/251‑1 (MQ=255)                                                                                                           ggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCgtg                                                                                                                                       >  2:206035/1‑251 (MQ=255)                                                                                                                         gTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcaca                                                                                                                         <  2:62573/251‑4 (MQ=255)                                                                                                                          tGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACCGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcaca                                                                                                                         >  2:225220/1‑247 (MQ=255)                                                                                                                           gAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGGAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAAAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAAc                                                                                                                       >  2:99797/1‑251 (MQ=255)                                                                                                                                 ggaaGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGTACTGCACCATTCAGCTGAAATGCAGA‑CGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTAAGAGAATCCGGATGAAGCCGGGCGTGACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGTCGGGGCCATCACCGTGTATGAAGATTCCCAAccgggga                                                                                                                 <  2:227830/248‑6 (MQ=255)                                                                                                                                   aaGTAAAGATTTCAGGTGTAAGGAAAGACGGCCAAGGAAAACGGTTACAGAACGGCACAATCCAGCTGAAA‑GCAAACGGTAAGAGCACCACGTTGGTGGTGAACGCGCTGGCCTCGGAAAATCCGGATGAAGCCGGGCGTTACAGCTTGGACGTTGAGTACGGTCAGTCCAGCGTTAGTCTGTTGGTGGAAGTGTTGCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg                                                                                                               <  1:109299/250‑1 (MQ=255)                                                                                                                                     gTAAAGATTTCAGGTGTACTGAAAGACGGCACAGAAAACCCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAACCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCt                                                                                                             <  1:55765/251‑1 (MQ=255)                                                                                                                                       aaaGATTTCACGTGTACTGCAAGACGGCACAGGAAAACCGGAACAGCTCTCCAAACTCAAGCTGAAA‑GCAAAACGTAACAGCACCACTGAGGTGGTGCACCCACTGGCCTCAGAACATCCTGATGAAGCCGGGCGCTCCAGCCTGGACGTTGCGTACGGTCCGTACCGCGTTCTTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGATGATTCTCAACCCGGTACGCTgc                                                                                                           >  2:42664/1‑250 (MQ=255)                                                                                                                                         aGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAAt                                                                                                         <  1:141663/251‑1 (MQ=255)                                                                                                                                              tcagtagtccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACCGCTCCACGGTGGTGGTGAACACGGTGGGCTCAGAGACTCCGGATGAAGCCGGGCTTTACa                                                                                                                                                                                                                          <  2:129091/124‑1 (MQ=255)                                                                                                                                              tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACa                                                                                                                                                                                                                          >  1:129091/1‑133 (MQ=255)                                                                                                                                                   tgtACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCg                                                                                               >  1:41676/1‑251 (MQ=255)                                                                                                                                                         gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGcc                                                                                          >  2:214565/1‑250 (MQ=255)                                                                                                                                                           aaGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATg                                                                                       >  2:136597/1‑251 (MQ=255)                                                                                                                                                             gACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGAc                                                                                     >  2:77606/1‑251 (MQ=255)                                                                                                                                                                ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                                                                                   >  1:124276/1‑250 (MQ=255)                                                                                                                                                                      ggAAAACCGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGgatga                                                                            <  1:90499/251‑1 (MQ=255)                                                                                                                                                                           aCCGGTACAAAACTTCATAATCCAGCTGAAC‑GCAAAACGTAACAGCACCACGGTGGTGTGGAACCCCCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTTTGGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGGTTATGAAGATTCTCAACCCGGTACGCTGAATGATTTGCTCGGTGCCATGACGGAGGATGATGccc                                                                       <  1:248010/251‑1 (MQ=255)                                                                                                                                                                             cGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCGCCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCgg                                                                     <  2:49299/251‑2 (MQ=255)                                                                                                                                                                             cGGTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCgg                                                                     >  1:92518/1‑250 (MQ=255)                                                                                                                                                                              ggTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggc                                                                    >  2:27293/1‑249 (MQ=255)                                                                                                                                                                               gTACAGAACTGCACCATTCAGCTGAAA‑GCCAGACGCAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggcc                                                                   >  2:199717/1‑248 (MQ=255)                                                                                                                                                                                     aaCTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGc                                                             >  1:216701/1‑251 (MQ=255)                                                                                                                                                                                            caATCTAGATGAAA‑GCAAAACGTAACAGCACCACGGGGGAGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCAGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGAGTCCCGTCCGGAGGCACTgcgc                                                      >  1:231342/1‑251 (MQ=255)                                                                                                                                                                                                  aGCTGAAA‑GCAAAACGTAACAGAACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGTTGCCATACCGGAGGATGATGCCCGTCCGGAGGCACTGCGCAGtttt                                                >  2:206302/1‑251 (MQ=255)                                                                                                                                                                                                   gCTGAAA‑GCATAACGTACCGGCACCACGGTGGGGGTAAACGCCCTGGCCTAAGAAAAGCCGGTGGAAGCCGGGGGTTACTGATTGGCCTTTGAGTACGGTCATTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTg                                               <  2:94410/251‑1 (MQ=255)                                                                                                                                                                                                    cTGAAA‑GCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCAGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgctgcgtcgtcttg                                               >  1:41709/1‑235 (MQ=255)                                                                                                                                                                                                             aaaaCGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGAtgg                                      >  1:132927/1‑251 (MQ=255)                                                                                                                                                                                                                   tAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGCTGGTGGa                                  >  2:160352/1‑249 (MQ=255)                                                                                                                                                                                                                   tAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGACCTGATGGTGGa                                  >  2:70744/1‑249 (MQ=255)                                                                                                                                                                                                                    aaCAGCACCACGGGGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAaga                               <  2:124276/251‑1 (MQ=255)                                                                                                                                                                                                                                  ggtggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGa                                  >  2:37208/1‑234 (MQ=255)                                                                                                                                                                                                                                  ggtggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGa                                  <  1:37208/234‑1 (MQ=255)                                                                                                                                                                                                                                  ggtggtgAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAAcg                 <  1:77606/251‑1 (MQ=255)                                                                                                                                                                                                                                                   ccTCAGAACATCCGCATGCAGCCGGGAGTTACAGCATGGACGTTGAGTACGGCCCGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGACGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGGGGcac  >  2:268057/1‑249 (MQ=255)                                                                                                                                                                                                                                                   ccTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac  >  1:41947/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                         GCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCAC  >  NC_000913/1428923‑1429412

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