Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,878,855 C→T A28V (GCT→GTT)  yoaK → putative membrane protein YoaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,878,8550CT100.0% 46.2 / NA 17A28V (GCT→GTT) yoaKputative membrane protein YoaK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (9/8);  total (9/8)

ATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTG  >  NC_000913/1878612‑1879082
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
atggtagAAGTCGTTTCTGAAACAACCGATGCATACAGGGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGAGTGGTGATAACATGCGAATCGGTTTTATTTTTCCCGTTGTAATCTTCATTACAGGGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCg                                                                                                                                                                                                                                <  2:176573/248‑1 (MQ=255)
                                                      gATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTg                                                                                                                                                                         <  1:442948/251‑1 (MQ=255)
                                                      gATATTGCAATGCATCAGGAGAAAAAACACAACCAGCAAACACCTTTTGTCGCGCATCCAGCGATACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTCTTGGCGGCTATGTTGCCCCGGGAGCATAAAAGTGAAAACAACCACGATTATTATGATTGGTGTTGCGattattt                                                                                                                                                                         >  2:565585/1‑250 (MQ=255)
                                                               aTGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGACCCGGGAGCATAAAGATGACAAAAACAACGATTATTATGTTGGGTGTGGCGAATATTGTCGTACTct                                                                                                                                                                >  1:351748/1‑250 (MQ=255)
                                                                      aGGAGAAAAAACACAACCAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCTCATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGGCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGa                                                                                                                                                         <  2:345054/251‑1 (MQ=255)
                                                                                acacAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGAtgg                                                                                                                                               >  1:133014/1‑251 (MQ=255)
                                                                                             ccacCTTTGGTCGCGCATCCAGCGCTACATTCCTGAGGCGTAGTTAGAGTCTGTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCTGTCGTATTTTTAGCATGGTGTTTTATTGGCGGCTATGTTGCACCGGGAGTATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCTATTATTGACGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCAcc                                                                                                                                  <  1:93787/250‑1 (MQ=255)
                                                                                                             aTCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCaa                                                                                                                   >  1:569736/1‑250 (MQ=255)
                                                                                                               ccAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGGGCATAAAAATGAAAAAAACAAAGCTTATTATGATGGGTGTGGCGAATAATGTCGTACTCGGCACTGAGCTGGGATGGGGGTAACGGCACCTCTAAAAAATAGCAAAgg                                                                                                               >  1:488243/1‑251 (MQ=255)
                                                                                                                         aTTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCtgtgt                                                                                                     <  1:571138/251‑1 (MQ=255)
                                                                                                                                                         gTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTCTTATGATGGGTGTGGCGATTATTGTCGTACTCGGCCCTGAGCTGGGATGGGGGTAACGTCACCTCTAAAAATTAGCAAAGGCTGCCTGTGTGCAGCCTTTGTTCAATTTAAGCGTTAACtttt                                                                     >  1:348629/1‑251 (MQ=255)
                                                                                                                                                                                 tttttCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCCATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGGTGGGGGAGGCGATTATTGTTGTAATAGGCACTGGGCTTGGATGTGGGTAACGCAACCTCTAAAAAAAAGAAAAGGCCGGCTGTGGGCTGCCTTTGTGCAAATTAAACGTTAACTTTTAATATTACTTTAGGTAAATaaaac                                             >  1:167782/1‑247 (MQ=255)
                                                                                                                                                                                       cGGTTGTAATCTTCATTACAGCGGGCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAAt                                       <  1:100823/251‑1 (MQ=255)
                                                                                                                                                                                        ggtggTAATCTTCATTACAGCGGCCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGGGCATAAAAATGAAAAAAACAAAGATTATTATGGTGGGGGGTGCCGTTATTGTTGTACTCGGCACTGAGCTGGGAGGGTGGGAACGTTCCCTCTAAAAAAAAGCAAAAGGATCCTGGGTGCAGCCTTTGTGCAATTTTAGCGTATACTTTTAAACTTCATGttaataaaaatcaagacaat                                       >  2:319186/5‑231 (MQ=255)
                                                                                                                                                                                                   tCATTTCAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATa                             <  2:133014/249‑1 (MQ=255)
                                                                                                                                                                                                    cATTAAAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGAGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACg                          <  2:103739/251‑1 (MQ=255)
                                                                                                                                                                                                                             tGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTAGGCACTGAGCTGGGGTGGGGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCGGCCAAAATTg  >  1:145810/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                     
ATAGTAGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTG  >  NC_000913/1878612‑1879082

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: