Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 46,790 A→G L328L (CTA→CTG yaaU → putative transporter YaaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091346,7900AG94.7% 37.0 / ‑3.9 19L328L (CTA→CTGyaaUputative transporter YaaU
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (8/10);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

GTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGG  >  NC_000913/46544‑47024
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
gtttgttGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATc                                                                                                                                                                                                                                         <  1:56748/250‑1 (MQ=255)
  ttgttgCCGCCATCTGGACCTGCCAGGTGATTCCAATGTTCGCCATTTACACCGTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTGGGGGCTGATccc                                                                                                                                                                                                                                       <  1:261656/250‑1 (MQ=255)
  ttgttgCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGAGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATccc                                                                                                                                                                                                                                       <  1:143940/250‑1 (MQ=255)
                                           gCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGgtgtt                                                                                                                                                                                             >  1:274083/1‑250 (MQ=255)
                                                        ttGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGGTCTGGCTGGGAGTGATGGCCTTTGCGGTGTATGCCttttttttt                                                                                                                                                                                >  1:205337/1‑247 (MQ=255)
                                                         tGGCCCGCAAAGCTTGGGGTTTTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGGGGTGATTATCCTTTTCTTTATGCTCGGCTGTAGTCCGCCGATGCTGTGGTTAAACACGGCCTGACGTCGTCCATTGTTGATTGGCAGCTTTGCTATGATGACGGTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTctct                                                                                                                                                                                <  2:178428/250‑1 (MQ=255)
                                                              cGCAAATCGTTGGTTTGTTGGGATTGGGTGTTGGCAAAAACGCGGCGCTGGGGAACGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGGTCCCGGATATTGGGGTCTGGCTGGTAGTGGTGGCCTTTGCGGTGTATGCCTTTTTCTTTTGCggg                                                                                                                                                                          >  1:18703/1‑251 (MQ=255)
                                                                                          tgTTGGCAAATACGCGGCTCTGGGGCACGTGGTGATTAGCCTGTTCGTTATGCGCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTGCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCt                                                                                                                                               <  2:272165/249‑1 (MQ=255)
                                                                                          ggTTTGCAATAACGCGGGGCTGGGGAACGGGGTGATTTGCCTGTTCTTTATGCTCGGCTGTATTCCTCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTTGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTa                                                                                                                                              <  2:274083/251‑1 (MQ=255)
                                                                                                                                         ttATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGCCGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATTCCGGGTATGGGGGTCTGGCTGGTAGTGGTGGCCTTTGCGGGGTATGCCTTTTTTTCTGgcggggcggg                                                                                                                                                                     >  1:332727/1‑181 (MQ=255)
                                                                                                                                                                   cTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGGGTATGCCTTTTTCTCTGGCGGGGCGGGGAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGGGGGCGGGATTTTGTCCTTAAGCCGTATTGGCACCAt                                                                     >  1:18192/1‑251 (MQ=255)
                                                                                                                                                                    tgtgGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGGTCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGGAATTTTCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCCTGGGCGGGATTATGTCCTTAAAGCGTATTGGCAccctt                                                                    >  1:245684/1‑248 (MQ=255)
                                                                                                                                                                    ggtgGTTAAACACGGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGTTATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCAtt                                                                    <  2:18703/250‑1 (MQ=255)
                                                                                                                                                                                 tGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCCCTGGCGGCTTTGGGGCTGAACCCGGGAATGGGGATCTGGCTGGGAGTGAAGGGCTTTTGGGGGTATGCCTTTTTCTCTGGCGGGGCGGGTAATTTCCCGTGGCTCTATTCTAATCGACTCTTCCCGCCAGATATACACCCCTCCCCCGCGGGGGGGATTATTTCCTTAAGACGTATTGGCACCATTGTTTTGTCCTggg                                                       >  1:34864/1‑251 (MQ=255)
                                                                                                                                                                                         ggCGTCCATTGTTTATTGGCAGCTTTGCCATGATGACGCTGGTGGTGGCGGTTTGGGGTCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGACTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTACTGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACtgcct                                               <  2:97677/251‑5 (MQ=255)
                                                                                                                                                                                            gTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGGGTATGCCTTTTTCTCTGGCGGGCCGGGGAATTTGCAGTGGGTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCCCTTCCGTGGGCGGGTTTATTTCCTTAAGCCGTATTGGCACCCTTGTTTTGGCCCGGGCACTGCCGATc                                            >  1:136059/1‑251 (MQ=255)
                                                                                                                                                                                              ccATTGTTGATTGGCAGCTTTGCCATGATGACGCGGTCGTTGGCGGTGTTGGGGCTGATCCCTGATATGGGGATCTGGCTGGTAGTTATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCtt                                          <  2:18192/251‑1 (MQ=255)
                                                                                                                                                                                                          ggCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTa                              <  1:235666/251‑1 (MQ=255)
                                                                                                                                                                                                                                       ggcgGTTTTGGGGCTGATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGCCGTATTGGCACCATTGTTTCGACCTGGGCACTGCCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATgggg  <  1:163277/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          
GTTTGTTGCCGCCATCTGGACCTGCCAGGTGATCCCAATGTTCGCCATTTACACCTTTGGCCCGCAAATCGTTGGTTTGTTGGGATTGGGGGTTGGCAAAAACGCGGCACTAGGGAATGTGGTGATTAGCCTGTTCTTTATGCTCGGCTGTATTCCGCCGATGCTGTGGTTAAACACTGCCGGACGGCGTCCATTGTTGATTGGCAGCTTTGCCATGATGACGCTGGCGCTGGCGGTTTTGGGGCTAATCCCGGATATGGGGATCTGGCTGGTAGTGATGGCCTTTGCGGTGTATGCCTTTTTCTCTGGCGGGCCGGGTAATTTGCAGTGGCTCTATCCTAATGAACTCTTCCCGACAGATATCCGCGCCTCTGCCGTGGGCGTGATTATGTCCTTAAGTCGTATTGGCACCATTGTTTCGACCTGGGCACTACCGATCTTTATCAATAATTACGGTATCAGTAACACGATGCTAATGGGG  >  NC_000913/46544‑47024

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