Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,943,869	0	A	G	50.0%	35.2 / 14.0	37	D178G (GAT→GGT)	csdA	cysteine sulfinate desulfinase
Reads supporting (aligned to +/- strand):  ref base A (2/16);  new base G (18/0);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 4.19e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGAC  >  NC_000913/2943731‑2944016                                                                                                                                           |                                                                                                                                                    aaaCTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGTCCCGGTCTGGcgcggg                                                                                                                                          >  2:459643/1‑148 (MQ=255)       gAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTAc                                                                                                                                   <  1:113420/151‑1 (MQ=255)          ccAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGGCTTGCGCGGGCGGTTTCCtt                                                                                                                                >  1:102205/1‑151 (MQ=255)           cAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCCGGTCTGGcgcggggggttaacct                                                                                                                                >  1:539459/1‑138 (MQ=255)               gtggtgAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGTCTGGCGCGGGCGGTTACCttt                                                                                                                               >  2:27292/1‑147 (MQ=255)                tggtgAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGTCTGGCGCGGGGGGTTACCttt                                                                                                                               >  2:231060/1‑146 (MQ=255)                 ggtgAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGGCTGGCGCGGGCGATTACCtttttttat                                                                                                                         >  2:123159/1‑145 (MQ=255)                       aTTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGcc                                                                                                                   <  1:158068/151‑1 (MQ=255)                             gCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGAtg                                                                                                             <  2:94025/151‑1 (MQ=255)                             gCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGAtg                                                                                                             <  2:51222/151‑1 (MQ=255)                                tAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGGTCTTGCGCGGGCGATTACCTTTTCTCCTTTAGCCggggggggg                                                                                                          >  2:637422/1‑145 (MQ=255)                                     cgcAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGTCTGGCGCGGGCGGTTACCTTTGTTCATTTTGCCCGGGtggggggg                                                                                                       >  2:371570/1‑147 (MQ=255)                                                  ggATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCt                                                                                        <  2:361287/151‑1 (MQ=255)                                                        cGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAgggg                                                                                  <  1:288834/151‑1 (MQ=255)                                                            ttgttgCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTTCCCCGGTCTGGGGCGGGCGGTTACCCTTTCTCATTTAGCCGGGGGGGGGGGGGTGGGTGTTGGGGCGCCGGGGGcgg                                                                              >  1:481356/1‑149 (MQ=255)                                                             tgttgCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGGTCTGGGGCGGGCGATTACCTTTTCTCATTTTGCCGGGGGGGGGGGGGTGGGTGATGGTGCTCAgggggggg                                                                              >  2:452261/1‑147 (MQ=255)                                                              gttgCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTg                                                                            <  2:259384/151‑1 (MQ=255)                                                                   cAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCAttt                                                                       <  1:422332/151‑1 (MQ=255)                                                                       aCTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGGCCTGGCGCGGGGGATTACCTTTTTTCCTTCCGCCCGGGGGGGGGGGGGGGGTGGTGGGGGTCAGGGGGGGGGGCGTTTcccc                                                                   >  1:147776/1‑151 (MQ=255)                                                                       aCTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTcccc                                                                   <  1:568289/151‑1 (MQ=255)                                                                          gATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGGTCTGGCGGGGGCGGTTTCCCTTTCTCCATCCGCCCGGGAGGGGGGGGTGGGTGGTGGGGGCCAGGGGGGAGGGCCTTTTCCcgc                                                                 >  1:545187/1‑150 (MQ=255)                                                                          gATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCcg                                                                  <  2:477556/149‑1 (MQ=255)                                                                          gATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCcg                                                                  <  1:602891/149‑1 (MQ=255)                                                                              aCTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGTCTGGGGCGAGGGATTACCTTTTCTCCTTCCGCCGGGGTGGGGGGGGTGGGTGGTGGTGGTCCGGGGGGGGGGCATTTtccc                                                                   >  1:296200/1‑144 (MQ=255)                                                                                    cGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCGGGTCTGGGGCGGGGGGTTTCCTTTTCTCATTTAGCCGGGGTGGGGGGGGTGGGTGGTGGTGCTCAGGGGGGGGGGCATTTTCCCGCGGGTGttta                                                        >  2:65357/1‑147 (MQ=255)                                                                                       aGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGGTCTGGCGCGGGCGGTTACCTTTTCTCATTCTGCCGGGGGGGGGGGGGTGGTTGGTGGGGGTCCGGGGGGGGTGCATTTCCCCCCGGGGGTTCTGCaa                                                    >  2:298561/1‑150 (MQ=255)                                                                                          cGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTgg                                                <  2:333867/151‑1 (MQ=255)                                                                                                 tGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGGCCCGGTCTGGCGCGGGCGGGTACCCTTTTTCATTTAGCCGGGGGGGGGGGGGTGGGTGGTGGGGCTCAGGGGGGGGTGGATTTTCCCCCGGGGGGTTTGCAACTGGGTAtt                                           >  1:212424/1‑149 (MQ=255)                                                                                                  ggCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGAt                                        <  2:569549/151‑1 (MQ=255)                                                                                                        gggTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTAt                                  <  2:615698/151‑1 (MQ=255)                                                                                                          gTCAGATGTCGAACGTTACTGGCGGTTGGCCCGGTCTGGGGCGGGCGGTTACCTTTTCTCTTTTTGCCGGGGGGGGGGGGGGGGGTGGTGGGGCTCAGGGGGGGGTGGGTTTTCCCCCGGGGGTTTAGCAAATGGGTATTGAttttttt                                  >  1:15172/1‑145 (MQ=255)                                                                                                          gTCAGATGTCGAACGTTACTGGCGGTTGCCCGGGTCTGGCGCGGGCGATTACCTTTTCTCATTTAGCCGGGGGGGGGGGGGGGGGTGGTGGTGCTCAGGGGGCGGGGCATTTTCCCCCGGGTGTTCAGCAACTGGGTAtttgttttttt                                  >  1:530607/1‑140 (MQ=255)                                                                                                            cAGATGTCGAACGTTACTGGCGGTTGCCCGGCTCTCGCGCGGGCGCTTACCTTTTCTTATTCTGTCGGGGGGGGGGGGATGGGTGGTGGTGGTCCGGGGGGGGCGCATTTTCCCCCGGGGGTTCTGCCAACGGGTATTGATTTTTTTGCtt                              >  1:417131/1‑151 (MQ=255)                                                                                                            cAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCtt                              <  2:132945/151‑1 (MQ=255)                                                                                                              gATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCGTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCttt                             <  1:587643/150‑1 (MQ=255)                                                                                                                             tGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTg             >  1:566580/1‑151 (MQ=255)                                                                                                                                        cGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGAc  >  1:140946/1‑151 (MQ=255)                                                                                                                                           |                                                                                                                                                    AAACTGGAGCCAAAGTGGTGAAATTGCCGCTTAATGCGCAGCGACTGCCGGATGTCGATTTGTTGCCAGAACTGATTACTCCCCGTAGTCGGATTCTGGCGTTGGGTCAGATGTCGAACGTTACTGGCGGTTGCCCGGATCTGGCGCGAGCGATTACCTTTGCTCATTCAGCCGGGATGGTGGTGATGGTTGATGGTGCTCAGGGGGCAGTGCATTTCCCCGCGGATGTTCAGCAACTGGATATTGATTTCTATGCTTTTTCAGGTCACAAACTGTATGGCCCGAC  >  NC_000913/2943731‑2944016

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