Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,171	0	A	G	55.6%	24.5 / 75.7	63	E41E (GAA→GAG)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (20/8);  new base G (18/17);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAG  >  NC_000913/1428926‑1429414                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    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                                    <  1:37208/234‑1 (MQ=255)                                                                                                                                                                                                                               ggtggtgAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGa                                    >  2:37208/1‑234 (MQ=255)                                                                                                                                                                                                                               ggtggtgAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAAcg                   <  1:77606/251‑1 (MQ=255)                                                                                                                                                                                                                                         ccggtgggcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGt         <  1:225220/243‑1 (MQ=255)                                                                                                                                                                                                                                                gcTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTGTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAACGATTTTCGCTGTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCATCGTCTTGAGCTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAg  <  1:37623/250‑1 (MQ=255)                                                                                                                                                                                                                                                ccTCAGAACATCCGCATGCAGCCGGGAGTTACAGCATGGACGTTGAGTACGGCCCGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGACGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGGGGcac    >  2:268057/1‑249 (MQ=255)                                                                                                                                                                                                                                                ccTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGcac    >  1:41947/1‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    CCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAA‑GCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGGTGGAAGAGGTGGCGCGTAACGCGTCCGCGGTGGCACAG  >  NC_000913/1428926‑1429414

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