Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,289	0	A	G	40.9%	54.4 / 17.4	45	G171G (GGA→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base A (8/18);  new base G (18/0);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 1.57e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATCAA  >  NC_000913/1157141‑1157436                                                                                                                                                      |                                                                                                                                                     aactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCCTACCGATACGCTCACAATTCTTCGCGAAGAAAAAGTCACGTATTGCAATGTCTTACTACACTTTTTTACAGACTACACAGAAACCGC‑ATATAAATTACTCGATCTCGGAt                                                                                                                                                     <  2:66940/150‑1 (MQ=255)     gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGAtttt                                                                                                                                                  <  2:220472/149‑1 (MQ=255)           gCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTcc                                                                                                                                          <  1:644856/151‑1 (MQ=255)               gTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCtttt                                                                                                                                      <  2:605037/151‑1 (MQ=255)               gTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGAGTGCGGTGGCGTACTACGCTGTTTTACAGAGGACAGAGACACGGC‑GGGTAAATTACTGGATCTCGGATTTTCCATCTCCtttt                                                                                                                                      <  2:334449/151‑1 (MQ=255)                      tCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGGTTTTACATCTCCTTTTTCGGCAt                                                                                                                               >  1:432131/1‑151 (MQ=255)                           aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGa                                                                                                                          <  1:412225/151‑1 (MQ=255)                           aCCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAATACTGGATCTCGGGGTTTAAATCTCCTTTTTCGGCATtgtg                                                                                                                           >  1:529169/1‑150 (MQ=255)                                 gACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTcc                                                                                                                    <  2:616850/151‑1 (MQ=255)                                   cGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGt                                                                                                                  <  1:734530/151‑1 (MQ=255)                                       cGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGGTTTCGGGTTTTAAATCTCCTTTTTCGGCCTTTTGAACTTCCGTaa                                                                                                                >  1:739589/1‑149 (MQ=255)                                            cGATACGCTGGCAATTCTGCGCGAAGAAAAAGGGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTATAATAATCGAGGTCGGGGTTTAAAGGTGCTTTTCCGAGCTTGTGGCGTTCCGTaaggacg                                                                                                           >  2:526989/1‑144 (MQ=255)                                                aCGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTGAAATAATGGAACTCTGGTTTTTCATTTTCTCTTCCCGCCTTGTGACCTTCCCTaaatctgtggaa                                                                                                      >  2:505555/1‑140 (MQ=255)                                                aCGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCaa                                                                                                      <  1:540007/150‑1 (MQ=255)                                                   cTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTACATTACTGGGAATCGGATTTTAAATTTCCTTTTTCGCGATTGTGGTCTTCCCTAATACgggggaaaat                                                                                                   >  2:200954/1‑142 (MQ=255)                                                    tGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGGTCTCGGGTTTTAAATCTCCTTTTTCGGCATTGTGACCTTTCGTAATGCGGGGCaa                                                                                                      >  1:9922/1‑146 (MQ=255)                                                     ggCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATTGGTCTCGGGTTTTACATTTCCTTTTTCGGCATTTTTACCTTTCGTCATGCGGAGGAAATg                                                                                                   >  1:375890/1‑148 (MQ=255)                                                        aaTTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGAt                                                                                             <  1:323680/151‑1 (MQ=255)                                                               cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTGGGTCTCGGGTTTTAAATCTCCTTTTTCGGCATTGTGGCCCTCCGTAAAGCGGGGCAACCGCGCGGcgctg                                                                                         >  1:807365/1‑148 (MQ=255)                                                                gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGGAATTTAATGAAT‑TCGGGTTTTACAAATCCCTTTTCCGCCTTGGGTCCCTCCGTTATGCGGAGTAACAGCGCGAtgctgc                                                                                        >  2:613617/1‑148 (MQ=255)                                                                 cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGGTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGGTCTCGGGTTTTTCATCTCCTTTTTCGGCATTGTGAGCTTCCGTAATGCGGGGGAACTGCGCGGTGCTGcgcg                                                                                     >  1:724857/1‑150 (MQ=255)                                                                    agaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTACTGGGTCTCGGGTTTTTCATCTCCCTTTTCGGCAATGTGACCTTCCGTAATGCGGAGCAAAAGAGCGATGTTGcgcggttt                                                                                 >  2:816948/1‑147 (MQ=255)                                                                         aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTATTTGATCTCGGGTTTTATATCTCCTTTTTCGGCATTGTGACCTTCCGTAATTTGGGGCAAAAACGCGATTGTGCGCGTTTTGt                                                                               >  2:382774/1‑148 (MQ=255)                                                                              aCGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTgg                                                                       <  2:431661/151‑1 (MQ=255)                                                                                ggATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATTATGGATCTCGGGTTTTAAATCTCCTTTTTCGGCATTGTGGACTTCCGTAATGCGGGGCAACTGCGCGATGCTGCGCGGTTTGTCCCCCCgg                                                                       >  1:755337/1‑149 (MQ=255)                                                                                 gATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATc                                                                    <  1:348340/151‑1 (MQ=255)                                                                                  aTTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATTTCTGTGTTTATAACTTCTTTTTTGGGATTTTGACCCTTTTTAATTCGGAGCAAAAGCGCGGTGCTGCGCGTTTTGTccccc                                                                          >  2:296979/1‑144 (MQ=255)                                                                                  aTTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCGCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCg                                                                   <  2:521211/151‑1 (MQ=255)                                                                                          ggCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATACTTGATCTCGGGTTTTAAATCTCCCTTTCCCGCCTTGTGACCCTCCCTAATGCGGGGCAACAGCGCGATGCTGCGCGGTATGTcccccctggccgggtaccgg                                                           >  1:107458/1‑136 (MQ=255)                                                                                            cGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTAATGGGTTTCGGGTTTTAAAACTCCCTTTTCGGCATTGTGACCTTCCCTAATGCGGAGCAACCGCGCGATGCTGCGCGTTATTTCCCCCCGGGTCGGTTAC‑tgggg                                                        >  1:660823/1‑149 (MQ=255)                                                                                            cGTACTACACTGTTTTACAGAGGACAGAGAAACGGA‑GGGCACAATACTGGAACTCCGGATTTATAATTACTTTTTCGCGGTTGTGTTCTTCCCTCAAACGAGGGAACTGCGCGAAGATTCGCGGTATGTACCCCCCGGTCGGGTAc                                                              >  2:558462/1‑146 (MQ=255)                                                                                                   cacTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATTCTGGGTCTCGGGATTTAAATATCCTTTTTCCGCCTTTTGACCTTCCGGAATACGGGGGAAATTGGCGAGGCTGCGCGGTTTTTTCCCCCCGATCGGTTTCTTGGGGGaaatg                                                  >  1:18816/1‑149 (MQ=255)                                                                                                       gTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTAATGGGTTTCGGGATTTATATCTCCTTTTTCGGCATTGTGGCCTTCCCTAATGCGGGGCAACCGCGCGATGCTGCGGGCTATTTCCCCCTGGGTCGGTTTC‑TGGTGGGAAATGAc                                               >  1:827720/1‑149 (MQ=255)                                                                                                                agGACAGAGAATCGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGTCCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCGGGATCGGTTAC‑TGGTGGAAACTGACTCaccttaccg                                    <  2:576425/151‑2 (MQ=255)                                                                                                                  gACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTg                                  <  1:244677/151‑1 (MQ=255)                                                                                                                        gaAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCTCCTTACCTTGCGCCgg                            <  2:81056/151‑1 (MQ=255)                                                                                                                              gc‑ggGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGc                      >  1:432686/1‑151 (MQ=255)                                                                                                                                 gggTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCAt                    >  2:504433/1‑151 (MQ=255)                                                                                                                                    tAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGc                      >  2:820384/1‑146 (MQ=255)                                                                                                                                            tGGATCTCGGATTTTACATCTCCTTTTCCGGCCTTGCGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGaaa            >  2:359518/1‑148 (MQ=255)                                                                                                                                            tGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGaaa            >  1:42887/1‑148 (MQ=255)                                                                                                                                              gATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAag          >  1:613261/1‑148 (MQ=255)                                                                                                                                                  tcGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAAt     >  2:627634/1‑149 (MQ=255)                                                                                                                                                   cGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAAt     <  1:817167/148‑1 (MQ=255)                                                                                                                                                    ggATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATCaa  <  2:485254/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTACTGGATCTCGGATTTTACATCTCCTTTTCCGGCATTGTGACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCGCGTTATGTCCCCCTGGATCGGTTAC‑TGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCGAGGAAAAGAGAATCAA  >  NC_000913/1157141‑1157436

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