Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,100,265	0	T	A	37.5%	23.2 / 32.3	32	Q335L (CAG→CTG)	gnd	6‑phosphogluconate dehydrogenase, decarboxylating
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base A (6/6);  total (16/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCGGCAACACGCTGATCTTTCAGAGAAGAGAT  >  NC_000913/2100122‑2100404                                                                                                                                                |                                                                                                                                            gtgGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCg                                                                                                                                        >  6:65180/1‑149 (MQ=255)            gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCGGCTGAGAGAAGCCCTg                                                                                                                                            >  8:150512/1‑134 (MQ=255)            gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCCGCTGAGAGAAGCCCTg                                                                                                                                            <  7:150512/134‑1 (MQ=255)            gCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGAt                                                                                                                             >  5:22304/1‑149 (MQ=255)                   cATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCTGCTGAGAGAAGCCCTgg                                                                                                                                           <  1:61883/128‑1 (MQ=255)                   cATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTgg                                                                                                                                           >  2:61883/1‑128 (MQ=255)                             tttCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaaa                                                                                                                                 >  7:93452/1‑128 (MQ=255)                             tttCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGaaa                                                                                                                                 <  8:93452/128‑1 (MQ=255)                                                 ggATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGAt                                                                                         >  3:67040/1‑148 (MQ=255)                                                    tgatgCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAAc                                                                                     <  6:22304/149‑1 (MQ=255)                                                                 aCGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGc                                                                        <  7:113906/149‑1 (MQ=255)                                                                 aCGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGc                                                                        >  4:83110/1‑149 (MQ=255)                                                                   ggAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTc                                                                          >  1:113340/1‑145 (MQ=255)                                                                   ggAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTc                                                                          <  2:113340/145‑1 (MQ=255)                                                                    gAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCa                                                                                                          >  1:81708/1‑112 (MQ=255)                                                                    gAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCa                                                                                                          <  2:81708/112‑1 (MQ=255)                                                                    gAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACACAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCtt                                                                             >  8:106475/1‑141 (MQ=255)                                                                    gAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACACAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCtt                                                                             <  7:106475/141‑1 (MQ=255)                                                                                 aTTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTg                                                        >  7:132711/1‑149 (MQ=255)                                                                                       ccGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGAcc                                                  <  4:67040/149‑1 (MQ=255)                                                                                        cGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCa                                                 <  4:63441/149‑1 (MQ=255)                                                                                          tAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCaga                                               <  5:65180/149‑1 (MQ=255)                                                                                               cAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAgagaga                                           <  3:83110/148‑1 (MQ=255)                                                                                                  aTCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACtt                                       >  7:44217/1‑149 (MQ=255)                                                                                                  aTCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACtt                                       >  4:138056/1‑149 (MQ=255)                                                                                                  aTCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACtt                                       >  4:135760/1‑149 (MQ=255)                                                                                                           gTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCg                               <  3:135760/148‑1 (MQ=255)                                                                                                           gTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCAGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCg                               <  3:138056/148‑1 (MQ=255)                                                                                                               tctTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCGGCaa                           >  6:114978/1‑148 (MQ=255)                                                                                                                       cgcagcacgcagcTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCGGCAACACGCTGAt                  <  5:114978/149‑1 (MQ=255)                                                                                                                          agcacgcagcTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCGGCAACACGCTGATCtt               >  8:148426/1‑149 (MQ=255)                                                                                                                                       gagaAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCGGCAACACGCTGATCTTTCAGAGAAGAGat  <  8:44217/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                            GTGGATTTTCGGCATAAGCATCGGTGATTTTCTGCAGGAACTGCGCACGGATGATGCAGCCAGCACGGAAAATCTTCGCGATTTCGCCGTAGTTCAGATCCCAGTTGTACTCTTCAGACGCAGCACGCAGCTGAGAGAAGCCCTGGGCGTAAGAAACGATTTTGCCCAGATACAGCGCACGACGAACTTTTTCGATGAACTCAGCCTTGTCGCCTGCTGGCTGTGCTTGCGGACCAGAGAGAACTTTAGATGCGGCAACACGCTGATCTTTCAGAGAAGAGAT  >  NC_000913/2100122‑2100404

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