Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,930,842	A→G	G84G (GGT→GGC)	ravA ←	regulatory ATPase RavA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,930,842	0	A	G	94.4%	106.8 / ‑4.2	36	G84G (GGT→GGC)	ravA	regulatory ATPase RavA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (18/16);  total (19/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCA  >  NC_000913/3930716‑3930977                                                                                                                               |                                                                                                                                        ggTATTAAGAATTGCGGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAActtctt                                                                                                                            <  1:301012/140‑1 (MQ=255) ggTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAActtctt                                                                                                                            <  2:61084/140‑1 (MQ=255) ggTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAActtctt                                                                                                                            >  1:527076/1‑140 (MQ=255)    aTTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCg                                                                                                                         <  2:624247/140‑1 (MQ=255)      tAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGc                                                                                                                       >  1:521184/1‑140 (MQ=255)      tAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTACCGCCTGACTAGGAAGGGGGCCAAAAACTTCTTCCGGc                                                                                                                       >  1:249008/1‑140 (MQ=255)      tAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCGTCATCTTTAAGCGCCTGAATAGAAAGGGGGGCAAAAACTTCTTCCGGc                                                                                                                       >  2:207655/1‑140 (MQ=255)            ttGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCAAAAAACTTCTTCCGGCGTGgag                                                                                                                 >  2:323289/1‑140 (MQ=255)                                     aTCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAAcg                                                                                        >  1:407519/1‑140 (MQ=255)                                     aTCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAAcg                                                                                        >  2:719214/1‑140 (MQ=255)                                     aTCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAAcg                                                                                        >  2:314853/1‑140 (MQ=255)                                     aTCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCAGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAAcg                                                                                        >  2:535031/1‑140 (MQ=255)                                            aaGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGc                                                                                 <  2:171011/140‑1 (MQ=255)                                                   tttCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGa                                                                          >  1:549326/1‑140 (MQ=255)                                                   tttCTGCTTCCGGCAAGTAACCACTCGTTTAACGTTCCTAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGa                                                                          <  1:207655/140‑1 (MQ=255)                                                         cTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCAt                                                                               <  2:109093/129‑1 (MQ=255)                                                         cTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCAt                                                                               >  1:109093/1‑129 (MQ=255)                                                                   gTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAgc                                                          >  1:677412/1‑140 (MQ=255)                                                                   gTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAgc                                                          >  1:53380/1‑140 (MQ=255)                                                                       ccGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCg                                                      <  2:677412/140‑1 (MQ=255)                                                                        cGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCgg                                                     <  2:549326/140‑1 (MQ=255)                                                                         gCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCggg                                                    <  2:6928/140‑1 (MQ=255)                                                                         gCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCggg                                                    <  2:527076/140‑1 (MQ=255)                                                                            ggTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATACTGAAAAGCGAATATTAAGCGCCGGGCGa                                                 >  1:359312/1‑140 (MQ=255)                                                                                               ccTTCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATAcc                              >  2:725715/1‑140 (MQ=255)                                                                                                 ttCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTg                            >  1:561430/1‑140 (MQ=255)                                                                                                  tCATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTgg                           >  2:524752/1‑140 (MQ=255)                                                                                                  tCATCTTTTAGCGCCAGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTgg                           <  1:314853/140‑1 (MQ=255)                                                                                                   cATCTTTTAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGc                          >  2:676344/1‑140 (MQ=255)                                                                                                         ttAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAGACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCa                    <  2:144890/140‑1 (MQ=255)                                                                                                         ttAGCGCCTGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCa                    >  1:715799/1‑140 (MQ=255)                                                                                                                 tGAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGaa            <  1:719214/140‑1 (MQ=255)                                                                                                                  gAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAac           <  2:16704/140‑1 (MQ=255)                                                                                                                  gAATAGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAac           <  1:46011/140‑1 (MQ=255)                                                                                                                      aGAAAGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAacacac       <  2:53380/140‑1 (MQ=255)                                                                                                                         aaGGGGGCCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACttt    <  2:155424/140‑1 (MQ=255)                                                                                                                           ggggGCCAAAAACTTCATCCGGAGTGGAGAAGCAGATCATCAGATATTAAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCAAAGGAGGAACACACTTTca  <  2:359312/140‑1 (MQ=255)                                                                                                                               |                                                                                                                                        GGTATTAAGAATTGCCGGGCCCGCTTTCCAGATCTCATCCAGAAAGACAATTTCTGCTTCCGGCAGGTAACCGCTGGTTAAACGTTCATAGCGCCCTTCATCTTTTAGCGCCTGAATAGAAAGGGGACCAAAAACTTCTTCCGGCGTGGAGAAGCGGGTCATCAGATATTCAAACGCGCGGGCATTCTGAAAGGCGAATTTTAAGCGCCGGGCGATCAAACTTTTGGCAATACCTGGCGGGCCAAGGAGGAACACACTTTCA  >  NC_000913/3930716‑3930977

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