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.9%	50.5 / 29.3	44	L12F (CTT→TTT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base G (5/21);  new base A (18/0);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 3.40e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTT  >  NC_000913/1312064‑1312352                                                                                                                                                       |                                                                                                                                             tCCTTCAACCCGCACAGGCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAg                                                                                                                                              <  2:527571/150‑1 (MQ=255)     tCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATc                                                                                                                                           <  2:132582/149‑1 (MQ=255)                aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                                                                                              <  2:606431/151‑1 (MQ=255)                aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                                                                                              <  1:591731/151‑1 (MQ=255)                aGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGAGCGTTAAAAACAAAAGCCCTGAggg                                                                                                                                 >  2:768048/1‑148 (MQ=255)                  tCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CCTCGCTAAAGTAACTTAAAAAAAAACGCCCTGTGGGGCggt                                                                                                                            >  1:12470/1‑149 (MQ=255)                   cACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTa                                                                                                                           <  1:735947/151‑1 (MQ=255)                      tACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAATAAGGAAAAAAAAAACGCCCCCAGGGGCGGTAtgt                                                                                                                        >  1:375734/1‑151 (MQ=255)                       aCGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAtgt                                                                                                                        <  2:404816/150‑1 (MQ=255)                           gACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                                                                                                   <  2:135305/151‑1 (MQ=255)                            aCCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGGGTAAAAACGTAAAAAAAAAAAGCCCCCAGGGGCGGTATGtgtt                                                                                                                     >  1:771715/1‑148 (MQ=255)                              cGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGa                                                                                                                <  1:281234/151‑1 (MQ=255)                                    cAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑GCTGGGCAAAAAAACTTAAACAAAAATACCCTGAGGGGCGGTTTGggttgt                                                                                                                   >  2:132905/1‑142 (MQ=255)                                     aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCGCAAAAACGTAAAAAAAAATCGCCCCGTGGGGCGGGTTGTGTTGTTGAAAAGGt                                                                                                          >  2:507987/1‑150 (MQ=255)                                             ttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCg                                                                                                   <  1:553800/149‑1 (MQ=255)                                             ttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGa                                                                                                  <  1:3975/150‑1 (MQ=255)                                               gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑GATGGCTAAAAAAAGTAAAAAAAAAACGCCCTGAGTGAGGGTATTTTTTGGATGCATGGTAAAAACAA‑Ctt                                                                                              >  2:190498/1‑151 (MQ=255)                                                ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTa                                                                                             <  2:376373/151‑1 (MQ=255)                                                    aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAg                                                                                         <  2:862274/151‑1 (MQ=255)                                                    aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CAAAGCTAAAAAAAGTAAAAAAAAAACGCCCCGAGGGGAGGTTTGTGTTGTTTACCAGGGAAAAACAA‑CTTtaaa                                                                                          >  2:72511/1‑150 (MQ=255)                                                       aGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGaa                                                                                      <  2:187574/151‑1 (MQ=255)                                                          gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGTAAAAt                                                                                   <  2:848977/151‑1 (MQ=255)                                                          gccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CAGGGCGAAAACAAGGAAAAAAAAAAAACCCCGAGGGGCGGTTTGGTTTGTAGAGAAGGGGAAACCCA‑GTTAAAAAGcaaa                                                                                    >  1:828794/1‑150 (MQ=255)                                                             gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑GATGGCTAAAAAAAGTAAAAAAAAAACGCCCTGAGGGGGGGTATTTGTTTTTGAAAAGGGAAAAACGACCTTAAAAGAaaaaaa                                                                                  >  1:582912/1‑148 (MQ=255)                                                              ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATcaca                                                                               <  2:243573/151‑1 (MQ=255)                                                               cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑GACGGGTCAAAAAAATAAAAAAAAAAAACCCCGGGGGGCGGTATGTGGTGTTGGAAAGGGAAAAAAAA‑CTTAAAAAGaaaa                                                                                    >  2:511619/1‑145 (MQ=255)                                                               cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCCAAAGAACGTAAAAAAAAATCGCCCCGAGGGGCGGTATGTGTTTTTGACATGGGAAAAACCA‑CTTAAAAAGaaaaaaac                                                                                >  1:665040/1‑145 (MQ=255)                                                               cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CAAGGCTAAAAAAAGTAAAAAAAAAAAGCCCTGCGGGGCGGTATTTGTTGTTGAAAAGGGAAAAACCA‑CCTAAAAAGaaaa                                                                                    >  1:371076/1‑145 (MQ=255)                                                                aaTATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCATGTC‑AAAAAAAGGAAAACAAAATAGCCCCTCGGGGGGGTGTGTGTTGTAGATATGGGAAAAAAGA‑ATTAAAAAGAAAAAAACAgg                                                                             >  2:434425/1‑151 (MQ=255)                                                                     ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CAAAGCCACAAAACGTGAAACAAAAACGCCCCCAGGGACGGTTTGGGATGTAGAAAAGGGCAAAACGA‑AATAAAAAGaaaaaaaacagcg                                                                           >  2:792039/1‑139 (MQ=255)                                                                            tgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CCGGGGTAAAGTACGGCAAAAAAAACCCCCCTGTGGGACGGTGTATTTTGTGGAAATGGGGAAAACCA‑CTTAAAAAGAAAAacacccccgggaaggcgc                                                                  >  1:255183/1‑136 (MQ=255)                                                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTAt                                                               <  1:8850/151‑1 (MQ=255)                                                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTAt                                                               <  2:582773/151‑1 (MQ=255)                                                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTAt                                                               <  1:301150/151‑1 (MQ=255)                                                                                     ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTc                                                        <  2:262823/151‑1 (MQ=255)                                                                                          caccaAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTgag                                                   <  2:333143/151‑1 (MQ=255)                                                                                            ccaAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTgagga                                                 <  1:201463/151‑1 (MQ=255)                                                                                               aaGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGGTAAAAAAAGGAAAACAAAATCGCCCCGAGGGGCGGTATGTGTTGTAGAAAAGGTAAAAAAAA‑CCTAAAAAGAAAATAAAAAGCGCGATGCTACTATGATTTTTCtttgg                                                  >  2:156294/1‑143 (MQ=255)                                                                                                          attggcattggTATCGGCGGG‑GATGGGTAAAGAAAGTAAAAAAAAAACGCCCTGAGGGGAGGTATGTGGTTTTGACATGGGAAAAACCA‑AATAAAAAGAAAAAAAAAGGCGGGGGGCTACTATTATTTTTCTGAGGGGAAAAgagg                                        >  1:849752/1‑146 (MQ=255)                                                                                                             ggcattggTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACa                                 <  1:555602/150‑1 (MQ=255)                                                                                                                        cggcggG‑CAATGCTAAAGAACGGAAAACAAGATCGCCCTGAGGGGCGGTATGGGTTGTAGACATGTTAAAAACCA‑ATTAAAAGGAAAAACACAGGCGGGGTGGTACTATGGTTTTACTGAGGAGAACAGAGGAGATACCACGCCAGcccc                      >  2:509243/1‑150 (MQ=255)                                                                                                                                   ggCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGc           >  2:819005/1‑151 (MQ=255)                                                                                                                                     cTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATATAGCGCCAGCCCCGAAGGACTGGc           >  1:104109/1‑149 (MQ=255)                                                                                                                                     cTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGc           >  2:24718/1‑149 (MQ=255)                                                                                                                                              gTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt  >  2:302252/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                             TCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGG‑CATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGA‑CTTAAAAGGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTT  >  NC_000913/1312064‑1312352

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