Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,345,374 C→T 24.3% G1156S (GGT→AGT)  yfaL ← adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,345,3740CT24.3% 69.0 / 19.1 37G1156S (GGT→AGT) yfaLadhesin
Reads supporting (aligned to +/- strand):  ref base C (15/13);  new base T (6/3);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTTGGTATGACATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCA  >  W3110S.gb/2345305‑2345442
                                                                     |                                                                    
aTGTTGGTATGACATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCAccc                                                                     >  1:2378279/1‑71 (MQ=255)
aTGTTGGTATGACATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCAcc                                                                      <  1:1092769/70‑1 (MQ=255)
   ttGGTATGACATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGc                                                                  >  1:626795/1‑71 (MQ=255)
         tGACATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACTCTGCGATtgt                                                            <  1:1752630/71‑1 (MQ=255)
           aCATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCAcc                                                                      >  1:2576155/1‑59 (MQ=255)
           aCATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCAcc                                                                      >  1:1094963/1‑59 (MQ=255)
           aCATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGa                                                          >  1:2478508/1‑71 (MQ=255)
           aCATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGa                                                          <  1:1305433/71‑1 (MQ=255)
           aCATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGa                                                          >  1:2319025/1‑71 (MQ=255)
            cATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGac                                                         <  1:668911/71‑1 (MQ=255)
             aTGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGa                                                          >  1:2552381/1‑69 (MQ=255)
                aaCAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACTCTGCGATTGTGACAcgc                                                     >  1:1208650/1‑71 (MQ=255)
                  cAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACTCTGCGATTGTGCCAcgcg                                                    >  1:2181742/1‑70 (MQ=255)
                   aGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACTCTGCGATTGTGACAcgcgcg                                                  <  1:1019541/71‑1 (MQ=255)
                                 cGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTa                                          >  1:1908708/1‑65 (MQ=255)
                                 cGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAgcggc                                     <  1:2772578/70‑1 (MQ=255)
                                      tgtAAACCCAGCCGCGTCTGAATATCATCACTCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTaa                                <  1:662788/70‑1 (MQ=255)
                                      tgtAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTa                                 <  1:837380/69‑1 (MQ=255)
                                          aaCCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatca                           >  1:2165466/1‑71 (MQ=255)
                                          aaCCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatca                           >  1:1871424/1‑71 (MQ=255)
                                          aaCCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatca                           <  1:2477416/71‑1 (MQ=255)
                                           aCCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatc                            <  1:1293929/69‑1 (MQ=255)
                                             ccAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATcc                        >  1:1887032/1‑71 (MQ=255)
                                             ccAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATcc                        >  1:1886221/1‑71 (MQ=255)
                                             ccAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATcc                        >  1:2391376/1‑71 (MQ=255)
                                                gCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatc                            <  1:2544076/64‑1 (MQ=255)
                                                    cgTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAat                             <  1:748557/59‑1 (MQ=255)
                                                    cgTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAat                             <  1:802224/59‑1 (MQ=255)
                                                      tCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACg               >  1:2591025/1‑71 (MQ=255)
                                                       cTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCctgc                     >  1:1137078/1‑64 (MQ=255)
                                                          aaTATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatca                           >  1:1044313/1‑55 (MQ=255)
                                                              tcatcaCTCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGccc            >  1:1157036/1‑66 (MQ=255)
                                                              tcatcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatc                            <  1:1870173/50‑1 (MQ=255)
                                                              tcatcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatc                            <  1:1826269/50‑1 (MQ=255)
                                                              tcatcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCctgctg                   <  1:2522099/59‑1 (MQ=255)
                                                                atcaCTCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatcatc                         >  1:2556034/1‑51 (MQ=255)
                                                                atcaCTCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatcatc                         >  1:743069/1‑51 (MQ=255)
                                                                atcaCTCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatcatc                         >  1:965554/1‑51 (MQ=255)
                                                                atcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatcat                          >  1:2647684/1‑50 (MQ=255)
                                                                 tcaCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAatca                           <  1:2694252/48‑1 (MQ=255)
                                                                   aCCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCa  >  1:1051115/1‑71 (MQ=255)
                                                                     |                                                                    
ATGTTGGTATGACATGAACAGCGGTACGCCATTCGCTGTGTAAACCCAGCCGCGTCTGAATATCATCACCCTGCGATTGTGACACGCGCGCACGGTTAGCGGCGGTAAAATCATCCTGCTGCACGCCCTGATAAATCA  >  W3110S.gb/2345305‑2345442

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: