Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	494,736	A→G	H111R (CAT→CGT)	recR →	DNA repair protein RecR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	494,736	0	A	G	96.8%	100.9 / NA	31	H111R (CAT→CGT)	recR	DNA repair protein RecR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (19/11);  total (19/12)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGCTACC  >  NC_000913/494610‑494863                                                                                                                               |                                                                                                                                gTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGACAGTTTTCGGGTGGTTATTTTGTGTGGATGGGGCGTCTGTc                                                                                                                           >  2:134040/1‑133 (MQ=255)    aCATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCAcc                                                                                                                        >  2:81255/1‑133 (MQ=255)    aCATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCAcc                                                                                                                        >  1:322622/1‑133 (MQ=255)    aCATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCATTATTTTTTGTTGATGGGGCGACTGTCAcc                                                                                                                        >  1:174699/1‑133 (MQ=255)    aCATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAACCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTCTGTGGAAGGGGCGCCTGTCAcc                                                                                                                        >  1:186829/1‑133 (MQ=255)     cATCTGTTCGTATCCCCGTCGTCACGAAGACGGTCAAACCTTCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCg                                                                                                                       <  2:186829/133‑1 (MQ=255)          gTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAATCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGt                                                                                                                            >  2:497315/1‑123 (MQ=255)          gTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAATCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGt                                                                                                                            <  1:497315/123‑1 (MQ=255)                                        aaaTCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGa                                                                                    <  2:322622/133‑1 (MQ=255)                                               cGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTg                                                                             >  2:304786/1‑133 (MQ=255)                                                              ggcggACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTGCGGATGATATCGAGCTGGATCGTCTTGAACAGCGTCTGGCa                                                              >  2:138940/1‑133 (MQ=255)                                                              ggcggACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCa                                                              >  1:297349/1‑133 (MQ=255)                                                              ggcggACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCa                                                              >  1:388520/1‑133 (MQ=255)                                                              ggcggACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTACGGATGATATCGGGCTGGATCGTCTAGAACAGCGTCTGGCa                                                              >  2:215540/1‑133 (MQ=255)                                                               gcggACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCag                                                             <  2:388520/133‑1 (MQ=255)                                                                      tCTACGCCATTGAGCAGACGGGGCTGTTCTCAAGTCGTTATGTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGaaaa                                                      <  2:148635/133‑1 (MQ=255)                                                                                       aCGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATcc                                     >  2:284951/1‑133 (MQ=255)                                                                                       aCGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATcc                                     >  1:319552/1‑133 (MQ=255)                                                                                        cGGGGCAGTTCTCAGATCGTTATTTTGTTTTGATGGGGCGTTTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCt                                    <  2:174699/133‑1 (MQ=255)                                                                                         ggggCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTc                                   >  1:481293/1‑133 (MQ=255)                                                                                                 tttCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCaa                           <  1:509949/133‑1 (MQ=255)                                                                                                 tttCAGGTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCaa                           >  1:387083/1‑133 (MQ=255)                                                                                                     aGGTCGTTATTTTGTGTTGATGGGGCGTTTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAAcccc                       >  2:514317/1‑133 (MQ=255)                                                                                                      ggTCGTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCa                      <  1:258408/133‑1 (MQ=255)                                                                                                          gTTATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGt                  >  2:40198/1‑133 (MQ=255)                                                                                                            tATTTTGTGTTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTg                <  1:429823/133‑1 (MQ=255)                                                                                                                ttGTGTTGCTGGGGCATCCGTCTCCGCTGGTCGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAgg            <  1:134040/133‑1 (MQ=255)                                                                                                                  gtgtTGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTg          <  2:319552/133‑1 (MQ=255)                                                                                                                      tGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGc      <  1:81255/133‑1 (MQ=255)                                                                                                                      tGATGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGc      >  1:349919/1‑133 (MQ=255)                                                                                                                        aTGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGCTa    >  2:488197/1‑133 (MQ=255)                                                                                                                        aTGGGGCGTCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGCTa    <  1:284951/133‑1 (MQ=255)                                                                                                                          ggggCGTCTGTCACCCCTAGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGCTAcc  <  1:215540/133‑1 (MQ=255)                                                                                                                               |                                                                                                                                GTAACATCTGTTCGAATCCGCGTCGTCAGGAAAACGGTCAAATCTGCGTGGTGGAGAGTCCGGCGGACATCTACGCCATTGAGCAGACGGGGCAGTTTTCAGGTCGTTATTTTGTGTTGATGGGGCATCTGTCACCGCTGGACGGCATCGGTCCGGATGATATCGGGCTTGATCGTCTGGAACAGCGTCTGGCAGAGGAAAAAATCACTGAAGTGATCCTCGCCACCAACCCCACGGTTGAAGGTGAAGCTACC  >  NC_000913/494610‑494863

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