Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	350,554	0	G	A	24.1%	122.7 / 11.5	59	K181K (AAG→AAA)	prpC	2‑methylcitrate synthase
Reads supporting (aligned to +/- strand):  ref base G (10/34);  new base A (14/0);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 1.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCTTACCGGCGGCGTCGCACCCGATGGATGTTATGCGCACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGTGGAA  >  NC_000913/350306‑350803                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                 gCCTTACCGGCGGCGTCGCACCCGATGGAGTTTATGCGCACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGcc                                                                                                                                                                                                                                                                <  2:204979/251‑1 (MQ=255)      cccGGCGGCGTCGCACCCGATGGATGTTATGCGCACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGGCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATGGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGc                                                                                                                                                                                                                                                           <  2:347324/250‑1 (MQ=255)                          tggatgTGATGCGCACCGGGGTTTCCGCGCTCGGCTGCACCCTGCCAGAAAAATAGGGACATACCGGCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCTCCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACATCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCAGCGGCTGTCACGTCCGGCATCTGCTGCATGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGAt                                                                                                                                                                                                                                      <  2:220574/250‑1 (MQ=255)                                  aTGCGCACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATAtctc                                                                                                                                                                                                                             <  2:133867/251‑1 (MQ=255)                                  aTGCGCACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATAtctc                                                                                                                                                                                                                             <  1:190217/251‑1 (MQ=255)                                       cACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAACAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGGCACTTCCCGCCACTGCTCCACGGC‑GAAAAAACACCGCCAAGCT‑GGG‑AAAAgcgaacgccaactccgcttg                                                                                                                                                                                                                        >  2:36888/1‑232 (MQ=255)                                         ccGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTg                                                                                                                                                                                                                      <  1:143064/251‑1 (MQ=255)                                              gtTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATGACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTa                                                                                                                                                                                                                 <  1:166054/251‑1 (MQ=255)                                                      gcTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑Aac                                                                                                                                                                                                       <  1:911/251‑1 (MQ=255)                                                           gCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCAACCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAg                                                                                                                                                                                                  <  2:248237/251‑1 (MQ=255)                                                                      ccaaaaaaaGAGGGGCATACGGTTGCTTGCTCGCGGTGTATTGCCGACAAACTGCTGGCGTCCCTGATCACGATTCTCCTTTAGGGTTAGCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCTGCGGTCACTTGCTGCAGCTGCTTCACGGC‑GAAAAGCCATCGCAAAGCT‑TGG‑AAAAGTCGATGCATACCTCGCCGGGGCTGTACGGC‑C‑AACATGAG‑TTTatggcctc                                                                                                                                                                                      <  2:142577/247‑8 (MQ=255)                                                                          aaaaaGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGCTCACCTTCTCCATATCCTGCATGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑GCAAGGCTAAGCATATCTCGCTGGTGCTGTACGCC‑ggatccccttttaccgtctccacc                                                                                                                                                                                    >  1:123890/1‑228 (MQ=255)                                                                            aaaGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCAGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACCTTCTGCCTCTGCTGCACCGC‑GGAAAGCCCACGCAAAACT‑GGG‑AAAAAGCGGTGCATATCCCGCTTGTGCTGTATCCC‑G‑CCCACGAG‑TTTAACGCCTCCCCC‑tt                                                                                                                                                                               >  2:114388/1‑251 (MQ=255)                                                                                agGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAACCTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTTCTGCCTCTGCTGCACGGC‑GAAAAACCCTCGCCAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCCGGGGCTGTACGCC‑G‑GACACGGG‑TTTAACCCCTCCCCC‑ttttcc                                                                                                                                                                           >  1:201816/1‑248 (MQ=255)                                                                                 ggggCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCGCCTCCTGCAACTGCCTCACGGC‑GAAAAACCATCGCCAAGCT‑GGG‑GAAAGGGGGTGCAAAACTCGCTGGTGCTGTACGCC‑G‑AACAAGAG‑TTTAAAGCCACAACC‑CTTAccc                                                                                                                                                                          >  2:193348/1‑250 (MQ=255)                                                                                      aTACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCCCTTCCTGCATCTGCTGCACGGC‑GAAAAACCAACGCAAAACT‑GGG‑AAAAAGCGATGCAAATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCCCCCCC‑TTTACCCGTCgg                                                                                                                                                                     >  2:148464/1‑251 (MQ=255)                                                                                      aTACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACCGC‑GAAAAACCATCGCCAAGCT‑GGG‑GAAAGGCGATGCATATCTCGCCGGGGCTGTACCCC‑C‑AACACCCG‑TTTaaacccccccccctttcccgtcgg                                                                                                                                                                      >  2:161647/1‑229 (MQ=255)                                                                                          cGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCCCTTTCTGCATCTCCCGCACGGC‑GAAAAGCCCTCGCAAAACC‑GGG‑GAAAGGCGATGCATATCTCGGTGGGGCTGGACGCCCG‑AA‑ACGAG‑GTTTACGCCTCCAccctttcccgccggggg                                                                                                                                                                   >  1:113593/1‑235 (MQ=255)                                                                                          cGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAACCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCgggggg                                                                                                                                                                 >  1:160419/1‑248 (MQ=255)                                                                                               ctGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGCCACTTCCCGCCTCTGCTGCCCGGC‑GGAAAACCATCGCCAAACT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑CCAAC‑CGAGGTTTAAC‑CCCCCCCC‑TTTACCCGTCGGGGGATTGCg                                                                                                                                                            >  2:243686/1‑251 (MQ=255)                                                                                                    gcgcgGGATATTGCCGAAAAACTGCTGGCGTCGCTTAGCTTGATTCTCCTTTATTGGAATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGCCTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGATGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCTCC‑TTTACCAGTCGGGTGATTGCGGGCAc                                                                                                                                                       <  2:190916/251‑1 (MQ=255)                                                                                                     cgcgGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCAcc                                                                                                                                                      <  1:210445/251‑1 (MQ=255)                                                                                                     cgcgGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAACCAACGCCAAGCTTGGG‑AAAAG‑CGATGCATATCTCGCTGGGGCTGTACGCC‑G‑AACACCCG‑GTTAACGCCCCCACC‑TTTTCCCGTCGGGGGATTGTGGGCAcc                                                                                                                                                      >  2:279096/1‑251 (MQ=255)                                                                                                       cgGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCgg                                                                                                                                                    <  1:205128/251‑1 (MQ=255)                                                                                                           atatTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAACCATCGCAAAACT‑GGG‑GAAAGGCGGTGCATATCTCGCTGGGGCTGTTCGCC‑C‑AACCCGCG‑TTTAACGCCCCCACC‑CTTTCCAGTCGGGTGATTGCGGGCACCGGgtct                                                                                                                                                >  2:18034/1‑251 (MQ=255)                                                                                                               tGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGat                                                                                                                                             <  1:312913/250‑1 (MQ=255)                                                                                                                 ccGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCCTCATCTTCTGCACCGC‑GAAAAAACATCCCCAAGCC‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACCCC‑G‑AACACGAG‑TTTAAAGCCTCCCCC‑TTTTCCCGGCGGGGGATTGCGGGCACCGGCTCTGacttg                                                                                                                                          >  2:208909/1‑247 (MQ=255)                                                                                                                 ccGACAAACTGCTGGCGTCGATTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTAACTTCCTTCATCTGCTGCACGGCCGCAAA‑CCATCCCAAAGCT‑GGG‑AAAAAGCGCTGCATATCTCTCTGGTGTTGGATCCC‑G‑AACACCAG‑GTTTACGCCCCCACC‑TTTACCCGTCCGGTGATTGCCGGCGCCGGCACTGGTATg                                                                                                                                          >  2:208894/1‑251 (MQ=255)                                                                                                                  cGACAAACTGCTGGCGTCGCTTAGCTCGATTCACCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAg                                                                                                                                                                                                                                           <  2:244739/157‑1 (MQ=255)                                                                                                                       agctgctgGCGTCGCTTAGCTCGATTCTCCTTTATTGTTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCTTCCGCC‑TTTACCAGTCGGGTGATTGCGGGAACCGGCTCTGATATGTATTcc                                                                                                                                    <  2:349036/249‑1 (MQ=255)                                                                                                                          tgctgGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCg                                                                                                                                                                                                                                        <  1:260401/152‑1 (MQ=255)                                                                                                                               gCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGattat                                                                                                                            <  1:189874/251‑1 (MQ=255)                                                                                                                                        aGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTTCTGCATCTGCTGCACGGC‑GGAAAACCATCGCAAAGCT‑GGG‑AAAAGGCGATGGATATCTCGCTGGGGGTGTTCGCC‑G‑AACACGAGGTTTAA‑GCCCCCACCCTTTCCAATTCGGG‑GAGTGCGGGCACCGGCTCTGGTATGTATTCCGCGATTTTTGGGGCGAt                                                                                                                   >  1:177317/1‑251 (MQ=255)                                                                                                                                                   ccTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTg                                                                                                         <  1:117975/250‑1 (MQ=255)                                                                                                                                                        aTTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCTCGATTGGCGCACTGCGCggg                                                                                                   <  2:24211/251‑1 (MQ=255)                                                                                                                                                                    aCAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACAcg                                                                                         <  1:267521/249‑1 (MQ=255)                                                                                                                                                                          acaACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCCTCCGCTGCCCGGC‑GAAAAACCATCGCAAAACT‑GGG‑AAAAGGCGGAGCATATCTCCCTGGTGCTGTACGCC‑G‑AACACGAG‑GTTAACGCCTCCACC‑CTTACCAGTCGGGTGATTGCGGGCACCCGCCCCGGTTTTGATTTCGCGATTATTGGCGCGGTTGGCGCCCTGCGCGGGCCCAAAAACGG‑Cgggggc                                                                                >  2:43251/1‑249 (MQ=255)                                                                                                                                                                           caACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGGCACATCCTCCATCTGCCGCACCGC‑GAAAAAACAACGCCAAGCT‑GGGGAAAAGG‑GATTACTATCTCGATGGTGCTGTCCGCC‑G‑AACACAAG‑GTTTACGGCACCACC‑ATTTCCCGTCGGGGGATTTCGGGCACCGGGTCTCAGAAGTATTCCCCGATTGTTGGCGCGATTGGCGCACCGGGCGGGACAAAAAACAG‑CGGcgcgg                                                                               >  2:32527/1‑250 (MQ=255)                                                                                                                                                                                gCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTTCTGCCTCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCCGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCCCCCCC‑CTTTCCAGTCGGGGGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTTGCGCGATTGGCGCACTGCGCGGGGCCAAACACGG‑Cgggggggaaga                                                                           >  2:161562/1‑243 (MQ=255)                                                                                                                                                                                  gAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACCTCCTGCCTCTGCTGCACGGC‑GGAAAACCATCGCCAAGCT‑GGG‑GAAAAGCGATGCATATCTCGCTGGGGGTGTACGCC‑G‑AAAACGAG‑TTTAACGCCTCCCCC‑CTTTCCAGTCGGGGGATTGCGGGCACCGGCCCTGATATGTATTCCGCGGTTTTTTGCGCGGTTGGCGCCCCGCGCGGGCCAAAAAACGGGCGGGGGGaaggagg                                                                         >  2:6011/1‑246 (MQ=255)                                                                                                                                                                                   aaCGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAgtg                                                                       <  1:294347/251‑1 (MQ=255)                                                                                                                                                                                    aCGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAgtgt                                                                      <  1:246131/251‑1 (MQ=255)                                                                                                                                                                                     cGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATGGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCTCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAAGGAAGTGTc                                                                     <  2:17614/251‑1 (MQ=255)                                                                                                                                                                                        aTCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCCGCATCTGCTGCACGGC‑GGAAAGCCATCGCACAGCT‑GGG‑AAAAAGCGATTCATATCTCGCTGGTGCTGTACGCC‑G‑GACACGAG‑TTTTACGCCTCCACC‑TTTACCAGGCGGGGGATTGCGGGCACCGGCTCTGGTATGTATTCCGCGATTTTTGGCGCGGTTTGCGCACCGCGCGGGCCAAAACAAGG‑CGGGGGGCATGAAGTTTCGCt                                                                  >  1:276662/1‑251 (MQ=255)                                                                                                                                                                                                ggAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATcc                                                          <  1:176658/251‑1 (MQ=255)                                                                                                                                                                                                           acgacTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTAc                                               <  1:256092/251‑1 (MQ=255)                                                                                                                                                                                                                 ccATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACg                                         <  1:84745/251‑1 (MQ=255)                                                                                                                                                                                                                        cggTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATGGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGAcg                                  <  2:287375/251‑1 (MQ=255)                                                                                                                                                                                                                            cACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGAcgaagc                              <  2:274115/251‑1 (MQ=255)                                                                                                                                                                                                                                tCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCTCGGTCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGAcgaagccga                           <  1:217305/250‑1 (MQ=255)                                                                                                                                                                                                                                      aTCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCCCTGGAGATCCAGCAACGCTATGAAACGCCGGAAGAAGCCGGAGCAGa                     >  2:260986/1‑250 (MQ=255)                                                                                                                                                                                                                                        ctgctgCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGatat                  >  1:167474/1‑251 (MQ=255)                                                                                                                                                                                                                                                     gAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGt      <  1:148464/251‑1 (MQ=255)                                                                                                                                                                                                                                                      aaaaGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGTg     <  1:210679/251‑1 (MQ=255)                                                                                                                                                                                                                                                       aaaGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGTgg    >  1:296949/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aaGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCCACGAAACGCAGGACGAAGCCGAAGCAGATATCCGCCAGCGCGTggg   >  2:32808/1‑250 (MQ=255)                                                                                                                                                                                                                                                         aGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTggcga                                                             >  2:212922/1‑189 (MQ=255)                                                                                                                                                                                                                                                         aGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGgaga                                                             <  1:212922/192‑1 (MQ=255)                                                                                                                                                                                                                                                         aGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGTGGaa  <  1:298167/251‑1 (MQ=255)                                                                                                                                                                                                                                                          gCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCCGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGTggga  >  1:249695/1‑248 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                 GCCTTACCGGCGGCGTCGCACCCGATGGATGTTATGCGCACCGGTGTTTCCGCGCTCGGCTGCACGCTGCCAGAAAAAGAGGGGCATACCGTCTCTGGCGCGCGGGATATTGCCGACAAACTGCTGGCGTCGCTTAGCTCGATTCTCCTTTATTGGTATCACTACAGCCACAACGGCGAACGCATCCAACCGGAAACCGATGACGACTCCATCGGCGGTCACTTCCTGCATCTGCTGCACGGC‑GAAAAGCCATCGCAAAGCT‑GGG‑AAAAGGCGATGCATATCTCGCTGGTGCTGTACGCC‑G‑AACACGAG‑TTTAACGCCTCCACC‑TTTACCAGTCGGGTGATTGCGGGCACCGGCTCTGATATGTATTCCGCGATTATTGGCGCGATTGGCGCACTGCGCGGGCCAAAACACGG‑CGGGGCGAATGAAGTGTCGCTGGAGATCCAGCAACGCTACGAAACGCCGGACGAAGCCGAAGCAGATATCCGCAAGCGCGTGGAA  >  NC_000913/350306‑350803

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