Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,943 T→C S379S (AGT→AGC)  yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,9430TC94.7% 44.0 / ‑5.5 19S379S (AGT→AGC) yaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (10/8);  minor base G (1/0);  total (11/8)
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.94e-01

CGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAAT  >  NC_000913/46702‑47170
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                   
cgATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATTCCGGGTATGGGGATCTGGCTGGGAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGGAATTTGCAGTGGGTCTTTCCTAATGAACTCTTTCCGGCAGAAATCCCCGCCTCTGCCGTGGGGGTGATTTTGTCCTTAAACCGTATTGGc                                                                                                                                                                                                                            >  2:133136/1‑251 (MQ=255)
          gTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCGGATCCCGGATATGGGGATCTGGCTGGTAGGG‑ATGGCCTTTGTGGGGTATGCCTTTTTCTCTGGCGGGCCGGGGGATTTGCAGTGGCTCTTTCCTAATGTACTCTTCCCCACAGAGATCCGCGCCCCTGCCGTGGGCGTGGTTATGTCCTTAAACCGTATTGGCACCATTGttt                                                                                                                                                                                                                  >  2:210059/1‑251 (MQ=255)
                              gTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTTGGGCTGATTCCGGATAGTGGGATCTGGCTGGGAGTGTATGGC‑TTTGCGGGGGTTGCCTTTTTCTCTGGGGGGGCGGGGGATTTGGCGTGGGTCTATACTTATGAACTCATCCCGACCGATATACGCGCCTCTGGCGGTGGCGGGATTGTGTCCTTAATGCGTATTGGCACCAtttttttggcctcggtacggcccatc                                                                                                                                                                                              >  2:130699/1‑231 (MQ=255)
                                     gtttaTTGGCAGCTTTGCCATGATGACGCTGGTGCGGGCGGTTTTGGGTCTGACCCCGGATTTGGGGATCGTGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGGGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTTGGCGTGATTATCTCCTTAAGCCGTATTGGCACCAATGTTTCGACCTGGGCACTGCAGATCTTTATCa                                                                                                                                                                                       <  2:165958/247‑1 (MQ=255)
                                                          gatgaCGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAaca                                                                                                                                                                  <  1:133136/251‑1 (MQ=255)
                                                                        tggcgGTTTTGGGGCTGATACCGGTTATGGGGATCTGGCTGGGAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTTCAGTGGCTCTATCCTTATGAACTCTTCCCGACAGATATCCGCGCCTCTTCCCTGGGCGTGATTATTTCCTTAAGCCGTATTGGCCCCCTTGTTTCGTCCTGGGCCCTGCCCGTCTTTTTCAATTATTACGGTTTCAGTAACCCGATGCTAAtcgtg                                                                                                                                                    >  2:9893/1‑247 (MQ=255)
                                                                            ggTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATg                                                                                                                                                       <  2:2754/244‑1 (MQ=255)
                                                                            ggTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGCTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATg                                                                                                                                                       >  1:2754/1‑244 (MQ=255)
                                                                                           cccGGATGTGGGGAGCTGGCTGGTAGTG‑TTGGCATTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGGGGCTCTATCCTAATGAACTCTTCCCTACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGATCTGGGCACTGCCGATCTTTATCAATAATTACGGTAACAGTTTCACGATGCTAAGGGGGGCTAGTATCTCGCTGttt                                                                                                                                  <  2:174278/250‑1 (MQ=255)
                                                                                              ggATATGGGGATCTGGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCt                                                                                                                              <  1:60832/251‑1 (MQ=255)
                                                                                                                             cTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGATCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCgg                                                                                                   >  2:161915/1‑248 (MQ=255)
                                                                                                                                                                                          gAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGt                                    >  1:211104/1‑250 (MQ=255)
                                                                                                                                                                                                      aCAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGGCCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGAt                       <  1:210059/251‑1 (MQ=255)
                                                                                                                                                                                                         gATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGCTTTCCGTAGCGTTTGCCCCGGAGTCTCGAGGGATGTCACTGGCGCACACCAGCAATATGACGATCCGCGTGCAGCGAATGGGGTAAATTGTTCAGCTTtc                    >  1:39320/1‑251 (MQ=255)
                                                                                                                                                                                                            atCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTtctct                 <  1:121719/251‑1 (MQ=255)
                                                                                                                                                                                                               cgcgCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCtt              <  1:27564/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGCGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATc     >  1:170565/1‑250 (MQ=255)
                                                                                                                                                                                                                          cGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCCCTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGCGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCa    >  2:222133/1‑250 (MQ=255)
                                                                                                                                                                                                                           gTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCCGTAACACGATGCTACTGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAAACCAGCAATATGACGATCCGCGGGCAGTGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATccat  >  1:259856/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                   
CGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTG‑ATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGGGCGGGTATCTCGCTGTTTGGCTTGTTGATTTCCGTAGCGTTTGCCCCGGAGACTCGAGGGATGTCACTGGCGCAGACCAGCAATATGACGATCCGCGGGCAGAGAATGGGGTAAATTGTTCAGATTTCTCTCTTTTCTGAATCAAT  >  NC_000913/46702‑47170

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