Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,777,419:1 +A 100% coding (3964/4581 nt) ypjA ← adhesin‑like autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,777,4191.A88.1% 107.3 / 11.9 42coding (3964/4581 nt)ypjAadhesin‑like autotransporter
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base A (14/23);  total (16/26)
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.02e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CACTTTCATGTTCCCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAATAT  >  W3110S.gb/2777350‑2777484
                                                                      |                                                                  
cACTTTCATGTTCCCAACTGGCATAGCCGCCCAGAGAATGACTGCCCACACTGCCATGTCCTCCGCGATc                                                                     >  1:2494885/1‑70 (MQ=255)
cACTTTCATGTTCCCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATc                                                                     >  1:1467499/1‑70 (MQ=255)
     tCATGTTCCCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGACATGTCCTCCGCGATC‑‑aaaaac                                                             >  1:1894528/1‑69 (MQ=255)
       aTGTTCCCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCCATC‑‑aaaaacc                                                            <  1:322469/70‑4 (MQ=255)
              cAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATAtgtg                                                     <  1:2597504/71‑1 (MQ=255)
              cAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATAtgtg                                                     <  1:1624919/71‑1 (MQ=255)
                    gCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:166099/71‑1 (MQ=255)
                    gCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:1638561/71‑1 (MQ=255)
                    gCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATgg                                                <  1:2823866/70‑1 (MQ=255)
                    gCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATgg                                                <  1:1706529/70‑1 (MQ=255)
                    gCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATgg                                                <  1:1213691/70‑1 (MQ=255)
                    gCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATgg                                                <  1:83914/70‑1 (MQ=255)
                         gccgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAg                                          >  1:1577670/1‑71 (MQ=255)
                         gccgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAg                                          >  1:1653245/1‑71 (MQ=255)
                         gccgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAg                                          >  1:585641/1‑71 (MQ=255)
                         gccgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAg                                          >  1:2338051/1‑70 (MQ=255)
                          ccgccCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAGcc                                        <  1:363223/71‑1 (MQ=255)
                                agagAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:2482132/59‑1 (MQ=255)
                                agagAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:2430397/59‑1 (MQ=255)
                                agagAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:3432996/59‑1 (MQ=255)
                                agagAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:3666899/59‑1 (MQ=255)
                                agagAATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:521574/59‑1 (MQ=255)
                                 gagaATAACTGCCCACACTGCCATGTCCTCCGCGATC‑‑aaaaa                                                              >  1:275616/1‑41 (MQ=255)
                                 gagaATAACTGCCCACACTGCCATGTCATCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATaa                                   >  1:2971531/1‑70 (MQ=255)
                                  agaATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATAtgtg                                                     <  1:2718460/51‑1 (MQ=255)
                                   gaATAACTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATAtgtg                                                     <  1:2060866/50‑1 (MQ=255)
                                        aCTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg                                >  1:3320128/1‑66 (MQ=255)
                                        aCTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg                                >  1:1017638/1‑66 (MQ=255)
                                        aCTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg                                >  1:1920622/1‑66 (MQ=255)
                                        aCTGCCCACACTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGcc                            >  1:2940566/1‑70 (MQ=255)
                                               acacTTCCATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGt                   <  1:612617/71‑1 (MQ=255)
                                               acacTGCCATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGt                   <  1:531338/71‑1 (MQ=255)
                                                  cTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATa                                           >  1:1144199/1‑45 (MQ=255)
                                                  cTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTaa                 >  1:208177/1‑71 (MQ=255)
                                                  cTGCCATGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTaa                 >  1:736545/1‑71 (MQ=255)
                                                     ccATGTCCTCCGCGATCAAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAAt                <  1:2031746/70‑1 (MQ=255)
                                                     ccATGTCCTCCGCGATCAAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAAt                <  1:446483/70‑1 (MQ=255)
                                                     ccATGTCCTCCGCGATCAAAAACCGATATGTGAATggg                                               <  1:183442/38‑1 (MQ=38)
                                                      cATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCt            >  1:3441373/1‑71 (MQ=255)
                                                        tGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAg                         <  1:842599/57‑1 (MQ=255)
                                                        tGTCCTCCGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAg                         <  1:1238095/57‑1 (MQ=255)
                                                          tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt           >  1:1002051/1‑69 (MQ=255)
                                                          tcctccGCGATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCtt           >  1:1148462/1‑69 (MQ=255)
                                                               cgcgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAg                                >  1:2268929/1‑43 (MQ=255)
                                                                 cgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGt                   <  1:3656260/54‑1 (MQ=39)
                                                                 cgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAatat  <  1:1129696/71‑1 (MQ=255)
                                                                 cgATCAAAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAatat  <  1:2729006/71‑1 (MQ=255)
                                                                      |                                                                  
CACTTTCATGTTCCCAACTGGCATAGCCGCCCAGAGAATAACTGCCCACACTGCCATGTCCTCCGCGATC‑‑AAAACCGATATGTGAATGGGAATAGCCCATAAAAGCGCCTAGCGTGGTAATTCCTTCAGGAATAT  >  W3110S.gb/2777350‑2777484

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: