Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,083,312 G→A 12.3% G325E (GGG→GAG)  rtcB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,083,3120GA12.3% 164.2 / 12.0 57G325E (GGG→GAG) rtcBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/33);  new base A (5/2);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 9.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

TCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAAAGCG  >  W3110S.gb/4083245‑4083382
                                                                   |                                                                      
tCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGc                                                                      <  1:3358569/70‑1 (MQ=255)
    cgcgTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTg                                                                    <  1:2156037/68‑1 (MQ=255)
       gTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAAt                                                              <  1:754418/71‑1 (MQ=255)
        tGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATg                                                             <  1:3833147/71‑1 (MQ=255)
        tGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATg                                                             <  1:218885/71‑1 (MQ=255)
          cTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGaa                                                           <  1:4057105/71‑1 (MQ=255)
              tCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGt                                                       <  1:3608614/71‑1 (MQ=255)
                aaTATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTc                                                      <  1:4218013/70‑1 (MQ=255)
                aaTATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCg                                                     >  1:1331771/1‑71 (MQ=255)
                  tatGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGtt                                                   <  1:4492400/71‑1 (MQ=255)
                  tatGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGt                                                    >  1:248188/1‑70 (MQ=255)
                   atGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGtt                                                   <  1:4497670/70‑1 (MQ=255)
                    tGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTggg                                                                  >  1:3978433/1‑54 (MQ=255)
                      gAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTg                                                >  1:4176012/1‑70 (MQ=255)
                      gAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGt                                               >  1:417638/1‑71 (MQ=255)
                       aaTTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGtt                                              <  1:1846679/71‑1 (MQ=255)
                       aaTTATTCCCGGTTCGATGCGAGCAAAAAGCTTTATCGTCCGTGGGCTggg                                                                  >  1:2194682/1‑51 (MQ=255)
                         ttattCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCg                                                     <  1:2478668/62‑1 (MQ=255)
                           attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTgg                                                                        >  1:1727973/1‑41 (MQ=255)
                           attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTgg                                                                        >  1:699777/1‑41 (MQ=255)
                           attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGTCTGGGAAATGAAGAGTCGTTCTGTTCGt                                           <  1:652732/70‑1 (MQ=255)
                           attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGt                                           <  1:3384823/70‑1 (MQ=255)
                           attCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTg                                          >  1:1394229/1‑71 (MQ=255)
                            ttCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGaa                                                           <  1:3496611/53‑1 (MQ=255)
                            ttCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAag                                                          >  1:1529088/1‑54 (MQ=255)
                            attccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTg                                          >  1:2332864/4‑70 (MQ=255)
                              cccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAg                                       >  1:1937086/1‑71 (MQ=255)
                               ccGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGc                                         <  1:3363190/68‑1 (MQ=255)
                                 ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa                                    <  1:326788/71‑1 (MQ=255)
                                 ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa                                    <  1:2138170/71‑1 (MQ=255)
                                 ggTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCa                                    >  1:3855970/1‑71 (MQ=255)
                                  gTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCAc                                   <  1:1106931/71‑1 (MQ=255)
                                  gTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCAc                                   <  1:2020348/71‑1 (MQ=255)
                                    tCGATGGGAGCAAAAAGCTTTATCGTCCGTGAGCTTGGAAATGa                                                            >  1:2366826/1‑44 (MQ=38)
                                     cGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGt                                           >  1:2610400/1‑60 (MQ=255)
                                      gATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTc                                             <  1:1627950/57‑1 (MQ=255)
                                      gATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGc                                      <  1:3030926/64‑1 (MQ=255)
                                       aTGGGAGCAAAAAGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGc                              <  1:3746883/71‑1 (MQ=255)
                                        tGGGAGCAAAAAGATTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGcc                             >  1:1731608/1‑71 (MQ=255)
                                          ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg                           <  1:4311830/71‑1 (MQ=255)
                                          ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg                           >  1:173958/1‑71 (MQ=255)
                                          ggAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg                           >  1:637433/1‑71 (MQ=255)
                                            aGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggc                         >  1:1038872/1‑71 (MQ=255)
                                             gCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcg                        >  1:279972/1‑71 (MQ=255)
                                               aaaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcgg                       <  1:3438424/70‑1 (MQ=255)
                                                aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcggg                      >  1:3798578/1‑70 (MQ=255)
                                                aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGt                     <  1:2683885/71‑1 (MQ=255)
                                                aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGt                     <  1:1752232/71‑1 (MQ=255)
                                                aaaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTGCTGTTCGTGCAGCCACGGTGCCGGGCGGGt                     <  1:2689315/71‑1 (MQ=255)
                                                aaaaGCTTTATCGTCCGTGAGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgg                           <  1:3063136/65‑1 (MQ=255)
                                                 aaaaCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGt                     >  1:4185992/1‑70 (MQ=255)
                                                 aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTg                                          <  1:4548841/49‑1 (MQ=255)
                                                 aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTa                    <  1:1334683/71‑1 (MQ=255)
                                                 aaaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCAGGCGGGTa                    <  1:3806061/71‑1 (MQ=255)
                                                  aaGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcggg                          >  1:2185653/1‑64 (MQ=255)
                                                   aGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTa                    >  1:203339/1‑69 (MQ=255)
                                                   aGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAAt                  <  1:3091469/71‑1 (MQ=255)
                                                        tATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGt                                               <  1:3519576/37‑1 (MQ=255)
                                                          tCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCcgggcggg                      <  1:1770192/60‑1 (MQ=255)
                                                                   ggCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAAAGCg  <  1:4320116/71‑1 (MQ=255)
                                                                   |                                                                      
TCTGCGCGTGCTGGTCAATATGGAATTATTCCCGGTTCGATGGGAGCAAAAAGCTTTATCGTCCGTGGGCTGGGAAATGAAGAGTCGTTCTGTTCGTGCAGCCACGGTGCCGGGCGGGTAATGAGCCGAACTAAAGCG  >  W3110S.gb/4083245‑4083382

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