Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,210	0	C	A	40.5%	83.0 / 22.1	42	V13F (GTT→TTT)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (5/20);  new base A (17/0);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTT  >  NC_000913/1312061‑1312352                                                                                                                                                      |                                                                                                                                                  cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAc                                                                                                                                                   <  1:107920/150‑1 (MQ=255) cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAc                                                                                                                                                   <  2:398133/150‑1 (MQ=255) cATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAc                                                                                                                                                   <  1:620095/150‑1 (MQ=255)    tCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGCGAAAGAACGTAaaaaaa                                                                                                                                                 >  1:776136/1‑149 (MQ=255)     ccTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGa                                                                                                                                               <  1:104301/150‑1 (MQ=255)              cGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAg                                                                                                                                    <  2:86157/151‑1 (MQ=255)              cGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAg                                                                                                                                    <  1:753895/151‑1 (MQ=255)               gCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCGAAAAGACGTAAAAAAAAATcac                                                                                                                                           >  2:566649/1‑142 (MQ=255)                 acaGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGCTAAAAAACCCAAAAAAAAATCGGCCC‑Ggggg                                                                                                                                  >  2:431014/1‑150 (MQ=255)                        cTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCTCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGCTGGCTGAAGGAAGTAAAAAAAGACCGgcccgtagggggggatg                                                                                                                           >  2:146331/1‑137 (MQ=255)                           cgcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTAtgtg                                                                                                                        <  1:294439/150‑1 (MQ=255)                            gcgACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGtgtt                                                                                                                      <  1:904547/151‑1 (MQ=255)                               aCCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTa                                                                                                                   <  1:769058/151‑1 (MQ=255)                                 cGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGCAGGCCAAAATAAGGAAAAAAAAAACGGCCT‑GAGGGGGGGTATTtgttgt                                                                                                                    >  2:479670/1‑148 (MQ=255)                                    gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACAt                                                                                                              <  1:686719/151‑1 (MQ=255)                                         aTTTTTTTTTCCGGAAGAGGGCAGCCACTATCCATTTGTGGCATTACCCAACCTCCAAAGATGTCTCCATTAGCTGTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAATATCGCCCT‑GAGGGACGTTAGGTGTTGCAGACATGGTaa                                                                                                         <  2:478361/151‑1 (MQ=255)                                          tttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGGGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATAGTaa                                                                                                         <  1:708244/150‑1 (MQ=255)                                           ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCCCGGCCAAAGAACCTAAAACAAAAACGGCCT‑GCGGGAGGGTATGTGTTGTGGAAATGGGaa                                                                                                         >  2:296456/1‑149 (MQ=255)                                           ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCGAAAAAACGGAAAAAAAAAAAGCCCC‑GAGGGGCGGTATGTGTTGTAGAAATGGGaaa                                                                                                        >  2:418572/1‑150 (MQ=255)                                                  gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAAAAAGTAAAAAAAAAACGCCCC‑GAGGGGCGGTATGTGTTGTAGACAAGGGAAAAACG‑ACtt                                                                                               >  1:68507/1‑151 (MQ=255)                                                   ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTa                                                                                              <  2:136342/151‑1 (MQ=255)                                                      gAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCAAGGGTAAAAAAAGTAAAAAAAAAACGCCCC‑GAGGGGGGGTTTGTGTTGTGGAGAAGGGAAAAACGAACTTaa                                                                                             >  2:341659/1‑150 (MQ=255)                                                            cgccgcCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGaaaa                                                                                    <  2:422444/151‑1 (MQ=255)                                                              ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCCAAAAGAAGTAAAAAAAAAACGCCCC‑CTGGGGGGGTTTGTGTTGTTGAAAAGGGAAAAACG‑ACCTAAAA‑AGaaaa                                                                                    >  1:440071/1‑149 (MQ=255)                                                               cgccAATGTCCATTTGTGACATTAACCAGCCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATc                                                                                  <  1:737249/150‑1 (MQ=255)                                                               cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGGATGGCCAAAATAAGTAAAAAAAAATCGGCCTCGGGGGAGGTT‑TGTGTTGGAGAAAAGGGAAAAACG‑ACTTAAAA‑GGaaaa                                                                                    >  1:835229/1‑148 (MQ=255)                                                                        ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAGCAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGcgcg                                                                         <  1:701019/150‑1 (MQ=255)                                                                          aTTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATg                                                                      <  1:461658/151‑1 (MQ=255)                                                                          aTTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATAGCTAAAAAAAGTAAAAAAAAAACGCCCT‑TGGGGGGGGGTTTTGTTGTTGAAAAGGTAAAAACC‑ACCTAAAA‑AGaaaaaaacaagggggg                                                                        >  2:458906/1‑137 (MQ=255)                                                                             tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCGAGGCGAAAAGAAGTAAAAAAAAAATACCCC‑GAGGGGCGTTATGTGTTTTTGAAATGGTAAAACAC‑CCGTAAAACGGAAAAT‑ACAGgcg                                                                           >  1:526299/1‑143 (MQ=255)                                                                                  cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATg                                                              <  1:677758/151‑1 (MQ=255)                                                                                        ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGGTAAAATACGTAAAAAAAAAATGCCCT‑GAGGGGGGGTATGTGTTGTTGACAAGGGAAAACCG‑ACTTAAAA‑AGAAAAACACAAGGGGGGTGGTACTATGAttttt                                                        >  1:408853/1‑150 (MQ=255)                                                                                            ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTga                                                    <  1:736442/151‑1 (MQ=255)                                                                                            ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGAATGGGTAAAAGACGTAAAAAAAAAAACCCCT‑GGAGGGGGTGATGTGGTGTTGAAATGGGAAAAACC‑ACTTAAAA‑AGAAAATTACAAAAGCGGTGCTACGATAAtttt                                                         >  2:451666/1‑146 (MQ=255)                                                                                             caccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAAAAAAATCGCCCC‑GAGGGGCGTTATGTGTTGTAGACAAGGGAAAACCG‑ACTTAAAA‑AGAAAAAAACAGGCGCGGTGCTACTATGAtttt                                                         >  1:611211/1‑145 (MQ=255)                                                                                              accaAAGATGTCACCATTGGCATTGGTATCGGCGGGCCCGGCCAAAAAAAGTAAAAAAAAAAAGCCCC‑GAGGGGCGGTATGTGTTTTAGAAAAGGGAAAAACC‑ACTTAAAA‑AAAAAAACACAAGGGCGGTGCTCCTATTAtttt                                                         >  2:878893/1‑144 (MQ=255)                                                                                                caAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTgaggag                                                <  2:383419/151‑1 (MQ=255)                                                                                                 aaaGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGgaga                                               <  2:528932/151‑1 (MQ=255)                                                                                                   aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTGAAAGAAGTAAAAAAAGAACGCCCT‑GAGGGAGGGTGTGTGTTGTAGAAAAAGGAAAAACG‑AACTAAAA‑AGAAAAAAACAGGGGCGGTGCTACTATGATTTTTCTtagggg                                                >  1:242659/1‑146 (MQ=255)                                                                                                   aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAAc                                             <  1:418589/151‑1 (MQ=255)                                                                                                   aGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAATACGTAAAAAAAGAACGCCCC‑GGGGGGCGGTATGTGTTTTAGAAAAGGGAAAAACG‑ACTTAAAA‑AGAAAATAACAGGCGCGGTGCTACTATTATTTTTCTtaggagga                                              >  1:810262/1‑148 (MQ=255)                                                                                                                          tCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGcccc                      >  1:61236/1‑151 (MQ=255)                                                                                                                            ggcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGa                    <  2:666059/151‑1 (MQ=255)                                                                                                                                    tGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTgg            >  2:273599/1‑151 (MQ=255)                                                                                                                                               cGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt  >  1:573457/1‑150 (MQ=255)                                                                                                                                                gTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGtt  >  1:496262/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CATTCCTTCAACCCGCACAGTCACTACGCGACCGTGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCT‑GAGGGACGTTATGTGTTGTAGACATGGTAAAACCG‑ACTTAAAA‑GGAAAATCACAGGCGCGATGCTACTATGATTTTCCTGAGGAGAACAGAGGAGATACAGCGCCAGCCCCGAAGGACTGGCAGTCTGGTT  >  NC_000913/1312061‑1312352

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