Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,213	0	G	A	40.0%	49.0 / 28.7	45	L12F (CTT→TTT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (6/21);  new base A (18/0);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 7.92e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTAT  >  NC_000913/1312079‑1312361                                                                                                                                        |                                                                                                                                                       aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTAAGTAAAAAAAAAACGCCCTGCgggggg                                                                                                                                        >  1:783482/1‑148 (MQ=255)  gTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑AGGCTAAAAAACCGAAAACAAAATCGccccccggggcg                                                                                                                                        >  1:401881/1‑141 (MQ=255)    cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTa                                                                                                                                     <  1:686755/151‑1 (MQ=255)       tACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGA‑ATGCACAAAAACGGAAAACAAAAACGCCCttgggggagct                                                                                                                                      >  2:917198/1‑137 (MQ=255)        aCGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAtg                                                                                                                                   <  1:667118/149‑1 (MQ=255)           cgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGC‑AGGCGAAAGAAAGTAAAAAAAAATAGCCCCCAGGGGCGGTTTGtttt                                                                                                                               >  1:403866/1‑150 (MQ=255)           cgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttg                                                                                                                              <  2:402755/151‑1 (MQ=255)            gACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGA‑TTGCTAAAATACCGAAAACAAGAATGCCCCGCCGGGCGTTATGagttgt                                                                                                                             >  1:902588/1‑151 (MQ=255)              ccGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGA‑AGACTAAAAAAAAAAAAACAAAAAAGCCCTGAGGGGCGGTATGtgttgttg                                                                                                                           >  2:578024/1‑149 (MQ=255)               cGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTTAAGTACGTAAAACAAAAACGCCCTGAGGGGCGGTATGtgtt                                                                                                                               >  1:195028/1‑146 (MQ=255)                gTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGA‑TGGCGAAAAAAAGTAAAAAAAAAACGCCCCGAgggggggtat                                                                                                                                    >  1:627170/1‑133 (MQ=255)                  gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACa                                                                                                                        <  1:182744/150‑1 (MQ=255)                   ggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                                                                                                      <  1:390123/151‑1 (MQ=255)                        tttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGGTAAAAGAAGTAAAAAAAAAAAGCCCTGAGGGGCGGTATGTGTTGTAGaaaagggaa                                                                                                                  >  1:840584/1‑142 (MQ=255)                                 ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGC‑TGGGTAAAAAAAGTAAAAAAAAAAAGCCCTGAGGGGCGGTATGTGTTGTTGAAAAGGG‑AAAAACGACtt                                                                                                        >  1:629127/1‑150 (MQ=255)                                   aGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGGTAAAAGACGTAAAAAAAAAACGCCCCGAGGGGCGGTATTTGTTGTAGACATGGG‑AAAAACCACTT‑aa                                                                                                     >  1:59081/1‑150 (MQ=255)                                        aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGa                                                                                                <  1:869291/150‑1 (MQ=255)                                        aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCCAAAATAAGGAAAAAAAAAACGCCCCGAGGGGCGGTATGTGTTGTAGACAAGGG‑AAAAACGACCT‑taaa                                                                                                   >  1:832107/1‑147 (MQ=255)                                         gcgcCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGa                                                                                                <  2:104069/149‑1 (MQ=255)                                           gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCC‑TGGCTAAAAAACGTAAAAAAAAATCGCCCCGAGGGGCGTTATGTGTTGTAGACAAGGG‑AAAAACCACTT‑aaaaaaaaaa                                                                                             >  1:629407/1‑150 (MQ=255)                                                  atatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGgc                                                                                     <  1:236866/151‑1 (MQ=255)                                                    atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGgcgc                                                                                   <  2:432985/151‑1 (MQ=255)                                                     tCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAGGGGGTAAAAAAGGTAAAA‑AAAAACGCCCCCAGGGGCGCTGTGAGTTTTAGGCATGGTGAAAACCATATT‑AAAGGGaaaaaaaaggt                                                                                      >  1:237033/1‑141 (MQ=255)                                                         tttGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGACA‑AAGCTCAAAAAAAAAAAAAAAAAACGCCCCGGGGGGCGGTGTTTGTTGTAGAAAAGGG‑AAAAAAGACTTAAAAAGGaaaaaaaa                                                                                         >  2:57927/1‑137 (MQ=255)                                                          ttGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCC‑TGGCTAAAAAAAACAAAAAAAAATCCCCCTCAAGGACGGTATTTGTTGTAGAAAAGGG‑AAAAACGACCT‑AAAAAGAAAAACACAAGCGGGATGct                                                                             >  2:735069/1‑151 (MQ=255)                                                                cATTAACCAACCGCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTAt                                                                        <  1:31514/150‑1 (MQ=255)                                                                cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATg                                                                       <  2:437/151‑1 (MQ=255)                                                                 aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGGATGGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGTGACGTTATGTGTTGTAGACATGGT‑AAATCCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTAt                                                                        <  2:453204/149‑1 (MQ=255)                                                                   tAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGAt                                                                     <  1:353095/150‑1 (MQ=255)                                                                     accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGAtttt                                                                  <  2:30915/151‑1 (MQ=255)                                                                        aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGGAAGGAAAAAAAAAACGCCCCGAGGGGCGGTTTTTGTTGTAGACATGGT‑AAAAACGACCT‑TAAAAGAAAAAAACAAGGGGGGTGCCACTATGAtttt                                                                  >  1:573166/1‑148 (MQ=255)                                                                                 aGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAAAAAAAACGCCCCGAGGGGCGTTATGTGTTGTAGACAAGGT‑AAAAACGACTT‑TAAAAGAAAAACACAAGGGCGACGCTACTATTATTTTTCTtagggg                                                         >  1:431291/1‑146 (MQ=255)                                                                                          cattggcattggTATCGGCGGGCC‑TGGCGAAAGTACGTAAAAAAAAAAAGCCCTTCGGGGCGGTTTGTGTTGTAGACAAGGT‑AAAAAAGACCT‑CAAAAGAAAAATACAAACGCGACGCTACTCTGAATTTTCTGGGGGGAAAAgagg                                                 >  1:870331/1‑147 (MQ=255)                                                                                          cattggcattggTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGAt                                             <  2:680782/151‑1 (MQ=255)                                                                                            ttggcattggTATCGGCGGGCA‑TAGCTAAAAGATGAAAAACAAGATCGGCCGGCGGGGCGGTATGTGTTTTTGTGAAGGG‑AAAAACGACCT‑CCAAAGAGAAACGCAAGGCCGGTGGCACAATGATTTTTCTCCGGAGAGCAGAGGAGAAAc                                           >  2:260138/1‑151 (MQ=255)                                                                                                       aTCGGCGGGCA‑TGGCTAAAGTACGTAAACCAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGcc                                 <  2:470072/150‑1 (MQ=255)                                                                                                             gcgCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTGATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAg                           <  2:595097/148‑1 (MQ=255)                                                                                                                cA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACt                       <  1:167473/151‑1 (MQ=255)                                                                                                                         aaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCt               >  2:634322/1‑151 (MQ=255)                                                                                                                         aaGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCt               >  1:78999/1‑151 (MQ=255)                                                                                                                          aGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTg              <  1:404708/151‑1 (MQ=255)                                                                                                                           gTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTgg             <  2:626801/151‑1 (MQ=255)                                                                                                                               gTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAg         <  1:529935/151‑1 (MQ=255)                                                                                                                                     cAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCAAGCCCCGAAGGACTGGCAGTCTGGTTAGGAttt    >  2:140820/1‑150 (MQ=255)                                                                                                                                      aaGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTAt  >  1:167553/1‑151 (MQ=255)                                                                                                                                        |                                                                                                                                                       AGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCA‑TGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGT‑AAAACCGACTT‑AAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTTAGGATTTAT  >  NC_000913/1312079‑1312361

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