Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,465,167	0	G	T	29.6%	94.2 / 18.5	54	C29F (TGT→TTT)	gspM	Type II secretion system protein GspM
Reads supporting (aligned to +/- strand):  ref base G (9/29);  new base T (16/0);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 9.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAA  >  NC_000913/3464922‑3465417                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                  cGAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGc                                                                                                                                                                                                                                                             <  1:34775/251‑1 (MQ=255)         aGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCTCTCCGCGGTTTTTTTTT‑GG‑‑TTgggggg                                                                                                                                                                                                                                                    >  1:688969/1‑249 (MQ=255)           ccAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTTTCGCTCGGCG‑GGTTTTGTT‑GGG‑TTgggggggt                                                                                                                                                                                                                                                  >  1:30348/1‑247 (MQ=255)            cAGAAGCAAAAGTCGCTGCCCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGGTTTCTCGGTCGGCG‑TTTTTTTTT‑TGG‑TTGGGGGGAt                                                                                                                                                                                                                                                  >  1:775764/1‑250 (MQ=255)                   aaaaGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCTCGGCG‑GTTTTTGTT‑GGC‑TTGGGGGGGTaaaaacg                                                                                                                                                                                                                                           >  1:458498/1‑247 (MQ=255)                             gACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCt                                                                                                                                                                                                                                 <  2:748615/250‑1 (MQ=255)                               ccTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCGCCCCGCG‑GTTTTTGTT‑GGG‑TTGGGGGGGTTAAAACGACAGATAcccacc                                                                                                                                                                                                                              >  1:302568/1‑247 (MQ=255)                                    aaaaTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCa                                                                                                                                                                                                                         <  2:773809/251‑1 (MQ=255)                                       aTGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCTCCCGGCG‑TTTTTTTTTTGCC‑TTGGGG‑GATAACAACGGCAGATACCTTCATTggcagg                                                                                                                                                                                                                      >  1:401124/1‑246 (MQ=255)                                                  tagTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCTGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGTCTGTTCTCTCGCTCTGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAGCCGATAGATACCTACAtt                                                                                                                                                                                                                            <  2:591747/233‑1 (MQ=255)                                                  gagTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGCTCTCTCGGTCGGCG‑GTTTTTGTT‑TGG‑TTGGGGGGGTAAAAACGATAGGTACCTACAttaggtcgggagggaaca                                                                                                                                                                                                            >  1:591747/1‑234 (MQ=255)                                                        aGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCGGCG‑GTTTTTTTT‑TGC‑TTGGGGGGGTAAAAACGATAAAAAACTACATTGCAGAGCATCAGGCCCACgcgc                                                                                                                                                                                                      >  1:536300/1‑250 (MQ=255)                                                         gCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGGTCTCTCGCTCGGCG‑GTTTTTTTT‑TGC‑TTGGGGGGGTAAAAACACTAGAAACCTACATTGCAGAGCATCAGTCCCACGCGCa                                                                                                                                                                                                     >  2:706675/1‑250 (MQ=255)                                                               atcgatcgTGTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAAtt                                                                                                                                                                                              <  2:96276/251‑1 (MQ=255)                                                                   aacgTTTCTGTGAGTTAACACAGTCCTGGCTGCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCCCATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGCCG‑TTTTTTGTT‑GGC‑TGGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTaaga                                                                                                                                                                                          <  2:648730/249‑1 (MQ=255)                                                                            gtgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACAt                                                                                                                                                                                   <  1:363439/249‑1 (MQ=255)                                                                            gtgAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCa                                                                                                                                                                                 <  2:372040/251‑1 (MQ=255)                                                                                           ccTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAt                                                                                                                                                                   <  2:458498/250‑1 (MQ=255)                                                                                              ggCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑g                                                                                                                                                             <  1:302244/251‑1 (MQ=255)                                                                                                    cGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGAATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTTCCGGTCGGCG‑TTTTTTTTT‑GGG‑TTTGGGTGGTTAAACCGATAGATACCTTCATTGCAGAGCATCAGGCCCCCCCGCCAAAAATTAAGAAGGAAATCAAAAGGaagaagggtcaatagaga                                                                                                                                                          >  1:561401/1‑233 (MQ=255)                                                                                                       tGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTttg                                                                                                                                                                                                                                                                   <  2:131482/144‑1 (MQ=255)                                                                                                         gAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTTTT‑GGC‑TTGGGGGGATAAAAACGCTAGAAAACTACATTTCAGAGCATTAGGCCCCCGCGCCAAAAAATAAAAAAGAAATCAAAAGGAAGCAGGAG‑CAA‑GCGAGCccccc                                                                                                                                                   >  1:246036/1‑249 (MQ=255)                                                                                                                  gAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑GTTTTTTTT‑TGC‑TTGGGGGGATAAAACCGATAGAAACCTACATTTCAGAGCATCAGGCCCACCCGCAAAAAATTAAAAAGGACATCAAAAGGAAGCAAGAT‑CAA‑ACGGGCACCC‑C‑TGGGTTACTc                                                                                                                                       >  1:180954/1‑251 (MQ=255)                                                                                                                    aaaaGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTTCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTGTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAC‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCgn                                                                                                                                     <  2:339897/251‑2 (MQ=255)                                                                                                                                cAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCt                                                                                                                         <  1:210004/251‑1 (MQ=255)                                                                                                                                 aGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCtt                                                                                                                        <  1:497384/251‑1 (MQ=255)                                                                                                                                       gtATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAg                                                                                                                  <  1:384255/251‑1 (MQ=255)                                                                                                                                            gACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCAt                                                                                                             <  1:672131/251‑1 (MQ=255)                                                                                                                                               aGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCTCCCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGGGGTAAAACCCATAGATACCTACATTTCAGAGCATCAGTCCCCCGCGCCAAAAAATAAAAAGGACATCAAAAGGAAGAAGGAA‑CAA‑GAGAGCAACC‑ACGGGTTTA‑TCGGTCATCCTGCGCTTACGCACCCCCTCa                                                                                                           >  2:462767/1‑250 (MQ=255)                                                                                                                                               aGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCaa                                                                                                          <  1:257765/251‑1 (MQ=255)                                                                                                                                                 tAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCCCGGCGGTTTTTTTTT‑GG‑‑TTGGGGGGGTAAAAACCATAGATACCTACATTTCAGAGCCTCAGGCCCCCGCGCAAAAAAATAAGAAAGAAATCAAAAAGATTCAGGAA‑CAA‑ACGGGCACCCCA‑TG‑TTTACTCGGGCCTCCTCCGCTTACGCAGCCCAACaaa                                                                                                         >  1:295880/1‑250 (MQ=255)                                                                                                                                                       aCTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTTTT‑GGG‑TTGGGGGGGTTAAAACGATAGATACCTAAATTGCAGAGAAACAGTCCCCCCCGCAAAAAAATAAAAAGGAAAACAAATGGATGCAGGAG‑CAA‑AAGAGCCCCC‑C‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCCTCAAAAATA‑Tcc                                                                                                 >  1:330068/1‑251 (MQ=255)                                                                                                                                                          cAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTtc                                                                                              <  1:462767/251‑1 (MQ=255)                                                                                                                                                           aGGGAAATGATGAAATCATGGTGGGCAGAAAAATCCCCATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCATAGCATCAGTCCCACGCGCAAAAAATTAAGTAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTGATCCTGCACTTACGCAACCCAATAAGTATA‑TCCTtct                                                                                             <  2:60175/251‑1 (MQ=255)                                                                                                                                                              gAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCa                                                                                                              <  2:569126/232‑1 (MQ=255)                                                                                                                                                              gAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGATCTTCCCCCCGGCG‑TTTTTTTTT‑GGCCTTGGGG‑GGTAAAACCCATAGATACCTACATTGCAGAGCATCAGTCCCACGCCCAAAAAATTTAGAAGGAAATCCAATGGATGCCGGGT‑CAA‑GCGAGCAACC‑A‑TGGTTTACTCGGTCATCCTTCGCTTTCGCAGCCCaaaatagggagaggaaccgc                                                                                           >  1:569126/1‑232 (MQ=255)                                                                                                                                                               aaaTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGATATCTCGCTCCGCG‑TTTTTTTTT‑TGG‑TTGGGGGGATAAAAACGATATAAACCTACATTGGCGAGCATCAGTCCCAAGCGCCAAAAATTAAAAAGGACAACAAAAAGATGCAGGAT‑CAAAGC‑AGCACCC‑A‑TGGTTTACTCGGTCATCCTTCGCTTACGCAGCCCCCCAAAAAAA‑TCCTtctcc                                                                                           >  1:788260/1‑248 (MQ=255)                                                                                                                                                                aaTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCtt                                                                                                                        <  2:664194/220‑1 (MQ=255)                                                                                                                                                                 aTGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACAGTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGa                                                                                       <  2:302568/251‑1 (MQ=255)                                                                                                                                                                   gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTTTT‑GGG‑TTGGGGGGGGAAAAACGATAGATACCTACATTTCAGAGCATCAAGCCCCCGCGCAAAAAATTAAAAAAGACAACAAATGGATGCAGGGT‑CAA‑GCGAGCACCC‑C‑TGGGTTACTCCGGCATCCCTCGCTTTCGCACCCCCTCAAAAAAA‑Tccctttccaggg                                                                                       >  1:388759/1‑239 (MQ=255)                                                                                                                                                                   gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAg                                                                                     <  1:450656/251‑1 (MQ=255)                                                                                                                                                                   gATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAg                                                                                     <  1:33223/251‑1 (MQ=255)                                                                                                                                                                        aaTCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCCGCG‑GTTTTTTTT‑TGC‑TTTGGGGGGTAAAAACCATAGATACCTACATTTCAGAGCATCAGTCCCCCGCGCAAAAAAATAAAAAAGACATCAAAAGGATGCCGGGTACAA‑‑CGAGCACCC‑A‑TGGTTTACTCGGTGATCCTTCCGTTACGCCGCCCCCCACAAAAA‑TCCTTTTCCAGGTAGCaaa                                                                                 >  2:531763/1‑250 (MQ=255)                                                                                                                                                                              ggtggGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATTGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑GGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTgag                                                                          <  2:561401/251‑1 (MQ=255)                                                                                                                                                                                 gggCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTTTCCCGCTCGGGG‑TTTTTTGGT‑GGC‑TTGGGGGGGTAAAAACGATAGATACCTACATTGCGAAGAATCAGGCCCACGCGCAAAAAATTAAAAAGGACATCAAAAGGGTGCAGGGT‑CAA‑GCGAGCACCC‑A‑TGGTTTTCTCGGGTATCCTGCGCTTTCCCAGCACATCAAAAATA‑TCCTTCTCGAGGAAACCAAACGTGAGaa                                                                        >  1:585896/1‑250 (MQ=255)                                                                                                                                                                                          aaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTTCCTCCGCG‑CTTTTTTTT‑GGC‑TTGGGGGGGTAAAAACGAAAAAAAACTACATTGTAGAGCATCAGGCCCCCGCGCCAAAAAATAAAAAGGAAATCAAAAGGATAGAAGAT‑CAA‑GCGAGCACCC‑C‑TTGGTTACCCCGGCATCCCGCGCCTACGCCACCCATCAAACAAA‑TCCATCTCCGGGGAGCAAAAcaagggaaattgggtacc                                                              >  2:716045/1‑234 (MQ=255)                                                                                                                                                                                                       aaaaCAGATCGTTGCGGCACTGGCGGTTCTCTCCCCCGGCG‑TTTTTTTTT‑GGG‑TTGGGGGGGTAAAAACGATAGGAACCTCCATTGCAGCGGATCAGTCCCACGCCCAAAAAATTAAAAAAGACCACAAATTGATGCAGGAT‑CAA‑CCGAGCACCC‑C‑TTGTTTAATCTGTCATCCTTCCCTTACGCAGCCCCTCAAAAATAATACT‑CTCCTGGTAGCAAAACGTGCGAATTTTGCtttctccccggtaaac                                                  >  2:268985/1‑235 (MQ=255)                                                                                                                                                                                                            aGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCg                                            <  1:531763/251‑1 (MQ=255)                                                                                                                                                                                                                    gCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATa                                     <  2:180954/250‑1 (MQ=255)                                                                                                                                                                                                                     cGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATAct                                   <  1:233431/251‑1 (MQ=255)                                                                                                                                                                                                                       gCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATActct                                 <  2:251487/251‑1 (MQ=255)                                                                                                                                                                                                                                 ttCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCAt                       >  1:659325/1‑251 (MQ=255)                                                                                                                                                                                                                                             gCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACg                                                 <  2:663679/213‑1 (MQ=255)                                                                                                                                                                                                                                             gCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACg                                                 >  1:663679/1‑213 (MQ=255)                                                                                                                                                                                                                                                 ttttttGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTccc        >  1:334770/1‑251 (MQ=255)                                                                                                                                                                                                                                                       gtt‑gGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGaa  <  2:388759/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                  CGAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCG‑TTTTTTGTT‑GGC‑TTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGAT‑CAA‑GCGAGCACCC‑A‑TGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATA‑TCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCACGCTGGAAAACGGCCCGGACAATACTCTGACAATTCATCCTGTCACAGCTCCCCTCGAA  >  NC_000913/3464922‑3465417

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