Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,563,534 A→G C24R (TGT→CGT)  yjiK ← uncharacterized protein YjiK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,563,5340AG96.4% 91.5 / NA 28C24R (TGT→CGT) yjiKuncharacterized protein YjiK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (12/15);  minor base C (0/1);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTA  >  NC_000913/4563405‑4563659
                                                                                                                                 |                                                                                                                             
ccGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTacgcac                                                                                                                            >  2:455666/1‑133 (MQ=255)
   aCCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCATTGCCGTCAATGGTGGCGTGGTAATTCTGGAATGATGCTTCATGATTGCGTTTATGGGCGCAAGATTGCACGAAGAAAGTACGCACttc                                                                                                                         >  2:509782/1‑131 (MQ=255)
        gTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACaa                                                                                                                    <  1:90352/133‑1 (MQ=255)
        gTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGGGGCGGGGTAATTCTGAAATGATGCTGCATGAGTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACaa                                                                                                                    >  2:507146/1‑133 (MQ=255)
                agagATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCaaaa                                                                                                            <  2:401656/133‑1 (MQ=255)
                     tattatTCCTTATGCATGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAAt                                                                                                       <  1:195007/133‑1 (MQ=255)
                          ttCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACaa                                                                                                  <  2:352709/133‑1 (MQ=255)
                           tCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACaat                                                                                                 >  1:480550/1‑133 (MQ=255)
                            cGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACaata                                                                                                >  2:472197/1‑133 (MQ=255)
                            cGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACaata                                                                                                >  1:83499/1‑133 (MQ=255)
                              ttATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAAAAAGTACGTACTCCCAGAATGGCAAAAAGAATCACaataat                                                                                              >  2:178937/1‑133 (MQ=255)
                               tATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGGAAGTACGCACTGACACGATGGCAAATAGAAACACAATAATg                                                                                             >  2:139129/1‑133 (MQ=255)
                                 tGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGga                                                                                           <  1:472197/133‑1 (MQ=255)
                                  gCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGgag                                                                                          <  1:95188/133‑1 (MQ=255)
                                    cTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGAt                                                                                        >  1:69277/1‑133 (MQ=255)
                                                          ggtggCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACtt                                                                  <  2:69277/133‑1 (MQ=255)
                                                               cGTGGTAATTCTG‑AATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTc                                                            >  1:375988/1‑133 (MQ=255)
                                                                gTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTc                                                            <  1:506824/133‑1 (MQ=255)
                                                                        tCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAGCGAAATACTTTTTGTCACGGTAAt                                                    <  2:480550/133‑1 (MQ=255)
                                                                             aaTGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATTCTTTTTATCACGGAAATCCCtt                                               >  1:165601/1‑133 (MQ=255)
                                                                                     tgcATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCaa                                       <  1:509782/133‑1 (MQ=255)
                                                                                       cATGATTGATTTTACCGGAGCACGATTGCATGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAAt                                     <  1:139129/133‑1 (MQ=255)
                                                                                         tGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAATTCa                                   <  2:5954/133‑1 (MQ=255)
                                                                                            ttGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTATCACGGTGATCGCTATAAGGCGAATCAATa                                >  1:185401/1‑133 (MQ=255)
                                                                                              gCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAg                              <  1:507146/133‑1 (MQ=255)
                                                                                                 tttACGGGCGCAAGATTGCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGata                           >  1:351199/1‑133 (MQ=255)
                                                                                                     cGGGCGCAGGATTGCCCGAAGACAGTACGCCCTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCga                       <  1:166988/133‑1 (MQ=255)
                                                                                                     cGGGCGCAAGATTGCCCGAAGAAAGTACCCACCGCCACAATCGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCga                       <  2:165601/133‑1 (MQ=255)
                                                                                                                  gCACGAAGAAAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAg          <  2:351199/133‑1 (MQ=255)
                                                                                                                        agaaAGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCt    >  2:25187/1‑133 (MQ=255)
                                                                                                                          aaaGTACGCACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTa  >  2:55772/1‑133 (MQ=255)
                                                                                                                                 |                                                                                                                             
CCGACCAGGTCAGCGAAGAGATATTATTCGTTATGCCTGCAATTTCTTTGCCGTCAATGGTGGCGTGGTAATTCTGAAATGATGCTGCATGATTGCTTTTACGGGCGCAAGATTGCACGAAGAAAGTACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTA  >  NC_000913/4563405‑4563659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑