Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	44.2%	70.2 / 13.9	43	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/2);  new base C (0/19);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 2.62e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCC  >  NC_000913/2984291‑2984586                                                                                                                                                       |                                                                                                                                                   cGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCcaca                                                                                                                                                    <  2:88574/151‑1 (MQ=255)    atatGTCATGAGAATGGAATAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                                                                 <  1:239279/151‑1 (MQ=255)              agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccgc                                                                                                                                      <  2:557976/151‑1 (MQ=255)                     aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATAc                                                                                                                                >  1:478125/1‑150 (MQ=255)                       aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACaa                                                                                                                              >  2:243411/1‑150 (MQ=255)                           gAGGGAGCCAAAAAAAGGTATGTCGACTTGCCAGCCATCATTATTAAACAGATAAAGGAGGGAAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATCCCACCCCCGGC‑CCCCCGCCGATACAAAGCGt                                                                                                                         <  1:618591/151‑1 (MQ=255)                             gggAGCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGtt                                                                                                                        <  1:310962/150‑1 (MQ=255)                                     aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                               >  1:16715/1‑151 (MQ=255)                                       aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                             >  2:335757/1‑151 (MQ=255)                                            gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                        >  1:361635/1‑151 (MQ=255)                                                     cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                               >  2:280108/1‑151 (MQ=255)                                                          cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                          >  2:67419/1‑151 (MQ=255)                                                            gccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                        >  1:530465/1‑151 (MQ=255)                                                              catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                      >  2:378525/1‑151 (MQ=255)                                                                   ttattCAACAAATTAATGACGGCAGATATTTTTGGGAAAGGAGGAAAGGCAGAGTTTTTCGCCACTTTAAAGGGTAAATCCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                      <  1:631741/146‑1 (MQ=255)                                                                    tattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                 >  1:743452/1‑150 (MQ=255)                                                                         aaCAGATAAATGACGGCAGATATTTATCGGTAAGGAGGAAATGAAGAGTTTTTCGTCACTTTCAATGGTAAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                               <  2:168402/147‑1 (MQ=255)                                                                               ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                     >  1:213709/1‑151 (MQ=255)                                                                                tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                     >  2:96055/1‑150 (MQ=255)                                                                                 aaTGACGGCAGTTATTTAGCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                     <  1:56941/149‑1 (MQ=255)                                                                                   tGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                  >  2:185294/1‑150 (MQ=255)                                                                                       ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                             >  2:613867/1‑151 (MQ=255)                                                                                          aGATATTTATGGGTAAGGAGAAAATGCAGGGTTTTTCGTCACTTTAAAGGGTAAATGCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                             <  2:203801/148‑1 (MQ=255)                                                                                          aGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                           >  1:126106/1‑150 (MQ=255)                                                                                             tatTTATCGGAAAGGAGAAAATGAAGAGTTTTTCGTCACGTTAATGGGTAAATGCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                          <  2:515370/148‑1 (MQ=255)                                                                                              atTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                      >  1:103533/1‑151 (MQ=255)                                                                                               tttATCGGTAAGGAGGAAAGGCAGAGTTATTCGTCCCGTTAAATGGCCAATGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                      <  2:241850/150‑1 (MQ=255)                                                                                                  agcGGTAAGGAGAAAATGGAGGTTTTTTCTTCACTTTCAATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                   <  1:117758/148‑1 (MQ=255)                                                                                                     ggTAGGAAGGAAATGCAGAGATATTCTTCACGTTAAAGGTTCAATCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                <  1:333024/150‑1 (MQ=255)                                                                                                     ggTAAGGAGAAAATGCAGGGCTTTTCTTCACGTAAAATGGTAAACGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                               <  1:697709/151‑1 (MQ=255)                                                                                                        aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                             >  1:395480/1‑150 (MQ=255)                                                                                                          ggaggaAATGCAGACCTTTTCGCACCGTTCAATGGTCAA‑GCCCCCCCCGGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                          <  2:298548/151‑1 (MQ=255)                                                                                                                 aTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                   >  1:401327/1‑151 (MQ=255)                                                                                                                    cagattttttCGCCACTTTAAATGTTCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                 <  2:105973/142‑1 (MQ=255)                                                                                                                      gacttATTCCTCCCTTTCAATGGTCACTCCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                              <  2:301106/147‑1 (MQ=255)                                                                                                                       agCTTTTCGTCACGTTAAAGGGTCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                             <  2:405406/151‑1 (MQ=255)                                                                                                                          tATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                           >  2:433957/1‑150 (MQ=255)                                                                                                                            ttCGTCACTTACAATGGTAAATGCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGTGTGACCGAGTGCGa                          <  1:235792/149‑1 (MQ=255)                                                                                                                                 cACGTTCAATGGCAAACCCCACCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                    <  2:633626/150‑1 (MQ=255)                                                                                                                                     tacagggttaattaccACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcca                  <  1:548720/134‑1 (MQ=255)                                                                                                                                        aaTGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg            >  2:199422/1‑151 (MQ=255)                                                                                                                                              cAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg      >  1:29777/1‑151 (MQ=255)                                                                                                                                               aaTGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATAGGGTGACCGAGTGCGATCGCGCCACcattgacattg      >  2:62851/1‑150 (MQ=255)                                                                                                                                                   ccACCCCCGGC‑CCCCCGCCGATGCAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc  <  1:147321/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CGAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCC  >  NC_000913/2984291‑2984586

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