Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,177,531	0	A	C	58.6%	14.0 / 28.5	30	V127G (GTG→GGG)	nudF	ADP‑ribose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base C (0/17);  total (5/25)
Fisher's exact test for biased strand distribution p-value = 6.67e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-02

TGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAA  >  NC_000913/3177383‑3177678                                                                                                                                                     |                                                                                                                                                    tGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAacac                                                                                                                                                    >  1:104938/1‑150 (MQ=255)     tGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGtt                                                                                                                                               <  1:283070/151‑1 (MQ=255)        tCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTgg                                                                                                                                            >  1:627029/1‑151 (MQ=255)                           cgcgAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTcc                                                                                                                         <  2:297027/151‑1 (MQ=255)                                       ctttttCATCAGCCAGCCCGTGAATCCCGCTTGCGGTCGTGGCGCCCACTCCCCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAAt                                                                                                             <  1:671036/149‑1 (MQ=255)                                       cGTTTCCACCAGCCAGACCGTGAATCCCGCTTGCGGTCGTGGCCCCCCCTCCCCCCACCATAATTGACGAACGCTCCCTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAAt                                                                                                             <  1:67555/151‑1 (MQ=255)                                                  cccAGACCGTAAATCCCCCTTGGGGTCGTGGCGCCCACTCCCCCCACCATAATTGACGAACGCTCACTGGGGCCCCCCGGGCTTCCCAGGAAACTTACCCCCGGTTTGTTCGGTTTGACTACCAGTCCCGCCTCTTCAATCGCTTCGCGAc                                                                                                  <  1:379693/150‑1 (MQ=255)                                                      gACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACggg                                                                                               >  2:422144/1‑150 (MQ=255)                                                             aaTCCCGCTTGCGGTCGTGGCGCCCCCTTCGCCCCCCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCt                                                                                       <  2:179087/151‑1 (MQ=255)                                                                    cgtGCGGTCGGGGCGTCCCCTCCGCCCACCAAAATTGCCGACCCCTCCCTGGTCCCCCCCGGGCTTGCCAGAAAACTAACCCCCGGTTTGGTCCTTTTGCCTACCACTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCacac                                                                                <  1:260221/149‑1 (MQ=255)                                                                            cggggCGCCCACTCCCCCCCCCATAATTGACGACCGCTCCCTGGTCCCCCCCGGGCTTGCCAGGAAACTTAACCCCGTTTGGGTCCGTTTGACTACCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCAccc                                                                        <  1:168229/148‑1 (MQ=255)                                                                            cGTGGCGTCCACTCCGCCCACCATAATTGACGAACGCTCCCTGGTGCCCCCCGGGCTTGCCAGGAAACTAAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCAccc                                                                        <  2:522817/151‑1 (MQ=255)                                                                            cGTGGCGTCCACTCCCCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCAccc                                                                        <  2:414727/151‑1 (MQ=255)                                                                                    cccctccGCCCACCATAATTGACGAACGCTACTTGGTCCCCCCCGGGTTTGCCAGGAACTTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAat                                                                <  2:389988/145‑1 (MQ=255)                                                                                    ccACTCCGCCCCCCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAat                                                                <  2:1033/151‑1 (MQ=255)                                                                                      cctccccccACCATAATTGACGACCGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAatca                                                              <  2:446474/145‑1 (MQ=255)                                                                                      cctccGCCCACCAAAATTGACGACCGCTCACTGTTCCCCCCCGGGCTTGCCAGGAAACTTACCCCCGGTTTGGTCCTTTTGACTACCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAatca                                                              <  1:133627/147‑1 (MQ=255)                                                                                         tCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCctc                                                                      >  1:551606/1‑140 (MQ=255)                                                                                         tCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGAAAACTTAACCCCGGTTTGTTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCctc                                                                      <  2:551614/140‑1 (MQ=255)                                                                                                ccaTAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAAc                                                    <  2:83347/151‑1 (MQ=255)                                                                                                  aTAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCa                                                  <  2:399535/151‑1 (MQ=255)                                                                                                              cGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACTTCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGcc                                      <  2:320336/151‑1 (MQ=255)                                                                                                              cGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGc                                       >  2:331545/1‑150 (MQ=255)                                                                                                                   aCTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAggg                                 <  2:201689/151‑1 (MQ=255)                                                                                                                          cccccgGGGCTTGCCAGGAAACTTAACCCCGGTTTGGTCCGTTTGCCTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGc                          <  2:149263/151‑1 (MQ=255)                                                                                                                          ccccccGGGTTTGCCAGGAAACTTACCCCCGGTTTGTTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGc                          <  1:97341/151‑1 (MQ=255)                                                                                                                          ccccccGGGCTTGCCAGGAAACTTAACCCCGTTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGc                          <  2:395078/151‑1 (MQ=255)                                                                                                                                 ggtttGCCAGAAACTTTACCGCCGGTTTGGTCCTTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGgtgt                    <  1:565777/147‑1 (MQ=255)                                                                                                                                   cTTGCCGGGAAACTTAACCCCGGTTGGTTCCGTTTGACTATCAGCCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCg                  <  1:403265/150‑1 (MQ=255)                                                                                                                                      gCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTAcg              >  1:446080/1‑151 (MQ=255)                                                                                                                                                  aaCACCGGTTTGGTCCTTTTGACTAGCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGaa  <  2:379932/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TGCCTGTTCCCGGCTTACCACATGAACGCGAATATCTTCGTTTTCATCAGCCAGACCGTGAATACCGCTTGCGGTCGTGGCGTCCACTTCGCCCACCATAATTGACGAACGCTCACTGGTGCCCCCCGGGCTTGCCAGGAAACTTAACACCGGTTTGGTCCGTTTGACTATCAGTCCCGCCTCTTCAATCGCTTCGCGACGGGCAACATCTTCCACACTTTCACCCTCTTCAATCATCCCGGCAACCATCTCCAGTAGCCAAGGGGTTTCGCTGGTGTCGTACGCGGCAATCCGAA  >  NC_000913/3177383‑3177678

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: ‑