Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 630,265 T→C C383C (TGT→TGC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb630,2650TC98.0% 169.1 / ‑4.6 51C383C (TGT→TGCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (25/25);  minor base G (1/0);  total (26/25)
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.97e-01

GCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGG  >  W3110S.gb/630198‑630327
                                                                   |                                                              
gCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCat                                                              >  1:4092012/1‑70 (MQ=255)
 cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCatca                                                            >  1:3827677/1‑71 (MQ=255)
 cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCatca                                                            >  1:3160188/1‑71 (MQ=255)
 cTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCatca                                                            >  1:2264064/1‑71 (MQ=255)
   ttATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCatcatc                                                          >  1:4415050/1‑71 (MQ=255)
           tACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCggg                                                  >  1:1312028/1‑71 (MQ=255)
           tACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCggg                                                  >  1:2136511/1‑71 (MQ=255)
               gTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCgtg                                              <  1:3332086/71‑1 (MQ=255)
               gTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCgtg                                              <  1:2707558/71‑1 (MQ=255)
                  ggATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTa                                            <  1:4473475/70‑1 (MQ=255)
                     tGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTAttt                                         >  1:1900714/1‑70 (MQ=255)
                     tGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCAGGGCGTGTAtt                                          >  1:31/1‑69 (MQ=255)
                       ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGc                                      >  1:1150615/1‑71 (MQ=255)
                       ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGc                                      >  1:3607811/1‑71 (MQ=255)
                       ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGc                                      >  1:2024816/1‑71 (MQ=255)
                       ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGc                                      >  1:1261737/1‑71 (MQ=255)
                       ttAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTAATTTg                                       >  1:1474774/1‑70 (MQ=255)
                         aaTGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTAttt                                         >  1:2802974/1‑66 (MQ=255)
                           tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATg                                  <  1:4326785/71‑1 (MQ=255)
                           tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATg                                  <  1:2707093/71‑1 (MQ=255)
                           tGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATg                                  <  1:3046771/71‑1 (MQ=255)
                            ggAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATg                                  <  1:798884/70‑1 (MQ=255)
                             gAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGa                                 >  1:2080170/1‑70 (MQ=255)
                              aaTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGa                                                        <  1:2101074/46‑1 (MQ=255)
                              aaTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGa                                                        <  1:2665176/46‑1 (MQ=255)
                               aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTAttt                                         >  1:1207193/1‑60 (MQ=255)
                               aTCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTTCCATGaa                                >  1:3420877/1‑69 (MQ=255)
                                tCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTa                                            <  1:155179/56‑1 (MQ=255)
                                tCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTa                                            <  1:1468678/56‑1 (MQ=255)
                                tCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATg                                  <  1:3457725/66‑1 (MQ=255)
                                  cTTCGTGGCGATTATGGCTCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGcc                           >  1:1788513/1‑71 (MQ=255)
                                  cTTCGTGGCGATTATGGCGCTGGTTTCCTCCTGCATCATCGATCCGGGCGTGTATTTTTcc                                     >  1:2828668/1‑61 (MQ=25)
                                     cGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGccc                          <  1:1093948/69‑1 (MQ=255)
                                     cGTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATAGATCCGGGCGTGTATTTTGCCATGAACAGccc                          <  1:1883934/69‑1 (MQ=255)
                                      gTGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGAt                       <  1:2153917/71‑1 (MQ=255)
                                       tGGCGATTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGc                                      <  1:3425648/55‑1 (MQ=255)
                                           gATTATGGCGCTGGTTTCCGCCTGCATCATAGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCgg                  <  1:1549784/71‑1 (MQ=255)
                                            aTTATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCgg                  <  1:813524/70‑1 (MQ=255)
                                             ttATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGaa                                <  1:2871852/55‑1 (MQ=255)
                                              tATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTAttt                                         >  1:3554165/1‑45 (MQ=255)
                                              tATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGcc                                     >  1:2842794/1‑49 (MQ=255)
                                              tATGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGc               >  1:227373/1‑71 (MQ=255)
                                                tGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCt              <  1:4595918/70‑1 (MQ=255)
                                                tGGCGCTGGTTTCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTg             <  1:1977126/71‑1 (MQ=255)
                                                           tCCGCCTGGATCATCGATCCGGGCGTGTATTTTGCCa                                    >  1:4255184/1‑37 (MQ=255)
                                                           tCCGCCTGCATCATCGATCCGGGGGTGTATTTTGCCa                                    >  1:3647908/1‑37 (MQ=38)
                                                           tCCGCCTGCATCATCGATCCGGGCGTGTATTTTGCCa                                    >  1:2122313/1‑37 (MQ=255)
                                                            ccGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAgg  <  1:2247330/70‑1 (MQ=255)
                                                            ccGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAgg  <  1:246059/70‑1 (MQ=255)
                                                            ccGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAgg  <  1:90984/70‑1 (MQ=255)
                                                            ccGCCTGCATCATCGATCCGGGCGTGTATTTTGCCATGAACACCCCGATGGCGGTGCTGGCTCCGGCAgg  <  1:4618335/70‑1 (MQ=255)
                                                                   |                                                              
GCTTTATCGGCTACGGTGGGATGTTAATGGAATCCTTCGTGGCGATTATGGCGCTGGTTTCCGCCTGTATCATCGATCCGGGCGTGTATTTTGCCATGAACAGCCCGATGGCGGTGCTGGCTCCGGCAGG  >  W3110S.gb/630198‑630327

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