Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,950,173	0	A	G	24.5%	101.0 / 11.0	53	E311G (GAG→GGG)	argA	N‑acetylglutamate synthase
Reads supporting (aligned to +/- strand):  ref base A (9/31);  new base G (13/0);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 5.91e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTT  >  NC_000913/2949923‑2950422                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                   aaCGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGTAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGa                                                                                                                                                                                                                                                                    <  1:65156/251‑1 (MQ=255)   cGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGt                                                                                                                                                                                                                                                                  <  2:468069/251‑1 (MQ=255)           aagcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACCACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCGGTTGCAAGAGGTGTTCTCACGCGACGGTATCGGTACGCAGAGTGTGGTGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATAGTGGCGGAATTCT‑CGAGA‑TGATTCGc                                                                                                                                                                                                                                                         >  1:16615/1‑251 (MQ=255)            agcgcGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGcc                                                                                                                                                                                                                                                        <  1:283670/251‑1 (MQ=255)              cgcgGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑a                                                                                                                                                                                                                                                     <  2:418476/251‑1 (MQ=255)                 ggggaGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGGCATTTACTCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGAGTAT‑GGAGT‑TGATTTGCCC‑ACtc                                                                                                                                                                                                                                                  >  2:85748/5‑250 (MQ=255)                 gggTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGGGG‑TGATTTGCCC‑ACTg                                                                                                                                                                                                                                                  >  1:132493/1‑251 (MQ=255)                   gTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCCGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGCTATTGGCGGGTTTCT‑GGAGG‑TGTTTCgtacgccggc                                                                                                                                                                                                                                                 >  2:215387/1‑242 (MQ=255)                       aaGCCCAGGAAGAGAAAGGCTATTCCAACCCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGagcag                                                                                                                                                                                                                                            <  2:225305/251‑1 (MQ=255)                               gAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTAt                                                                                                                                                                                                                                    <  1:199989/251‑1 (MQ=255)                                            aTTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGGGG‑TGGTTTGCCC‑CCCTGAGCAGCGAGGTATTATGG‑TAcacccg                                                                                                                                                                                                                      >  1:471456/1‑246 (MQ=255)                                                     ccGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTT‑GGGGT‑TGATTTGCCCCAC‑GGGGCAGCAAGGTAATCttggcagcccgttctccggg                                                                                                                                                                                                              >  1:332783/1‑232 (MQ=255)                                                        gTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑TGAGG‑TTATTCGCCC‑CCCGCAGCAGACAGGGATTCTTG‑TACGCCCGTCTCGCGagctg                                                                                                                                                                                                          >  1:176459/1‑249 (MQ=255)                                                                       tGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAt                                                                                                                                                                                           <  1:353071/250‑1 (MQ=255)                                                                         cGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGCAATTT‑GGGGG‑TGGTTCGCCC‑ACCGGGGCAGCAAGGGATTCTGG‑GACGCCGTTTTCGCGAGCAGCTGGGG‑ATGGGAAT‑Cg                                                                                                                                                                                        >  1:53035/1‑250 (MQ=255)                                                                             gcgcAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCCTTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAACCAACGATATTGGCGGAATTCT‑GGGGG‑TGATTTGCCC‑ACCGGGGCAAGAAGGGATTTTTG‑GACGGCGGTTTCGCGAGCCGCTGGGG‑ATGGGAAT‑CGACaa                                                                                                                                                                                    >  1:482360/1‑250 (MQ=255)                                                                               gcAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAAtt                                                                                                                                                                                 <  2:282483/251‑1 (MQ=255)                                                                                      aaGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTT‑GGGGG‑TGGTTCGCCC‑CCCGGAGCCGCCAGGGATTTTTG‑GAAGCCGGTTTCGCGGGCAAGTGGTG‑GTGGAAAT‑CGTCAAATTCAccctt                                                                                                                                                                          >  1:389503/1‑248 (MQ=255)                                                                                         ccggCCGCATCGGCGGGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACTCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGAGTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCattatt                                                                                                                                                                       <  2:178598/248‑1 (MQ=255)                                                                                              cGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAgcg                                                                                                                                                                  <  2:168947/251‑1 (MQ=255)                                                                                                      cgtgcgtCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTACCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCT‑GGGGTGTG‑TTCGCCC‑CCTGGGGCAGCAAGGGATTTTGG‑TACGCCGTTCTCGCGAGCCGCCGGGG‑GTGGAAATAGGACAA‑TTCACCATTATTCACCGCGA‑Taaa                                                                                                                                                          >  2:408324/1‑249 (MQ=255)                                                                                                       gtgcgtCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTTT‑GGGGT‑TGATTCGCCC‑ACTGGAGCAGACAGGGATTCTGG‑TACCCCGTTCTCCCGGGCGGCTGGGG‑GGGGAAAT‑CGACAAATTAACCATTTTTCAgcggtgaaaca                                                                                                                                                          >  1:416580/1‑242 (MQ=255)                                                                                                             cGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTa                                                                                                                                                  <  1:327946/250‑1 (MQ=255)                                                                                                                      ttAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTgcg                                                                                                                                         <  2:287738/250‑1 (MQ=255)                                                                                                                       tAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTAATCT‑GGGGG‑TGGTTAGCCC‑ACTGGGGCAGCAAGGGATTCTGG‑GAAGCCGTTCTCGCGGGGAGCCGGGG‑AGGGAAAA‑CGGAAAAAACACCATTATTCAGCGCGA‑TAACA‑AGACTATTGCCTGggcc                                                                                                                                       >  1:280854/1‑251 (MQ=255)                                                                                                                                aTCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTAt                                                                                                                              <  2:389503/251‑1 (MQ=255)                                                                                                                                aTCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTAt                                                                                                                              <  2:207705/251‑1 (MQ=255)                                                                                                                                aTCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTCT‑GGAGT‑TGATTCGCCC‑CCTTGCGCAGCACGGTATTCTGG‑TACGCCGTTCTCGAGAGCTGATGGGG‑ATGGAAAT‑AGACAAATTCAACATTATTTGTCGCGAATAAC‑‑CGACTATTTCCTGCGCCGCGcttcac                                                                                                                              >  1:409589/1‑247 (MQ=255)                                                                                                                                   aGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGGATTTT‑GGGGT‑TGATTCGCCC‑ACCGGAGCAGAAAGGTATTCTGG‑AACGCCGGTCTCGCGCGCAGCTGGAG‑ATGGGAAT‑CGAAAAAATCACCATTATTTAGCGCGA‑TAAAC‑CGGCAATTGCCTGATTCGCGCTCTAtacc                                                                                                                           >  1:380908/1‑248 (MQ=255)                                                                                                                                    ggAAGATGGCGCGCTGTTGCAAGAGTTGTTCGTTCGCGCCGGTATCGGTACGCTGATTGGTATGGAAAGCGCCGAGCTGATTCGTCGCGCAACAATCAACGATATTGTCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCTTTATTCAGCGCGA‑TCATA‑CGACTATTTCCTGCGCAGCGCTCTTTCCGt                                                                                                                          <  2:317771/251‑1 (MQ=255)                                                                                                                                            gcgcgcGGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaa                                                                                                                  <  2:273512/251‑1 (MQ=255)                                                                                                                                              gcgcTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGAGTCGTCGCGCAACAATCAACGATATTGGCTGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAG‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGgaaga                                                                                                                <  1:329045/251‑1 (MQ=255)                                                                                                                                                     tGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGAt                                                                                                          <  1:469703/250‑1 (MQ=255)                                                                                                                                                                           cGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGAGTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAg                                                                                 <  1:308227/251‑1 (MQ=255)                                                                                                                                                                               aTCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTc                                                                              <  2:438058/250‑1 (MQ=255)                                                                                                                                                                                  ggTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGGGT‑TGGTTTGCCC‑ACCGGAGCCGCAAGGGATTCTGG‑TACGCCGTTCTCGCGAGCAACTGGAG‑ATGGAAAT‑AGACAAATTCACCATTATTCAACGCGA‑TAAAA‑AGGCCATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGGAGGTTGGGGAAATGGCC‑TGTGTGG‑CGGTTCAccc                                                                          >  1:402709/1‑251 (MQ=255)                                                                                                                                                                                         aGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTAc                                                                   <  1:96840/251‑1 (MQ=255)                                                                                                                                                                                                 agggAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTtc                                                           <  2:547829/249‑1 (MQ=255)                                                                                                                                                                                                   ggAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGGTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTtca                                                          <  2:145616/250‑1 (MQ=255)                                                                                                                                                                                                    gAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCTTGGGGT‑T‑ATTTGCCC‑CCTGGGGCAGCAAGGGATTTTGGGAACGC‑GTTCTCGCCAGCCGCTGGAGGAGGGAAA‑‑CGGAACAATCACCATTATTTAGCGCGG‑TAAAAACAAC‑ATTGGCCGCGGCGCGGCCTATCCGTTCCCGGGAGGGAAGGTTGGGGGAAAGGCC‑T‑TGTGGGCAGTTCACCCggggttccacaggtcatc                                                        >  2:494698/1‑235 (MQ=255)                                                                                                                                                                                                         cgcCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTATCCC‑CCTGGAGCAGCCAAGTATTCTTG‑GACGCCGTTCTCGCGGGCCGCTGGGG‑ATGGGAAA‑CGACCAATTCACCCTTATTCAGCGCGA‑TAACA‑CGCATATTTCCTCCGCCGCGCCCTTTCCCGTCCCGGAAGAGCAGATTGGGGGAATGGCC‑TGTGTGG‑CAGTTCCCCCCGATTACCCCAgcttcacacggg                                                   >  2:101354/1‑240 (MQ=255)                                                                                                                                                                                                          gcCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGGGT‑TGATTTGCCC‑ACTGGGGCAGCCAGGGATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGGG‑ATGGGAAT‑CGACAAATTTTCCCTTTTTCAGCGCGG‑TAACA‑CGACTATTGCCTGCCCCGCGCCCTATCCCTTCCCGGGAGAGGAGATTGGGGGAATGGCCCTGGGGGG‑C‑GTTTACCCGGGATACCGCCGTTCATCACgggg                                                  >  2:315581/1‑251 (MQ=255)                                                                                                                                                                                                           ccGAGCAGATTCGGCGCGCAACAATCAACGATACTGGCGGTATGCT‑GGAGG‑TGATTCGCCC‑ATTGGAGCAGCAAGGTATTCTGG‑TACGCCGCTCTCGCGAGCAGCGGGAG‑ATGGAAAT‑CGACAAATCCATCATTAGTCAGCGCGA‑TAACG‑CGACTACTGCCTGCGTCGCGCTCTAGCTGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGCTCACCCGGATTACCGCAGATCATGAAGGGGt                                                 <  2:279588/251‑1 (MQ=255)                                                                                                                                                                                                             gAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGa                                               <  1:315581/251‑1 (MQ=255)                                                                                                                                                                                                                      cgtcgCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgc                                       <  1:14777/250‑1 (MQ=255)                                                                                                                                                                                                                      cgtcgCGCAACAATCAACGATATTGGCGGTAGTCT‑GGAGG‑TGAGTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTGCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATGCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTGCTCCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTctgc                                       <  2:356737/250‑1 (MQ=255)                                                                                                                                                                                                                            gcAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAATT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACg                                <  2:371502/251‑1 (MQ=255)                                                                                                                                                                                                                                        gATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGTCAAATTCCCCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTTCCGCAGTTCATCAAGGGGTTAAGTTCTGCTGGAACGCATTGCCGCTCa                    <  2:179812/251‑1 (MQ=255)                                                                                                                                                                                                                                           agtGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGc                 <  2:132493/249‑1 (MQ=255)                                                                                                                                                                                                                                               gCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAACTGGCC‑TGTGTGG‑CAGTTCACCCGGCTTACCGCCGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCCTTGCCGCTCAGGCGa               >  2:426166/1‑249 (MQ=255)                                                                                                                                                                                                                                                 ggTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTACATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCa           <  1:512701/251‑1 (MQ=255)                                                                                                                                                                                                                                                 ggTATTCG‑TGAGG‑TGATTCGGCC‑GCTGGAGCAGGAAGGTATTCTGG‑TATGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAAGTCGTCATTATTCAGCGGGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTAGCCGGTCCCGGAAGAGAAGATTGGGGAAGTGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCATGGGGTGAAGGTCTTCTTGAATGCATTGCCGCTCAGGCGATGCa           <  2:71875/251‑1 (MQ=255)                                                                                                                                                                                                                                                   tATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCTCTCAGGCGAAGc            >  2:100170/1‑248 (MQ=255)                                                                                                                                                                                                                                                   tATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGc            >  1:536664/1‑248 (MQ=255)                                                                                                                                                                                                                                                            aGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAACTGGCC‑TGTGTGG‑CAGTTCACCCGGCTTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCCGGCGAAGCAGAGCGGCtt  >  2:484313/1‑250 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                   AACGAAGCGCAAGCGCGGGTAGAAGCCCAGGAAGAGAAAGGCGATTACAACTCCGGTACGGTGCGCTTTTTGCGTGGCGCAGTGAAAGCCTGCCGCAGCGGCGTGCGTCGCTGTCATTTAATCAGTTATCAGGAAGATGGCGCGCTGTTGCAAGAGTTGTTCTCACGCGACGGTATCGGTACGCAGATTGTGATGGAAAGCGCCGAGCAGATTCGTCGCGCAACAATCAACGATATTGGCGGTATTCT‑GGAGT‑TGATTCGCCC‑ACTGGAGCAGCAAGGTATTCTGG‑TACGCCGTTCTCGCGAGCAGCTGGAG‑ATGGAAAT‑CGACAAATTCACCATTATTCAGCGCGA‑TAACA‑CGACTATTGCCTGCGCCGCGCTCTATCCGTTCCCGGAAGAGAAGATTGGGGAAATGGCC‑TGTGTGG‑CAGTTCACCCGGATTACCGCAGTTCATCAAGGGGTGAAGTTCTGCTGGAACGCATTGCCGCTCAGGCGAAGCAGAGCGGCTT  >  NC_000913/2949923‑2950422

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