Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,070,763 G→T R65S (CGT→AGT)  rutE ← putative malonic semialdehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,070,7630GT96.7% 99.9 / ‑2.8 30R65S (CGT→AGT) rutEputative malonic semialdehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (17/12);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCAT  >  NC_000913/1070686‑1070838
                                                                             |                                                                           
gTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTc                                                                          >  1:217451/1‑81 (MQ=255)
gTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTc                                                                          >  1:217444/1‑81 (MQ=255)
     aGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTc                                                                     >  1:217445/1‑81 (MQ=255)
         gacgaTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTc                                                                 >  2:217439/1‑81 (MQ=255)
            gaTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGc                                                              >  1:217461/1‑81 (MQ=255)
             aTAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGcc                                                             >  1:217446/1‑81 (MQ=255)
              tAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCg                                                            <  1:217408/81‑1 (MQ=255)
                 cGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGc                                                         <  1:217409/81‑1 (MQ=255)
                 cGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGc                                                         >  2:217440/1‑81 (MQ=255)
                    tGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCggg                                                      >  2:217458/1‑81 (MQ=255)
                      aCGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTa                                                    <  2:217410/81‑1 (MQ=255)
                       cGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTaa                                                   <  2:217411/81‑1 (MQ=255)
                         ggCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTaaa                                                  >  1:217449/1‑80 (MQ=255)
                          gcgcGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAaca                                                >  1:217462/1‑81 (MQ=255)
                              ggTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGAt                                            <  1:217412/81‑1 (MQ=255)
                                tCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATcc                                          >  1:217454/1‑81 (MQ=255)
                                    ggTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGc                                      <  1:217413/81‑1 (MQ=255)
                                         tttGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCg                                    >  1:217456/1‑78 (MQ=255)
                                             cAGATTGCCGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGaa                               >  1:217448/1‑79 (MQ=255)
                                                     cGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCt                      >  2:217460/1‑80 (MQ=255)
                                                     cGCTGGAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTg                     <  2:217414/81‑1 (MQ=255)
                                                         ggAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGt                 >  2:217468/1‑81 (MQ=255)
                                                         ggAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGt                 >  2:217472/1‑81 (MQ=255)
                                                          gAAAGTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTc                >  1:217463/1‑81 (MQ=255)
                                                              gTGCCGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGcc            <  2:217415/81‑1 (MQ=255)
                                                                  cGGGCGCAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCa        <  1:217416/81‑1 (MQ=255)
                                                                   gggCGCAGACTTTCTGTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCAt       <  2:217417/81‑1 (MQ=255)
                                                                     gcgcAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCAttt     >  1:217457/1‑81 (MQ=255)
                                                                        cAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTtcat  <  1:217418/81‑1 (MQ=255)
                                                                        cAGACTTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATACGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTtcat  <  1:217419/81‑1 (MQ=255)
                                                                             |                                                                           
GTCCCAGGCGACGATAGCGGTGACGGGCGCGGTCAGGGTTTTTTGCAGATTGCCGCTGGAAAGTGCCGGGCGCAGACGTTCTTTTCCTTCTGCCGTGCGGGTAAACACGATCCGTGCCGGAGAACAGTTAGCTGATGTCGGCCCCCATTTCAT  >  NC_000913/1070686‑1070838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: