Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,616 C→A 75.4% R2S (CGT→AGT)  ydiE → PF10636 family protein YdiE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,6160CA75.4% 37.2 / 5.1 25R2S (CGT→AGT) ydiEPF10636 family protein YdiE
Reads supporting (aligned to +/- strand):  ref base C (2/4);  new base A (13/5);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGC  >  NC_000913/1789480‑1789751
                                                                                                                                        |                                                                                                                                       
gAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTtata                                                                                                                                   <  1:157106/143‑1 (MQ=255)
         cccATTTCTGACGGAATGGGCATTTGTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCa                                                                                                                          <  2:30215/143‑1 (MQ=255)
          ccATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCGATTATTATTAATTTTTATGAGTTATACGGATAGCAg                                                                                                                         <  1:314904/143‑1 (MQ=255)
              tCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGaaaa                                                                                                                     <  2:178330/143‑1 (MQ=255)
                  gACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTATTTTTTATGAGTTATACGGATAGCCGAAAACTcc                                                                                                                 >  1:216690/1‑142 (MQ=255)
                          gggCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGaa                                                                                                         <  2:73351/143‑1 (MQ=255)
                               tgtCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGa                                                                                                    <  2:370937/141‑1 (MQ=255)
                                    gcgcAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAAAGGATGcga                                                                                               >  1:400964/1‑141 (MQ=255)
                                    gcgcAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAAATTTCAGTTATCGGACGGTTAGCAGAGGACTCACGCCTGAAACGGATGCCa                                                                                               >  1:220735/1‑143 (MQ=255)
                                           ttgttgTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGCATGCCAATCACaa                                                                                        >  2:79027/1‑143 (MQ=255)
                                               tgtCTTATCATCATATTTCTGTTGGCTCTGGTCTGCTAGAATATTGATAATATGAATCATTGTTATATTAATTTTTATTAATTTTAATGCGTCATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCTAATCACAAGAcc                                                                                    <  2:26504/143‑1 (MQ=255)
                                                   ttctCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGCAACGGAGGCCAATCACAAGACCGCtc                                                                                >  2:284609/1‑142 (MQ=255)
                                                   ttctCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCtt                                                                                >  2:116038/1‑143 (MQ=255)
                                                      tcAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTccc                                                                             >  2:290132/1‑143 (MQ=255)
                                                      tcAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAAAAAGAATCATTGTTATAACAATTACTAATAAGTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCAGAAACGGATGCCAATCACAAGACCGCTAccc                                                                             >  1:25611/1‑143 (MQ=255)
                                                        aacaaATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAAGTAGAATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGAATTCCAATCACAAGACCGATCCCCCg                                                                           >  2:194565/1‑143 (MQ=255)
                                                                     ttgctCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGc                                                              <  2:48014/143‑1 (MQ=255)
                                                                                  ccATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTccagcca                                                 >  1:7492/1‑143 (MQ=255)
                                                                                          ttGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACGCTGt                                         >  1:45214/1‑143 (MQ=255)
                                                                                            gATAATAAGAATCATTGTTATAAACATTGTTCGCAAGCAGTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGCATTTCCAGCCAGACACCGtga                                       >  2:217595/1‑141 (MQ=255)
                                                                                                  aaGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCGATGCGGCGAATTGCCAGCCAGAGACTGGTAGGTGCg                                 >  2:75145/1‑143 (MQ=255)
                                                                                                  aaGAATCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAAATCATCAGTGAGACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCg                                 >  1:68469/1‑143 (MQ=255)
                                                                                                      aTCATTGTTATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATg                             <  1:414586/143‑1 (MQ=255)
                                                                                                             ttATATCAATTATTATTAATTTTTATGAGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTCTTAGCTCCGGATGGCAAACt                      >  2:284402/1‑143 (MQ=255)
                                                                                                                                 ttttATGAGTTATACGGAGAGCAGAAAACTCACGCCTTAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACGCGCTACGGTCCGGAAGGCGAACGGATCATCGGTCGTGACGGGc  >  1:388921/1‑143 (MQ=255)
                                                                                                                                        |                                                                                                                                       
GAATGCGTGCCCATTCCTGACGGAATGGGCATTTCTGCGCAACTTGTTGTCTTCTCAACAAATTACTGCTTGCTCTGGTCAGCCATAATATTGATAATAAGAATCATTGTTATATCAATTATTATTAATTTTTATGCGTTATACGGATAGCAGAAAACTCACGCCTGAAACGGATGCCAATCACAAGACCGCTTCCCCGCAGCCTATTCGGCGAATTTCCAGCCAGACACTGTTAGGTCCGGATGGCAAACTGATTATCGATCATGACGGGC  >  NC_000913/1789480‑1789751

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