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,8550CT96.3% 81.5 / ‑5.9 27A28V (GCT→GTT) yoaKputative membrane protein YoaK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/12);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGACT  >  NC_000913/1878617‑1879094
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
agaaagcGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTTCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTGTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAg                                                                                                                                                                                                                                      <  1:431203/244‑1 (MQ=255)
 gAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAAGCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCa                                                                                                                                                                                                                                    <  2:214325/251‑1 (MQ=255)
       gATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTTCCCCGGGAGCATaa                                                                                                                                                                                                                                 >  2:371264/1‑248 (MQ=255)
              aaaCAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCTCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGaaaaa                                                                                                                                                                                                                       <  1:297249/251‑1 (MQ=255)
                    ccGATGCAGACAGGGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACAGGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGACGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAaacaac                                                                                                                                                                                                                 <  1:201609/251‑1 (MQ=255)
                     cGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCTCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGGTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACg                                                                                                                                                                                                                <  2:231921/251‑1 (MQ=255)
                                    ccggCATGATGACGATATTGCGATGTATCAGGAGATAAAACGCAAACATGAAACACCTCTTGGCGCGCATCCAGCGCGACATTCCTGAGGCGTATGCACATCCTGATGATTGGTTAGAACATGAGAATGGGTATTATTTTTCCGGTTGTAATCGTGATTACATCGGTCGTATGTTTAGCATCGTTTTTTATTGGCGGCGGAGCTTCCCCGGGCTCATAAAGATGAAAAAAACAACGATTATTATGATGGgt                                                                                                                                                                                                 <  1:363973/248‑1 (MQ=255)
                                                         gATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTc                                                                                                                                                                            >  1:262698/1‑251 (MQ=255)
                                                                      aaaaaaCACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGGGCATAAAGATGAAAAAAAACACGATTATTATGATGGGGGGTGCGTTTATTTGCGTACTCGGCACTggggtgg                                                                                                                                                               >  1:15314/1‑245 (MQ=255)
                                                                                                 gTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGAACTCGGCACTGAACTGGGGTGGGGGGAACGTCACCTCTaaaaa                                                                                                                                     >  2:107205/1‑250 (MQ=255)
                                                                                                    gcgcATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTTAAAAATa                                                                                                                                  >  1:242126/1‑250 (MQ=255)
                                                                                                      gcATCCAGCGATACATTCCTGAGCCGTATTCACATCCTTTTGATTGGCGATAACATGCGCATCGGTATTATTTTTCCGGTTGAAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGTCCCGGGGGCATAAAGATGCAAAAAACAACAGTTATTATGATGGGGGGGGGGCTTTTTTTCGTACTCGGCCCTGAGATGGGATGGTGGTAACGGCCCCTCTAAAAAAtcac                                                                                                                                >  2:256449/1‑247 (MQ=255)
                                                                                                                                 aTTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGTGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCTCCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTtgtg                                                                                                     <  1:312802/250‑1 (MQ=255)
                                                                                                                                      caTCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGTCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCGCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAAttt                                                                                               <  1:256449/251‑1 (MQ=255)
                                                                                                                                      caTCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAAttt                                                                                               <  2:262698/251‑1 (MQ=255)
                                                                                                                                             ttGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGt                                                                                         <  1:186926/250‑1 (MQ=255)
                                                                                                                                               gATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGGGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCCGCCTTTGTGCAATTTAAGAGTTa                                                                                       >  2:419224/1‑250 (MQ=255)
                                                                                                                                                           cATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTc                                                                          <  1:571448/251‑1 (MQ=255)
                                                                                                                                                           cATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGAATAAAGATGAAAAAAACAACGATTATTTTGATGGGGGTGGCGATTTTTGTCGTACTCGGCCCTGAGCTGGGAGGGGGGTAAAGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGAAGCCTTTGTGACATTTAAACGGTAACTTTTAATCttt                                                                          >  2:571448/1‑250 (MQ=255)
                                                                                                                                                             tGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGGGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCt                                                                        >  1:404115/1‑251 (MQ=255)
                                                                                                                                                                  aTCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCTGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGAGTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAAGCTTCCTGTAGa                                                                   <  2:80198/251‑1 (MQ=255)
                                                                                                                                                                                     ttGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGc                                                >  1:437902/1‑251 (MQ=255)
                                                                                                                                                                                       gTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGGTGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCAc                                              >  1:78909/1‑251 (MQ=255)
                                                                                                                                                                                        tAATCTACATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTAAGTTGTACCGTGATTATAAAGATGTAATAAATACCGTTTATTTTAATGGTGGTGGCGAATAATGTCGTACACGGCACTGAACTGGGGTGGTGGTAACGTGACCTCTATAAAATAGCCAAGGCTGTCCGTGGGCAGCCTTTTTGCCATTTAAACATTAACTATTAGTCTTTCTGTCGATAAATataccaaaaagctcaac                                             >  2:171987/1‑235 (MQ=255)
                                                                                                                                                                                                           cGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGaa                           <  1:86780/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ttGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGAc   >  2:247160/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tGGCGGCTATGTTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGACt  >  1:433019/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
AGAAGTCGATCCTGAAACAACCGATGCAGACAGTGCCCTGCATGCTGCCGATATTGCGATGTATCAGGAGAAAAAACACAAACAGAAAACACCTTTTGTCGCGCATCCAGCGCTACATTCCTGAGGCGTATTCACATCCTTTTGATTGGTGATAACATGCGAATCGGTATTATTTTTCCGGTTGTAATCTTCATTACAGCGGTCGTATTTTTAGCATGGTTTTTTATTGGCGGCTATGCTGCCCCGGGAGCATAAAGATGAAAAAAACAACGATTATTATGATGGGTGTGGCGATTATTGTCGTACTCGGCACTGAGCTGGGATGGTGGTAACGTCACCTCTAAAAAATAGCAAAGGCTGCCTGTGTGCAGCCTTTGTGCAATTTAAGCGTTAACTTTTAATCTTCCTGTAGATAAATAGCACGACAATCGCACCAATAACGGCAACCACGAAGCTGCCAAAATTGAAGCCATCGACT  >  NC_000913/1878617‑1879094

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