Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,964	0	A	C	47.8%	60.4 / 18.9	46	G332G (GGT→GGG)	yaaJ	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (20/4);  new base C (0/22);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.94e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACG  >  NC_000913/6818‑7103                                                                                                                                                    |                                                                                                                                             ttCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCtt                                                                                                                                          >  2:740275/1‑151 (MQ=255)    aGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCa                                                                                                                                       <  1:80164/151‑1 (MQ=255)        gTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCgg                                                                                                                                   >  2:667510/1‑151 (MQ=255)          aaaaaGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGc                                                                                                                                  >  2:1007070/1‑150 (MQ=255)           aaaaGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCAt                                                                                                                                >  2:365482/1‑151 (MQ=255)           aaaaGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCAt                                                                                                                                <  2:16454/151‑1 (MQ=255)               gaagaGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGt                                                                                                                              >  1:604800/1‑149 (MQ=255)                   agaTTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTt‑g                                                                                                                       >  2:40365/1‑151 (MQ=255)                      ttGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGc                                                                                                                    >  2:91519/1‑151 (MQ=255)                                cATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGc                                                                                                          >  2:305678/1‑151 (MQ=255)                                                aaCGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATg                                                                                          >  2:1042977/1‑151 (MQ=255)                                                aaCGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATg                                                                                          >  2:895941/1‑151 (MQ=255)                                                   gATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGgcg                                                                                       >  1:918699/1‑151 (MQ=255)                                                        aGCTGAAGAGAAGCAAAATATCCAAAGGGGAAACAACCTCAGCCCCCCAGCAACCC‑ATTCCCCCCATCTTCCCTTACTGGGTAAGCTGAATCCCTCCCGGCGTCGTGTGGATG‑CAGCACTTACCCGCCAGTAATATCAGCATGgcgc                                                                                      <  2:16233/147‑1 (MQ=255)                                                          cggAAGGCAAAAAAAATAACCACAGGGGTAACAAACCCAGCCCCCCAGGACCCC‑ATTACCCCCCGCATCGCCTTCTGAATAAGCTAAATCCCTCCCGGCGGCATTTAGTTT‑GTGCCTTTACCCGCCAGTAATATCAGCATGGCGCTTGCCg                                                                                <  2:123530/149‑1 (MQ=255)                                                           tGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGt                                                                               >  2:616010/1‑151 (MQ=255)                                                                  aaaaAGAATACCCACAAGGGTAACAAACTCAGCCCCCCAGGAACCC‑ATTACCCCCCGCATCGCCTTCTGGATAAGCTAAATCCCTCCCAGCGGCATGTGGTTT‑GTGCCTTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGAt                                                                         <  2:43885/150‑1 (MQ=255)                                                                   aaCAGAATACCCACAAGGGTAACAACCTCACCCCCCCAGGACCCC‑ATTACCCCCCCCATCGCTTTCTGGATAAGCTAAATCCCTCCCAGCGGAATGTGGTTT‑GGCCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGa                                                                       <  2:993440/151‑1 (MQ=255)                                                                   aaCAAAAAACCCAAAAGGGTAACAAACCCAGCCCCCCAGGACCCC‑ATTAGCCCCCGCATCGCCTTTTGGATAAGCTAAATCCCTTCCAGCGGCATGTGGGTT‑GTGCCTTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGa                                                                       <  2:288054/151‑1 (MQ=255)                                                                        aaaaaCCACAAGGGAAACAAACTCAGCCCCCAAGGAACCC‑ATTAGCACCCGCACCGCCTTCTGGAAAAGCTAAATCCCTCCCAGCGGCATG‑AGGTTGGTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAgg                                                                  <  2:13102/148‑1 (MQ=255)                                                                         aTAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAggg                                                                 >  2:516406/1‑151 (MQ=255)                                                                                    ggTAACAAACTCAGCACCCCGGGACCCC‑ATGAGCCCCCCCATCGCCTTCGGGATAACCTAAATACCTCCCAGCGGCATGTAGGTT‑CTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATaa                                                        <  2:684440/149‑1 (MQ=255)                                                                                    ggTAACAAACCCAGCACCCCAGGACCCC‑ATTACCCCCCGCATCGCCTTGTGGATAAGCTAAATCCCTTCCAGCGGCATGTAGGTT‑GTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATaaaa                                                      <  2:86850/151‑1 (MQ=255)                                                                                     gTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATaaaa                                                      >  1:193498/1‑150 (MQ=255)                                                                                       aaaaacccaaGCACCCC‑GGAACCCAATAACCCCCCCCATCGCCTTCTGGAAAGGTTAAATCCTTTCCAGGGGTATGTAGGTT‑GGCCCTTAACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATg                                                   <  2:606669/142‑1 (MQ=255)                                                                                         caaacccggacccccAGCACCCC‑ATTAGCACCCGCATCCCTTTCTGAATAAGATGAATCCTTTCCGGCGGAATGGGGGTT‑GCCCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGc                                                  <  2:488866/139‑1 (MQ=255)                                                                                         cAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGcc                                                 >  1:1050208/1‑151 (MQ=255)                                                                                              tCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatc                                             >  2:1051363/1‑150 (MQ=255)                                                                                               cacccccccAGAACCCC‑ATTACCCCCCGCATCGCCTTTGGGATAAGATAAATCCCTCCCAGCGGCATGGAGGTG‑GTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAatc                                             <  1:831414/144‑1 (MQ=255)                                                                                                  ccccccggtacccc‑TTTAGCCCCGGCATCGCTTTCTGGATAAGCTAAACCCCTTCCAGGGGCATGTAGGTT‑GGGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCt                                         <  2:833730/139‑1 (MQ=255)                                                                                                  cccccaAGGAACCC‑ATTCCCCCCCGCATCCCCTTTTGGAAAAGCTAAATCCCTCCCAGGGCCTTGTGGGTT‑GTGCCTTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTg                                        <  1:1055618/149‑1 (MQ=255)                                                                                                   aCCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTgg                                       >  1:55181/1‑151 (MQ=255)                                                                                                   aCCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTgg                                       >  1:509079/1‑151 (MQ=255)                                                                                                     cccAGGAACCC‑ATTAGCCCCCACATCGCCTTCTGAATAAGCTGAATCCCTCCCAGCGGCATGTAGGTT‑GGGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGAc                                     <  1:71213/151‑1 (MQ=255)                                                                                                        agaacccc‑ATTAGCACCCGCATCCCCTTCTGGATAAGCTGAATACCTCCCAGCGGCATGTAGTTT‑GTGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACaa                                   <  2:838934/145‑1 (MQ=255)                                                                                                         ggAACCC‑ATTACCCCCCGCTCCGCCTTCTGGATAAGCTAAATCCCTCCCGGCGGCATGTGGTTT‑GTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAAt                                  <  1:967052/150‑1 (MQ=255)                                                                                                           aaCCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGCGTCGGTAAAAATGCCAATCATCTGGACCATccc                               >  1:519316/1‑151 (MQ=255)                                                                                                             ccc‑ATTACCCCCCGCCTCGCCTTTTGGATAAGCTAAATCCCTTCCAGCGGTATGTGGGTT‑GGGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCt                              <  2:430934/150‑1 (MQ=255)                                                                                                                 aTTAGCCCCCGCATCGCCTTCTGGATAAGCTGAATCCCTCCCAGCGGAATGTAGGTT‑GTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTgc                           <  2:779096/150‑1 (MQ=255)                                                                                                                  ttAGCCCCCGCCTCCCTTTTTGAATAAGTTAAATCCCTTCCAGGGGCATGAGGGTT‑GGGCTGTGACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTg                            <  2:550698/148‑1 (MQ=255)                                                                                                                      cccccGCTTCCCCTTCTGGATAAGCTGAATCCCTCCCAGCGGCATGGAGGTT‑GGCCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTg                       <  2:71570/147‑1 (MQ=255)                                                                                                                       ccccGCATCGCCTTCTGGATAAGCTAAATCCCTTCCGGGGGTAGGTGGGTT‑GGCCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTgcgc                         <  1:881240/145‑1 (MQ=255)                                                                                                                         ccGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCgg                   >  1:478901/1‑150 (MQ=255)                                                                                                                                    tgtGGATAAGCTAAATCCCTCCCAGCGGCATGTAGGTT‑GTGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGcc        <  1:932070/148‑1 (MQ=255)                                                                                                                                         aTAAGCTAAATCCCTCCCGGCGGCTTGTGGGTG‑GGGCCGTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGAc   <  1:28034/150‑1 (MQ=255)                                                                                                                                           aaGCTAAATCCCTCCCGGCGGTATGTAGGTT‑GGGCCTTTCCCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGAcg  <  2:604814/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             TTCAGGCGTAAAAAGAAGAGATTGTTTTCGGCATAAATGTAGTTGGCAACGATGGAGCTGAAGGCAAACAGAATAACCACAAGGGTAACAAACTCAGCACCCCAGGAACCC‑ATTAGCACCCGCATCGCCTTCTGGATAAGCTGAATACCTTCCAGCGGCATGTAGGTT‑GTGCCGTTACCCGCCAGTAATATCAGCATGGCGCTTGCCGTACAGATGACCAGGGTGTCGATAAAAATGCCAATCATCTGGACAATCCCTTGCGCTGCCGGATGCGGAGGCCAGGACG  >  NC_000913/6818‑7103

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