Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,909	C→A	R149R (CGC→CGA)	cysZ →	sulfate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,909	0	C	A	91.7%	71.9 / 1.5	24	R149R (CGC→CGA)	cysZ	sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (11/11);  total (12/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

AACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGT  >  NC_000913/2531772‑2532045                                                                                                                                          |                                                                                                                                         aaCGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCgcgagc                                                                                                                                        <  1:29875/140‑1 (MQ=255)       ttATTGGCTGAACAACTGGAAGCACGATTGACTGGCGGTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGt                                                                                                                                  <  1:210856/140‑1 (MQ=255)            ggCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTAc                                                                                                                             <  2:335599/140‑1 (MQ=255)             gCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACt                                                                                                                            <  2:290525/140‑1 (MQ=255)                    aacTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGC‑‑‑‑CGCAATTGTATTACTAATTCTTTACt                                                                                                                 >  1:417040/1‑140 (MQ=255)                           aGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCa                                                                                                              >  1:38162/1‑140 (MQ=255)                               cGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAGAATTTGTCTGGTATCTGCCGCGAGCAATTGTATTTGCAATTCTTTCCTTCATccc                                                                                                          <  2:339134/140‑1 (MQ=255)                                                 ccGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGt                                                                                        <  2:38162/140‑1 (MQ=255)                                                      aGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCgcgc                                                                                   <  1:208548/140‑1 (MQ=255)                                                              ggATTTTCGATATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACtg                                                                           >  2:186163/1‑140 (MQ=255)                                                                          tCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGttt                                                               >  1:238526/1‑140 (MQ=255)                                                                                 aGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCt                                                        <  2:238526/140‑1 (MQ=255)                                                                                    tgtgCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGa                                                     >  1:424099/1‑140 (MQ=255)                                                                                    tgtgCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGa                                                     >  1:195097/1‑140 (MQ=255)                                                                                     gtgCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGAt                                                    <  2:195097/140‑1 (MQ=255)                                                                                      tgCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATg                                                   >  1:59414/1‑140 (MQ=255)                                                                                         cgcgAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTa                                                >  2:187382/1‑140 (MQ=255)                                                                                          gcgAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAg                                               >  1:461350/1‑140 (MQ=255)                                                                                              aTCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCAt                                           <  2:438863/140‑1 (MQ=255)                                                                                                        gcgAATGGCAAAAATTTGCCTGGTATCTGC‑‑‑‑CGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGAtt                             <  2:417040/140‑1 (MQ=255)                                                                                                           aaTGGCAAAAATTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGAt                              >  1:274574/1‑140 (MQ=255)                                                                                                                  aaaaTTTGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCt                       <  2:461350/140‑1 (MQ=255)                                                                                                                       ttGCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGAt                  <  2:59414/140‑1 (MQ=255)                                                                                                                         gCCTGGTATCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATaa                >  2:264642/1‑140 (MQ=255)                                                                                                                                 tCTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAg        >  1:19415/1‑140 (MQ=255)                                                                                                                                  cTGCCGCGAGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGt       <  1:186163/140‑1 (MQ=255)                                                                                                                                       gcgagcAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGt  <  1:187382/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         AACGGTTTATTGGCTGAACAACTGGAAGCACGATTGACTGGCGCTACACCGCCAGATACCGGGATTTTCGGTATCATGAAAGATGTGCCGCGAATCATGAAACGCGAATGGCAAAAATTTGCCTGGTATCTGCCGCGCGCAATTGTATTACTAATTCTTTACTTCATCCCTGGTATTGGGCAAACCGTCGCGCCGGTACTGTGGTTCCTGTTTAGCGCCTGGATGTTAGCCATCCAGTATTGCGATTACCCCTTCGATAACCACAAAGTGCCGT  >  NC_000913/2531772‑2532045

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