Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,601,170	0	G	C	51.4%	42.6 / 68.8	35	A306G (GCA→GGA)	ftsX	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (14/3);  new base C (0/18);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 2.93e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-08

CGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACG  >  NC_000913/3601020‑3601319                                                                                                                                                          |                                                                                                                                                      cGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑t‑gtg                                                                                                                                                       >  1:953102/1‑151 (MQ=255)    ttttATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGcc                                                                                                                                                     >  1:794642/1‑150 (MQ=255)       tattatTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACtt                                                                                                                                                 >  2:902560/1‑151 (MQ=255)        attattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTc                                                                                                                                                <  1:728753/151‑1 (MQ=255)          tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTcc                                                                                                                                               <  2:846672/150‑1 (MQ=255)          tattCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCg                                                                                                                                              >  1:740050/1‑151 (MQ=255)               aGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAAcc                                                                                                                                         <  2:380003/151‑1 (MQ=255)                      aaGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATg                                                                                                                                  >  2:471451/1‑151 (MQ=255)                            cGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATc                                                                                                                            >  1:630118/1‑151 (MQ=255)                                                 aGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGAc                                                                                                        >  2:638933/1‑150 (MQ=255)                                                    cccgcccccccccAGTCAATCATCGACAACCCCAGCAGAAATAGAGGCCATCCAGAAAATAATACGCCTTTGAT‑CCAAACCTCCCCCCCAAAA‑CCTT‑GTCCCTCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACaa                                                                                                      <  1:346108/139‑1 (MQ=255)                                                     aCGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAAt                                                                                                   >  2:218915/1‑151 (MQ=255)                                                        ctccccccaacccAATCACCGACCATCCCAACAGAAATAGCAGCCATTCACCAAATAAAAAGCCATTGATA‑AAAACCTCCTTCCCAAAACCC‑T‑GTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTa                                                                                                 <  1:870487/139‑1 (MQ=255)                                                         tGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt                                                                                               >  2:145179/1‑151 (MQ=255)                                                          cccccccACCAAATCATCGAACAAACCAGCAGAAATAGCAGCCATTCATAGAATGATAAGCCTTTGAAAAAAATCTTCG‑TCCAAAAA‑CCGT‑GCCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAAt                                                                                               <  2:823562/146‑1 (MQ=255)                                                           ccccccacccaATCATCGAGCATACCAGCACCATTAGCAGGCTTTCATCAAATAATAAGCATTTAAAAACAACCTCCTTCCCAAAAA‑CC‑T‑GTCCCCCTCCCGAATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGa                                                                                             <  2:841340/148‑1 (MQ=255)                                                           ccccccACCAAATAATCGAGCATACCAGAAGAAATAGAAGCCATTAAAAGAATAAAAAGCAATAAATACAAATCTTCGTCCCAAAAC‑CC‑T‑GTCCCTCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGa                                                                                             <  2:583323/148‑1 (MQ=255)                                                                    cAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATgc                                                                                    >  2:601141/1‑151 (MQ=255)                                                                         accGACCATACCAGAAAAAATAGCAGCCATTAATAAAATAAAAAGCATTTAAAAAAAAACTTCGTCCCAAAAC‑CC‑T‑GTCCCCCTCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGacaaca                                                                                        <  1:708244/140‑1 (MQ=255)                                                                         accGACCATACCACCAAAAAAAGCAGGCTTTCATCAAATAATAACCTTTTAATATCAATCTTCTTTCCAAAAC‑CC‑T‑GTCCCCCTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCa                                                                                <  2:857133/148‑1 (MQ=255)                                                                             aGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGaaa                                                                            >  1:461096/1‑150 (MQ=255)                                                                                                     aTTAATCAAAGGATATCCCATTGAAATAAATCTTCTCCCCGAAAC‑CC‑G‑GCCCCCCTTCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATa                                                      <  1:594565/148‑1 (MQ=255)                                                                                                         aTCAAATGATAAGCCTTTGATAAAAACCTTCGTCCAAAAAC‑CC‑T‑GTCCCCCTCCCGAATCCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCAt                                                       <  2:687736/143‑1 (MQ=255)                                                                                                            gAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACgg                                             >  1:798/1‑150 (MQ=255)                                                                                                             aataatAAGCTATTGATATAAAACTTCTTCCCAAAAA‑CC‑T‑GTCCCACTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACg                                              <  1:114521/148‑1 (MQ=255)                                                                                                             aaTGAAAAGCCATTAATACAAAACTCCGAACCGAAAC‑CC‑T‑GACCCACCTCCAAACCCGCCGATGACAATCGCGGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACggg                                            <  1:930096/150‑1 (MQ=255)                                                                                                               tGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgcgc                                         >  2:608238/1‑151 (MQ=255)                                                                                                                 aTAAGCAATGGATACAAAACTCCCTCCAAAAAC‑CC‑T‑GTCCCACCCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGgc                                           <  1:751523/147‑1 (MQ=255)                                                                                                                    agcatttGATAAAAACCTTC‑TTCCCAAAA‑CC‑TGGTCCCCCCCCCACAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGAt                                    <  2:738100/146‑1 (MQ=255)                                                                                                                            ataaCAACCTTCGTCCAAAAAC‑CC‑T‑GCCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCAt                            <  1:780654/151‑1 (MQ=255)                                                                                                                                 aaaCTTCTTCCCAAAAA‑CC‑T‑GTCCCACCCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGt                        <  1:543374/150‑1 (MQ=255)                                                                                                                                   ccttcttccaaaaac‑CC‑T‑GTCCCCCTCCCCCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGc                     <  2:588246/140‑1 (MQ=255)                                                                                                                                             gAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGt           >  2:560677/1‑151 (MQ=255)                                                                                                                                               aaa‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt         >  2:254182/1‑151 (MQ=255)                                                                                                                                                acccc‑T‑GTCCCACTCCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt         <  1:572768/148‑1 (MQ=255)                                                                                                                                                aaccc‑T‑GTCCCCCTCCCCAAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGttt         <  2:707698/148‑1 (MQ=255)                                                                                                                                                      t‑gtgCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGt     >  2:169009/1‑150 (MQ=255)                                                                                                                                                         tgCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACg  >  1:646847/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                      CGCTTTTATTATTCAGGCGTAAAGTGGCGTAAATGTTGTACCGTGGCAAGCCACGCTGCCACCCAGCCAATCATCGAGCATACCAGCAGCAATAGCAGGCATTCATCGAATGATAAGCCATTGATATCAAACTTCGTTCCGAAAA‑CC‑T‑GTGCCACTTCCGCAACCGCCGATGACAATCGCAGCACCAGAATTTCTGACAAAATTAATGACAACAATGCGCCAGAAAATCCCAGCAGTGCGCCACCATACAGGAACGGGCGCAGGATGAATCCATCTGTCGCACCAATCAGTTTCTGTACG  >  NC_000913/3601020‑3601319

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