Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	62.1%	19.9 / 31.2	30	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (7/4);  new base C (0/18);  total (7/23)
Fisher's exact test for biased strand distribution p-value = 2.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-05

AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984293‑2984585                                                                                                                                                      |                                                                                                                                                 aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCacac                                                                                                                                                 <  2:301737/150‑1 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCacac                                                                                                                                                 <  2:21876/150‑1 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACAcc                                                                                                                                                <  2:489727/151‑1 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACAcc                                                                                                                                                <  1:454585/151‑1 (MQ=255)       tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGcc                                                                                                                                          <  2:576629/151‑1 (MQ=255)        aaTGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGcct                                                                                                                                         <  1:662969/150‑2 (MQ=255)              aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACTCCCTGGCCCccgc                                                                                                                                     <  2:142511/149‑1 (MQ=255)                       acaaGGGACCACAAAAAAGGTATGTCGATTTCCCACCCATCATTATTCAACAGATTAATGCCGGCACTTATTTATCGTTAAGGAGAAAATGCAGATCTTTTCGTCACGTTAAAT‑GTTCACCCCCACCCCCGGCCCCCCGCCGATACAAAGc                                                                                                                          <  1:716777/148‑1 (MQ=255)                               gCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                  >  1:526959/1‑151 (MQ=255)                                cccAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATACCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                 >  2:677119/1‑151 (MQ=255)                                  aaaaaaaTGTATGTCGACTTGCCAGCCATCATTATTAAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAAT‑GGTAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                               <  2:403955/148‑1 (MQ=255)                                  aaaaaaaTGTATGCCGTCTTGCCACCAATCATTATTCAACAGATTAAGGACGGAAAATATTTATCGGAAAGGAGAAAATGCAGAGCTATTCGTCACTTTAAAT‑GGTAAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                               <  1:268528/148‑1 (MQ=255)                                                     tGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                             >  2:7823/1‑150 (MQ=255)                                                           ccaTCATTATACAAAAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTGGTCACTTTCAAT‑GGTCATTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                         <  1:617023/148‑1 (MQ=255)                                                           ccaACATTATTCAAAAGATTAAGGGCGGCAAATTTTTATCGGAAGGGAGAAAAGGCAGAGCTATTCGTCACGTTCATT‑GTTAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                       <  1:120892/150‑1 (MQ=255)                                                              tcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCa                                                                                   >  1:6153/1‑151 (MQ=255)                                                                   attAAACAGATTAATGACGGCAGATATTTACCGGTAAGGAGGAAGTGCAGAGCTTTTCGTCGCTTTAAGG‑GGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                              <  1:228570/151‑1 (MQ=255)                                                                                                    gTAAGGAGGAAATCCAGAGCTATTCGTCACTTTCAAT‑GGTCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                              <  1:152678/150‑1 (MQ=255)                                                                                                    gTAAGGAGAAAATGCAGGCTTTTTCGTCACTTTCAAT‑GGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                              <  1:717927/150‑1 (MQ=255)                                                                                                     aaaGCAGGAAATGCAGCGCTTGCCGTC‑TATTCAAGGGGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                            <  1:389999/150‑1 (MQ=255)                                                                                                             aaaggaagacttATTCGTCACTTTCAAT‑GGTCAAGGCCGCCCCCTGGCCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                      <  1:539941/138‑1 (MQ=255)                                                                                                                  cAGAGCTTTTCGTCACGTTAAAG‑GGTAAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                <  1:644282/150‑1 (MQ=255)                                                                                                                         tttCGTCACGTTCAAT‑GGTCAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                        <  2:535710/150‑1 (MQ=255)                                                                                                                         tttCGTCACGTTCAAT‑GGTCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                        <  2:526971/150‑1 (MQ=255)                                                                                                                                        t‑GGTCAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac         <  2:39145/151‑1 (MQ=255)                                                                                                                                         tgttCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac         <  1:48471/148‑1 (MQ=255)                                                                                                                                          ggTCAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac         <  1:325478/150‑1 (MQ=255)                                                                                                                                          ggTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgaca        >  1:295202/1‑151 (MQ=255)                                                                                                                                            taaaTGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat       <  2:612760/148‑1 (MQ=255)                                                                                                                                            tCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat       <  2:706273/150‑1 (MQ=255)                                                                                                                                               aTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattga    >  1:206572/1‑150 (MQ=255)                                                                                                                                                 accACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGCTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc  >  1:167387/2‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAAT‑GGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACC  >  NC_000913/2984293‑2984585

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