Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,457,604 A→G 8.7% L60P (CTG→CCG)  yfcS ← predicted periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,457,6040AG 8.7% 185.3 / 8.9 69L60P (CTG→CCG) yfcSpredicted periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base G (4/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CACTGGCGGCAGCACAGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCA  >  W3110S.gb/2457534‑2457673
                                                                      |                                                                     
cACTGGCGGCAGCACAGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCa                                                                       <  1:2123516/71‑1 (MQ=255)
cACTGGCGGCAGCACAGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCa                                                                       <  1:1474247/71‑1 (MQ=255)
       gccagcaCAGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAgg                                                                <  1:238563/69‑1 (MQ=255)
        gcagcaCAGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAggg                                                               <  1:1719883/71‑1 (MQ=255)
             acaGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGa                                                              >  1:3991330/1‑67 (MQ=255)
             acaGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGttt                                                          >  1:2442708/1‑71 (MQ=255)
             acaGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGttt                                                          >  1:1682115/1‑71 (MQ=255)
              caGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCGGATAAGGGAGTTTc                                                         >  1:382590/1‑71 (MQ=255)
              caGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTc                                                         >  1:2544145/1‑71 (MQ=255)
              caGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCAAGCTTTGAGCCAGATAAGGGAGTTTc                                                         >  1:3520465/1‑71 (MQ=255)
               aGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATaa                                                                  >  1:2684230/1‑61 (MQ=255)
               aGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCg                                                        <  1:2605059/71‑1 (MQ=255)
                gTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCAAGCTTTGAGCCAGATAAGGGAGTTTCg                                                        >  1:104816/1‑70 (MQ=255)
                   aaCGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCGGATAAGGGAGTTTCGGATc                                                    >  1:3595520/1‑71 (MQ=255)
                     cGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGt                                                  <  1:1178704/71‑1 (MQ=255)
                        cGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATa                                                                   <  1:1214452/51‑1 (MQ=255)
                        cGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGttat                                               <  1:3982399/71‑1 (MQ=255)
                        cGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGttat                                               <  1:2451686/71‑1 (MQ=255)
                        cGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGttat                                               <  1:3464235/71‑1 (MQ=255)
                            gTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATa                                                                   <  1:765709/47‑1 (MQ=255)
                              gATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAg                                                                 >  1:3598079/1‑47 (MQ=255)
                              gATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAg                                                                 >  1:3552915/1‑47 (MQ=255)
                               aTCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAg                                        >  1:886804/1‑71 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCGGATAAGGGAGTTTCGGATCGTTATTACGCAGCg                                      >  1:3756101/1‑71 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCGGATAAGGGAGTTTCGGATCGTTATTACGCAGCg                                      >  1:1533951/1‑71 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAg                                                                      >  1:3381108/1‑39 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAg                                                             >  1:2376523/1‑48 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGc                                       >  1:1208754/1‑70 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCg                                      >  1:1720602/1‑71 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCg                                      >  1:73355/1‑71 (MQ=255)
                                 cTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCg                                      >  1:2547180/1‑71 (MQ=255)
                                    gTTACCCTTTTCATCTTCAATCCAGCTTTGAGGCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGt                                     <  1:2527957/69‑1 (MQ=255)
                                    gTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGt                                     <  1:1479049/69‑1 (MQ=255)
                                    gTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGt                                     <  1:1102476/69‑1 (MQ=255)
                                    gTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGt                                     <  1:2117792/69‑1 (MQ=255)
                                    gTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGt                                     <  1:1653838/69‑1 (MQ=255)
                                    gTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGt                                     <  1:3903506/69‑1 (MQ=255)
                                     ttACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCAc                                  >  1:10247/1‑71 (MQ=255)
                                       aCCCTTTTCATCTTCTATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGt                                                  <  1:2879273/53‑1 (MQ=255)
                                         ccTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGTGTTTCGGATCGTTATTACGCAGCGTCACGCtg                              <  1:3759680/71‑1 (MQ=255)
                                         ccTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCt                               <  1:3707069/70‑1 (MQ=255)
                                           ttttCATCTTCAATCCATCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCt                               <  1:2583691/68‑1 (MQ=255)
                                           ttttCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGttat                                               <  1:528936/52‑1 (MQ=255)
                                           ttttCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGttat                                               <  1:3463114/52‑1 (MQ=255)
                                           ttttCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCt                               <  1:3730384/68‑1 (MQ=255)
                                            tttCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGt                                                  <  1:3362797/48‑1 (MQ=255)
                                             ttCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCtgatt                           >  1:3455661/1‑70 (MQ=255)
                                              tCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGtta                                                >  1:2256227/1‑48 (MQ=255)
                                                 tCTTCAATCCAGCTTTGAGCCGGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCt                               <  1:824521/62‑1 (MQ=255)
                                                 tCTTCAATCCAGCTTTGAGCCGGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCt                               <  1:2415878/62‑1 (MQ=255)
                                                 tCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCAc                                  <  1:2432808/59‑1 (MQ=255)
                                                  cTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTg                     >  1:1008629/1‑71 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGt                                                  <  1:336629/41‑1 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGt                                                  <  1:274629/41‑1 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCtgatt                           <  1:3225498/64‑1 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCtgatt                           <  1:3205487/64‑1 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCtgatt                           <  1:3970295/64‑1 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTgattgatt                       >  1:263500/1‑68 (MQ=255)
                                                   ttCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGt                    >  1:1378080/1‑71 (MQ=255)
                                                     cAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGAttt                      <  1:374539/67‑1 (MQ=255)
                                                          cAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTc             <  1:1760602/71‑1 (MQ=255)
                                                          cAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTc             <  1:1210191/71‑1 (MQ=255)
                                                            gccttGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCAt           >  1:1203376/4‑71 (MQ=255)
                                                            gCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCAt           >  1:3172550/1‑71 (MQ=255)
                                                            gCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCAt           >  1:1567158/1‑71 (MQ=255)
                                                               ttGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGc                                >  1:1475658/1‑47 (MQ=255)
                                                               ttGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCAt           >  1:1486187/1‑68 (MQ=255)
                                                               ttGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTa         >  1:201910/1‑70 (MQ=255)
                                                                     cAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCa  <  1:3221829/71‑1 (MQ=255)
                                                                     cAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCa  <  1:271423/71‑1 (MQ=255)
                                                                     cAGATAAGGGAGTTTCCGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCa  <  1:1643040/71‑1 (MQ=255)
                                                                      |                                                                     
CACTGGCGGCAGCACAGTTAACGGCGAGGTGATCTTGTTACCCTTTTCATCTTCAATCCAGCTTTGAGCCAGATAAGGGAGTTTCGGATCGTTATTACGCAGCGTCACGCTGATTGATTTGTCGCTTTCATTAAATATCA  >  W3110S.gb/2457534‑2457673

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: