Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,200	0	G	A	42.9%	47.6 / 16.7	42	T16I (ACT→ATT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (6/18);  new base A (18/0);  total (24/18)
Fisher's exact test for biased strand distribution p-value = 5.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGG  >  NC_000913/1312051‑1312350                                                                                                                                                      |                                                                                                                                                         gTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGt                                                                                                                                                         <  2:735037/151‑1 (MQ=255)      aaaGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGt                                                                                                                                                     <  2:37454/150‑1 (MQ=255)          tCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTTAAGTACGGaaaa                                                                                                                                                 >  2:175899/1‑150 (MQ=255)               ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGCAAAAAAAGTAaaaaaaaa                                                                                                                                            >  1:206957/1‑148 (MQ=255)                cTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGAAAGTAaaaaaaaa                                                                                                                                            >  1:722873/1‑147 (MQ=255)                          cacaGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAggg                                                                                                                                <  2:431833/151‑1 (MQ=255)                             aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAATGCCCTGCGGGGCg                                                                                                                             >  2:453291/1‑151 (MQ=255)                                      gcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAACGTAAAAAAAAAACGCCCCGGGGGGGGGTATGtgtt                                                                                                                    >  1:38235/1‑151 (MQ=255)                                         aCCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                                                                                                 <  1:347183/151‑1 (MQ=255)                                                gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                                                                                          <  1:44660/151‑1 (MQ=255)                                                gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGACTGGCCAAAATAAGTAAAAAAAAATAGCCCTGAGGGGGGTTGTGTGTTTTAGAAATg                                                                                                           >  2:619237/1‑150 (MQ=255)                                                     ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTCAAAAACGGAAAAAAAAATAGCCCCGGGGGGGGGTCTTTGTTGTTGACATGGG‑aaa                                                                                                     >  2:186066/1‑150 (MQ=255)                                                            gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCAAAAAACGTAAAACAAAAACGCCCTTCGGGGGGGTATGTGTTTTTGACAAGGT‑aaaaaacacc                                                                                              >  2:526481/1‑144 (MQ=255)                                                             ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTa                                                                                            <  1:53513/151‑1 (MQ=255)                                                                gAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTaaa                                                                                          <  1:291180/150‑1 (MQ=255)                                                                    aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTCGCCAAAAGACGTAAAAAAAGATAGCCCTGCGGGGGGGTATGTGTTTTCGATAAGGTAAAAACGGACTTaaa                                                                                          >  1:249550/1‑147 (MQ=255)                                                                    aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGAAAAAAAAAAACGCCCTGAGGGGGGTTGTGTGGTGTAGACATGGT‑AAAAAACACTTAAAAAGa                                                                                      >  2:637534/1‑150 (MQ=255)                                                                    aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAGGGCTCAAAAACGTAAAAAAAAAACGCACCTCGGGAGGCTATGGGGTGTTGAGATGGG‑AAAAACAACTTaaaaa                                                                                        >  2:640318/1‑147 (MQ=255)                                                                       gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAAt                                                                                  <  1:760333/151‑1 (MQ=255)                                                                             aaTATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTAAAAGACGTAAAAAAAAAACGCCCCGCGGGGCGGTTTGTGTTGTAGAAAAGGG‑AAAAACGCAGTAAAAAGaaaa                                                                                   >  1:640436/1‑144 (MQ=255)                                                                              atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATAAGTAAAAAAAAAAAGCCCTGAGGGGCGGTATGTGTTGTAGAAAAGGG‑AAAAACGACTTAAAAAGAAAATcaca                                                                              >  1:622438/1‑148 (MQ=255)                                                                                   cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAAAAAAAACGCCCCGAGGGGGGGTTTGTGTTTTTGAAAAGGG‑AAAAACGACTTAAAAAGAAAAAcacaag                                                                            >  1:481398/1‑143 (MQ=255)                                                                                    aTTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATAAGGAAAAAAAAAACGCCCCGAGGGGGGGTATGTGTTGTAGAAATGGG‑AAAACCGACTTAAAAAGAAAAAcacaaggg                                                                          >  1:600983/1‑142 (MQ=255)                                                                                         tgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGACTAAAGTAAGAAAAAAAAAAAAAGCCTGCGGGGCGGTATGTGTTGTAGAAAAGGT‑AAAAACGACATTAAAAGAAAAACAACGGCGCGGTGg‑tac                                                                >  2:132121/1‑150 (MQ=255)                                                                                           aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTAt                                                             <  2:498078/151‑1 (MQ=255)                                                                                                   caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAAGACGGTAAAAAAAAACGCCCCCAGGGGGGTTTTGTGTTTTAGAAAAGGG‑AAAAACCAGTTAAAAAGAAAAATACAAGCGCGCTGCTTACT‑TGAtttt                                                       >  1:48464/1‑149 (MQ=255)                                                                                                   caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCGAAAGTACGGAAAAACAAAACGCCCTGCGGGGGGGGGTGTGTTGTAGCAAAGGG‑AAAAACAACCTAAAAGGAAAAAAACAGGGGCGAGGG‑TACTAAGAatttt                                                      >  2:425417/1‑149 (MQ=255)                                                                                                    aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAAAAAAATCGCCCCGAGGGACGGTTTTTGTTGTAGACAAGGG‑AAAACCGACTTAAAAAGAAAAAAACAGGGGCGATGC‑TACTATGAtttt                                                       >  1:724895/1‑148 (MQ=255)                                                                                                     accaccAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCTGGCTAAAAAACGGAAAAAAAAAATGCCCCGGGGGGCGGTATGTGTTTTAGACAAGGG‑AAAAACCACTTTAAAAGAAAAAAACAAGGGCGATGC‑TACTATGAtttt                                                       >  1:289185/1‑147 (MQ=255)                                                                                                      ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGCTAAAGAAAGCAAAACAAAAACGGCCTTAGGGGGGGTTTGTGTTGTTGACAAGGT‑AAAAACGACCTTAAAAGAAAAAAACAAGGGCGGGGC‑CTCTATGAtttt                                                       >  2:40514/1‑146 (MQ=255)                                                                                                        accaAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCTAAAAGAAATAAAACAAAAATGCCCCGAGGGGGGGGTTGTGTTGTTGAAAAAGT‑AAAAACAACTTTAAAAGAAAAAAACAAGCGCGGTGC‑CACAATGAtttttttttgg                                                >  2:605766/1‑144 (MQ=255)                                                                                                             aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAAc                                           <  1:611899/151‑1 (MQ=255)                                                                                                                  tCACCATTGGCATTGGTATCGGCGGGCGTGGCTAAAATACGTAAAAAAAAAACGCCCTTCGGGGCGGTATGTGTTGTTGACAAGGG‑AAAAACGACTTAAAAAGAAAAACACAGGCGCGGGGC‑TACTATGGTTTTTCTGAGGAGGAAAgagg                                      >  2:49738/1‑151 (MQ=255)                                                                                                                   cACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAgagg                                      <  1:798419/150‑1 (MQ=255)                                                                                                                   cACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAgagg                                      <  2:496416/150‑1 (MQ=255)                                                                                                                          ggcattggTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAg                              <  2:837659/151‑1 (MQ=255)                                                                                                                                   aTCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGccc                     <  1:434240/151‑1 (MQ=255)                                                                                                                                           gCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGAc             <  2:564999/151‑1 (MQ=255)                                                                                                                                             aTGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTg           <  2:820423/151‑1 (MQ=255)                                                                                                                                                   aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGt      <  1:790795/150‑1 (MQ=255)                                                                                                                                                   aaaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTc     <  1:662258/151‑1 (MQ=255)                                                                                                                                                     aGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTg   >  1:596094/1‑151 (MQ=255)                                                                                                                                                      gTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTgg  <  1:835335/151‑1 (MQ=255)                                                                                                                                                      gTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTgg  >  1:220638/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                         GTAAGAAAGTCATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGG  >  NC_000913/1312051‑1312350

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