Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,213	0	G	A	46.3%	39.6 / 23.5	41	L12F (CTT→TTT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (3/19);  new base A (19/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 8.01e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCT  >  NC_000913/1312076‑1312363                                                                                                                                          |                                                                                                                                                         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)                                                                                                                                      aCAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGAttt      <  1:656298/151‑1 (MQ=255)                                                                                                                                       cAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTa     <  2:841018/151‑1 (MQ=255)                                                                                                                                          gATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCt  >  2:533055/1‑151 (MQ=255)                                                                                                                                          |                                                                                                                                                         CACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTTAAAAGGAAAATCACAGGCGCGATGC‑TACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTATCT  >  NC_000913/1312076‑1312363

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