Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,405,552	0	A	G	27.6%	41.7 / 14.8	29	E289G (GAG→GGG)	purA	adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base G (6/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGACCAGGCGCGTCAGCGTGGCGATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGA  >  NC_000913/4405313‑4405800                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         cGACCAGGCGCGTCAGCGTGGCGATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGc                                                                                                                                                                                                                                                <  2:265621/250‑1 (MQ=255)         cgcgTCAGCGTGGCGATTTCGTCATGTTTGAAGATGCTCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTCCCTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTCTGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCCAGCCgggg                                                                                                                                                                                                                                       >  2:174237/1‑250 (MQ=255)            gTCAGCGTGGCGATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAAc                                                                                                                                                                                                                                    <  1:192064/251‑1 (MQ=255)                    ggCGATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACGTATCCGTACGTATCTTCTTCCAACTCCTCTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCTCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTTCAGGTCCGGTCCCGACCGAACTGTTTGATGAAACGGTCGAGTTCCTCTGCAAGCAGAGTAACGAGTTCgg                                                                                                                                                                                                                            <  2:245461/251‑1 (MQ=255)                                aTGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCTCGTTATGTTGATTGCGTTCTGGGTATCCTCTAATCTTACTCCACTCGTGTATGGGCAGGTCCTTTCCCGACCGAGCTGTTTGATGAAACTGGCTAGTTCCTCTGCTAGCAGGGTAACGAATTCGGCGCAACTACggg                                                                                                                                                                                                                <  1:177383/251‑1 (MQ=255)                                                aGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGGATCCTCAAAGCTTACTCCACCCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCCGCAAGAAGGGTAAAGAATTTGGCGGAACAAAGGGGGGACGTCGTCGTaac                                                                                                                                                                                                >  2:4521/1‑249 (MQ=255)                                                                 atCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGt                                                                                                                                                                               >  1:262423/1‑251 (MQ=255)                                                                 atCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGt                                                                                                                                                                               <  2:194956/251‑1 (MQ=255)                                                                        aCGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGCCCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTATGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGtt                                                                                                                                                                        >  1:40953/1‑251 (MQ=255)                                                                                tATCCGTACGTAAATTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGATGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCTGGGTATCGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGc                                                                                                                                                                >  2:240125/1‑251 (MQ=255)                                                                                     gtacgtaACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTAc                                                                                                                                                           >  1:108371/1‑251 (MQ=255)                                                                                                          cacTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCAAGCAGGGTAAAGAATTCGGCGCAACTACGGGGGGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCCTCGTGCGGTAAAGCTGAACTCCCTGTCTggat                                                                                                                                      >  1:179289/1‑249 (MQ=255)                                                                                                               ctgGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCCAGCAGGGTAACg                                                                                                                                                                                                                                   >  2:178988/1‑153 (MQ=255)                                                                                                               ctgGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACGCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCAAGCAGGGTAACg                                                                                                                                                                                                                                   <  1:178988/153‑1 (MQ=255)                                                                                                                   tggcgtggcgACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTg                                                                                                                              >  1:83245/1‑250 (MQ=255)                                                                                                                          gcgACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACt                                                                                                                      <  1:74311/251‑1 (MQ=255)                                                                                                                              ccGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGAc                                                                                                                  <  2:81035/251‑1 (MQ=255)                                                                                                                                      ggCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTGCAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGa                                                                                                          <  2:259496/251‑1 (MQ=255)                                                                                                                                          tGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGc                                                                                                      >  2:146536/1‑251 (MQ=255)                                                                                                                                                                     tGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCt                                                                           >  1:254678/1‑251 (MQ=255)                                                                                                                                                                      gggTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTGCCTCTGCAAGCAGGGTGACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTgg                                                                             <  2:159567/248‑1 (MQ=255)                                                                                                                                                                          aTCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCg                                                                      >  1:74356/1‑251 (MQ=255)                                                                                                                                                                                    cTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGAt                                                            >  1:80632/1‑251 (MQ=255)                                                                                                                                                                                            aCTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTGTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGa                                                    <  1:51355/251‑1 (MQ=255)                                                                                                                                                                                               cGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGt                                                 >  1:93153/1‑251 (MQ=255)                                                                                                                                                                                                                gTTCCCGACCGAACTGTTTGATGAAACTGGCGGGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGGGCGGTACAGCCGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTTGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATTCCGGATTGGCGCGAAGTGACTACCACCCCGcggg                                >  2:206908/1‑248 (MQ=255)                                                                                                                                                                                                                   cccGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTAGGCCGTATGGTCGCGAAGTGACTACCACTCCGCTGGCa                              <  1:13014/250‑1 (MQ=255)                                                                                                                                                                                                                      gaccgaACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGAAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAATAGACTACCACTCCGATGGCAGCTg                          >  2:224769/1‑251 (MQ=255)                                                                                                                                                                                                                                              cGAGTTCCTCTGCAAGCAGGGTAACGAATTCGTCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACGGGACGTTCTGGATGGCCTGAAAGAGGTTAATCTCTGCGTGGCTTACCGTATGCCGGAGGGTCTCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGa  <  1:105104/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         CGACCAGGCGCGTCAGCGTGGCGATTTCGTCATGTTTGAAGGTGCGCAGGGTACGCTGCTGGATATCGACCACGGTACTTATCCGTACGTAACTTCTTCCAACACCACTGCTGGTGGCGTGGCGACCGGTTCCGGCCTGGGCCCGCGTTATGTTGATTACGTTCTGGGTATCCTCAAAGCTTACTCCACTCGTGTAGGTGCAGGTCCGTTCCCGACCGAACTGTTTGATGAAACTGGCGAGTTCCTCTGCAAGCAGGGTAACGAATTCGGCGCAACTACGGGGCGTCGTCGTCGTACCGGCTGGCTGGACACCGTTGCCGTTCGTCGTGCGGTACAGCTGAACTCCCTGTCTGGCTTCTGCCTGACTAAACTGGACGTTCTGGATGGCCTGAAAGAGGTTAAACTCTGCGTGGCTTACCGTATGCCGGATGGTCGCGAAGTGACTACCACTCCGCTGGCAGCTGACGACTGGAAAGGTGTAGAGCCGA  >  NC_000913/4405313‑4405800

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