Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 949,276 T→C 20.0% N98S (AAT→AGT)  ycaN ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913949,2760TC20.0% 38.2 / 6.6 20N98S (AAT→AGT) ycaNLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (5/11);  new base C (2/2);  total (7/13)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

GAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGTTGTCCGAT  >  NC_000913/949136‑949417
                                                                                                                                            |                                                                                                                                             
gAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGa                                                                                                                                      <  2:226981/150‑1 (MQ=255)
gAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGa                                                                                                                                      <  2:211736/150‑1 (MQ=255)
gAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGa                                                                                                                                      <  2:211738/150‑1 (MQ=255)
                              tatCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAAtcat                                                                                                        >  1:139400/1‑150 (MQ=255)
                               atCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAAtcatt                                                                                                       <  1:191161/150‑1 (MQ=255)
                                                           tCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCACTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTtcg                                                                           >  1:182224/1‑150 (MQ=255)
                                                           tCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCACTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTtcg                                                                           >  1:48635/1‑150 (MQ=255)
                                                                        cAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGgc                                                              <  1:121345/150‑1 (MQ=255)
                                                                          ggATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGgcgc                                                            >  1:84821/1‑150 (MQ=255)
                                                                          ggATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGgcgc                                                            >  1:84823/1‑150 (MQ=255)
                                                                                    gtgtGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCAt                                                  >  2:1191/1‑150 (MQ=255)
                                                                                                                        gTGCAGCTACCCGGGCGGCACTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTc              <  2:48635/150‑1 (MQ=255)
                                                                                                                        gTGCAGCTACCCGGGCGGCACTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTc              <  2:182224/150‑1 (MQ=255)
                                                                                                                         tGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCg             <  2:139191/150‑1 (MQ=255)
                                                                                                                             gCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTAATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGtt         >  1:26856/1‑150 (MQ=255)
                                                                                                                                aCCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGTTGTc      <  1:141064/150‑1 (MQ=255)
                                                                                                                                    gggCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGTTGTCCGAt  <  2:66961/150‑1 (MQ=255)
                                                                                                                                    gggCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGTTGTCCGAt  <  2:66963/150‑1 (MQ=255)
                                                                                                                                    gggCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGTTGTCCGAt  <  2:9035/150‑1 (MQ=255)
                                                                                                                                    gggCGGCATTGATTGTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCCTCTCATCTACCCTGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATGAGAACCCGCTTCCGTTAACGCTACGCTTCGGGTTGTCCGAt  <  2:95518/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GAACGCCAGCGTCAAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAACCAGTAAAGACATTAAAAATATACGTGCAGCTACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCATAGAGATTAGAACCCGCTTCCGTTAACGATACGCTTCGGGTTGTCCGAT  >  NC_000913/949136‑949417

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