Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,429,708	0	A	G	57.8%	27.7 / 78.9	65	T220T (ACA→ACG)	stfR	putative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (19/18);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAA  >  NC_000913/1429465‑1429943                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            tGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCAt                                                                                                                                                                                                                                      >  1:160349/1‑251 (MQ=255)          ggcccttgtgactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACcgcgg                                                                                                                                                                                                                             >  1:111836/23‑250 (MQ=255)          ggcccttgtgactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACcgcgg                                                                                                                                                                                                                             <  1:222138/229‑2 (MQ=255)             ccgtgtgactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCa                                                                                                                                                                                                                          <  1:80344/232‑1 (MQ=255)              cATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCAc                                                                                                                                                                                                                         >  2:140376/1‑251 (MQ=255)                 gtgactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGa                                                                                                                                                                                                                       >  1:100070/16‑250 (MQ=255)                   gactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGaa                                                                                                                                                                                                                      >  1:78670/14‑249 (MQ=255)                   gactgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGaa                                                                                                                                                                                                                      >  1:195212/14‑249 (MQ=255)                    actgatgcaactgACTCAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGaaagc                                                                                                                                                                                                                   >  1:132306/13‑248 (MQ=255)                               ggACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTg                                                                                                                                                                                                        >  2:80076/1‑251 (MQ=255)                                       cACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCg                                                                                                                                                                                                >  2:196311/1‑251 (MQ=255)                                         cgcgCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGaggccgccacttcagc                                                                                                                                                                                              <  1:6870/251‑16 (MQ=255)                                           cgcAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGccc                                                                                                                                                                                            <  1:143722/251‑1 (MQ=255)                                                         cAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCggcggc                                                                                                                                                                               <  2:11943/250‑1 (MQ=255)                                                          aGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCt                                                                                                                                                                             <  1:179884/251‑1 (MQ=255)                                                            ccGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTc                                                                                                                                                                            >  1:113242/1‑250 (MQ=255)                                                              ggACAGGCTGCATCGTCAGCTCAGGAAGCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCaaa                                                                                                                                                                         <  1:138029/251‑1 (MQ=255)                                                                         atcgtcaGCTCAGGAAGCGTCCGCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAACCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCaa                                                                                                                                                              <  2:79325/250‑1 (MQ=255)                                                                             tcagctcaggaagCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatc                                                                                                                                                          >  2:48083/12‑251 (MQ=255)                                                                              cagctaaggaagCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAatc                                                                                                                                                          >  2:211076/11‑250 (MQ=255)                                                                                  tcagtcagCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAatcatc                                                                                                                                                       >  2:190421/1‑249 (MQ=255)                                                                                     gtcagCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGaaa                                                                                                                                                  <  1:121172/251‑1 (MQ=255)                                                                                     ggaagCGTCCTCCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGaaa                                                                                                                                                  <  2:111836/248‑1 (MQ=255)                                                                                        agCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCCGAaacg                                                                                                                                                >  2:4886/1‑250 (MQ=255)                                                                                          cGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAaacgaac                                                                                                                                             >  1:172334/1‑251 (MQ=255)                                                                                               tCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGca                                                                                                                                          <  2:199121/249‑1 (MQ=255)                                                                                                  aGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCatcatc                                                                                                                                     <  1:80076/251‑1 (MQ=255)                                                                                                       atacgcGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGATTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGGGCGGCGAAAACGTCAGAAACGAATGCGGCAGCGTCACAGCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAgtg                                                                                                                                <  2:172436/246‑3 (MQ=255)                                                                                                               cATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCagtagg                                                                                                                        >  2:9807/1‑250 (MQ=255)                                                                                                                     aaaaGGCCACTGAAGCGGAAAAAAGGGCCGCAGCCGAAGAGTCCTCAAAAAACGCGGCGGCCACCAGGGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCt                                                                                                                  <  2:73099/250‑1 (MQ=255)                                                                                                                      aaaGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCtt                                                                                                                 >  1:79246/1‑251 (MQ=255)                                                                                                                       aaGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTc                                                                                                                <  1:60967/251‑1 (MQ=255)                                                                                                                         ggCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCt                                                                                                              >  1:7525/1‑251 (MQ=255)                                                                                                                          gCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCCCCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTc                                                                                                             <  1:67775/251‑1 (MQ=255)                                                                                                                                 aaGCAACAAAAAGTGCGGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAAcg                                                                                                      <  2:81363/251‑1 (MQ=255)                                                                                                                                       cAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACggcggc                                                                                                 <  1:180651/250‑1 (MQ=255)                                                                                                                                        aaaaaGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGGAGTGCAGCTTCCTCGGCAACggcggccg                                                                                               >  2:31921/1‑249 (MQ=255)                                                                                                                                           aaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACggcggccgaa                                                                                             >  2:183677/1‑246 (MQ=255)                                                                                                                                             gTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCACGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGaaa                                                                                           >  2:130073/1‑250 (MQ=255)                                                                                                                                               gccgcagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTc                                                                                        >  1:129975/10‑251 (MQ=255)                                                                                                                                                ccgcagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTct                                                                                       <  2:100070/243‑2 (MQ=255)                                                                                                                                                 cgcagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctg                                                                                      <  2:3378/244‑3 (MQ=255)                                                                                                                                                 cgcagccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCCGCCACGAAAGCGTCAGAAGCCGCCACTTCAGCACGAGATGCGGAGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGGGCCGTTCGGGCAGCTTCCTCGGCAACGGTCGCAGAAAAttacc                                                                                      >  2:235487/8‑248 (MQ=255)                                                                                                                                                    acccggagagTCCGCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGGCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCTGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgcca                                                                                   <  1:233262/245‑6 (MQ=255)                                                                                                                                                        gCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctgccagggc                                                                               >  2:128970/1‑242 (MQ=255)                                                                                                                                                            agTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTCGTGCAGCTTCCTCGGCCACGGCGGCAGGCAATTCCGCGAAggcggc                                                                            >  2:203467/1‑250 (MQ=255)                                                                                                                                                                  cAAAAAACGCGGCGGCCACCAGGGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCaaaaa                                                                       <  1:227330/249‑1 (MQ=255)                                                                                                                                                                     aaaGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTcc                                                                  <  1:59776/251‑1 (MQ=255)                                                                                                                                                                                          ccGGGGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGtcatct                                             <  2:33972/251‑1 (MQ=255)                                                                                                                                                                                           cGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTg                                            >  1:35670/1‑251 (MQ=255)                                                                                                                                                                                              tGCGGCGAAAACGTCAGAAACGAACGCTTCAGCGTCACTACAATCAGCAGCAACCTCTGCATACACCGCGACCACGAAGGCATCAGAAGCTGCGAACTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGaa                                          >  1:215935/1‑250 (MQ=255)                                                                                                                                                                                                 ggcgAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAc                                        >  1:54264/1‑249 (MQ=255)                                                                                                                                                                                                  gcgAACACGTCAGAAACAACTGCTGCAGAGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCCGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcagc                                     >  1:208110/1‑248 (MQ=255)                                                                                                                                                                                                          gTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGCATTCTGAAACAGAAGCggggc                             >  2:73191/1‑249 (MQ=255)                                                                                                                                                                                                          gTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaac                             <  1:12630/251‑4 (MQ=255)                                                                                                                                                                                                           tCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCATCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaacg                            >  1:215142/1‑247 (MQ=255)                                                                                                                                                                                                            cAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaacg                            <  1:248645/250‑5 (MQ=255)                                                                                                                                                                                                              gAAACGAATGCTTCCGAGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACag                           >  1:47675/1‑249 (MQ=255)                                                                                                                                                                                                              gAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaacgg                           >  1:182880/1‑244 (MQ=255)                                                                                                                                                                                                                      ggctgcAGCGTCACAAAAAACAGCCGCCACGTCTGCCTCCCCCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTctg                 <  2:197806/246‑3 (MQ=255)                                                                                                                                                                                                                       gCTTCAGCGTAACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCggg                >  1:45424/1‑250 (MQ=255)                                                                                                                                                                                                                           cAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCg            <  2:201249/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGccgccgcc   >  2:65341/2‑244 (MQ=255)                                                                                                                                                                                                                                     aCAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTc    <  1:178087/249‑1 (MQ=255)                                                                                                                                                                                                                                      cAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAgacgcaa  >  1:74599/1‑244 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            TGAGGCGGCAACCCATGCGGCTGATGCTGCGGACTCAGCACGCGCAGCCAGCACGTCAGCCGGACAGGCCGCGTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAA  >  NC_000913/1429465‑1429943

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