Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,893,190 A→G 5.2% S113S (TCA→TCG ygcP → predicted anti‑terminator regulatory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,893,1900AG 5.2% 210.6 / 2.6 77S113S (TCA→TCGygcPpredicted anti‑terminator regulatory protein
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base G (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCT  >  W3110S.gb/2893127‑2893260
                                                               |                                                                      
aaTGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTAATTTCg                                                                  <  1:4388546/69‑1 (MQ=255)
  tGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGttt                                                               <  1:3692789/71‑1 (MQ=255)
     tGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTca                                                             >  1:4126442/1‑70 (MQ=255)
      gAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAAt                                                                      >  1:1595618/1‑60 (MQ=255)
      gAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTcac                                                            >  1:3590478/1‑70 (MQ=255)
      gAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTc                                                              <  1:1105622/68‑1 (MQ=255)
           cGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACAtt                                                      >  1:4255328/1‑71 (MQ=255)
            ggCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTg                                                     <  1:1355817/71‑1 (MQ=255)
            ggCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTg                                                     <  1:528554/71‑1 (MQ=255)
            ggCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTg                                                     <  1:4093589/71‑1 (MQ=255)
             gCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGa                                                    >  1:380587/1‑71 (MQ=255)
             gCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGa                                                    >  1:2448097/1‑71 (MQ=255)
                agagCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATaa                                                 <  1:1973996/71‑1 (MQ=255)
                       aaGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGtt                                          <  1:3925979/71‑1 (MQ=255)
                        aGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTg                                         >  1:2893280/1‑71 (MQ=255)
                         ggTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCacaa                                                          >  1:3259596/1‑53 (MQ=255)
                           ttttttCTGTATTCATCGCCTGTTTATTGTTGATTCAAtt                                                                     <  1:3371729/40‑1 (MQ=255)
                           ttttttCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTca                                                             <  1:273237/48‑1 (MQ=255)
                            tttttCTTTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCaa                                             >  1:2089105/1‑63 (MQ=255)
                            tttttCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTgcgc                                      >  1:2936460/1‑70 (MQ=255)
                              tttCTGTATTCATCGCGTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCaa                                    <  1:2208324/70‑1 (MQ=255)
                                  tttATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCa                                              <  1:3779451/54‑1 (MQ=255)
                                    tATTCATCGCCTGTTTATTGTTGATTCAATTTCGt                                                                 <  1:1992606/35‑1 (MQ=255)
                                    tATTCATCGCCTGTTTATTGTTGATTCAATTTCGt                                                                 <  1:3213452/35‑1 (MQ=255)
                                    tATTCATCGCCTGTTTATTGTTGATTCAATTTCGt                                                                 <  1:3208588/35‑1 (MQ=255)
                                    tATTCATCGCCTGTTTATTGTTGATTCAATTTCGt                                                                 <  1:2396270/35‑1 (MQ=255)
                                    tATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGa                                                    <  1:1640906/48‑1 (MQ=255)
                                    tATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAAt                              <  1:3393626/70‑1 (MQ=255)
                                     aTTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATcc                            >  1:125441/1‑71 (MQ=255)
                                       tCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCa                                              <  1:1739758/51‑1 (MQ=255)
                                       tCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCa                                              <  1:3681778/51‑1 (MQ=255)
                                        cATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGAt                                                   <  1:2244729/45‑1 (MQ=255)
                                        cATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGAt                                                   >  1:3676156/1‑45 (MQ=255)
                                        cATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTCCGCAATCGa                                >  1:2401080/1‑64 (MQ=255)
                                         aTCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGa                         <  1:860687/70‑1 (MQ=255)
                                         aTCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGa                         <  1:1453643/70‑1 (MQ=255)
                                          tCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCaa                                             <  1:3980659/49‑1 (MQ=255)
                                          tCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCaa                                             <  1:890327/49‑1 (MQ=255)
                                          tCGCCTGTTTATTGTTGATTCAATTCCGTTTCACAACATTg                                                     <  1:2731821/41‑1 (MQ=255)
                                           cGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg                      >  1:2235593/1‑71 (MQ=255)
                                            gCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg                      <  1:132501/70‑1 (MQ=255)
                                            gCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGt                     <  1:3614814/71‑1 (MQ=255)
                                              cTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCTCAATCGAATCCGGATTGTAt                   <  1:1292301/71‑1 (MQ=255)
                                               tGTTTATTGTTGATTCGATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg                      <  1:58777/67‑1 (MQ=255)
                                               tGTTTATTGTTGATTCGATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg                      <  1:4422084/67‑1 (MQ=255)
                                                gTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTAtt                  >  1:2797780/1‑70 (MQ=255)
                                                 tttATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCa                                              <  1:3184333/41‑1 (MQ=255)
                                                 tttATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCa                                              <  1:744481/41‑1 (MQ=255)
                                                 tttATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGt                                           <  1:2691241/44‑1 (MQ=255)
                                                 tttATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCgg                          <  1:1088823/61‑1 (MQ=255)
                                                 tttATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGa                         <  1:4407446/62‑1 (MQ=255)
                                                 tttATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTAtt                  >  1:1000982/1‑69 (MQ=255)
                                                 tttATTGTTATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTAtt                  >  1:1228342/1‑68 (MQ=255)
                                                   tATTGTTGATTCGATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGt                     >  1:1265457/1‑64 (MQ=255)
                                                   tATTGTTGATTCGATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGt                     >  1:974383/1‑64 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGtt                                          >  1:3514640/1‑43 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGGAATCGAATCCGGAt                        >  1:462649/1‑61 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAAt                                   >  1:2031044/1‑50 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGa                         >  1:237710/1‑60 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTg                      <  1:3105194/63‑1 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTAtt                  >  1:1877125/1‑67 (MQ=255)
                                                   tATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTgag               <  1:2559812/70‑1 (MQ=255)
                                                     ttgttgATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTgag               <  1:717822/68‑1 (MQ=255)
                                                     ttgttgATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGAt             <  1:3477785/70‑1 (MQ=255)
                                                         tgATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGc        >  1:1677367/1‑71 (MQ=255)
                                                          gATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATc                                  >  1:2909057/1‑44 (MQ=255)
                                                          gATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTa                    >  1:1993675/1‑58 (MQ=255)
                                                          gATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTAt                   >  1:977766/1‑59 (MQ=255)
                                                          gATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCttc        >  1:1021062/1‑68 (MQ=255)
                                                          gATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGc        <  1:1206282/70‑1 (MQ=255)
                                                            ttCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTgaga              >  1:1753451/1‑62 (MQ=255)
                                                            ttCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAg     >  1:2337562/1‑71 (MQ=255)
                                                            ttCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAg     >  1:1251989/1‑71 (MQ=255)
                                                             tCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAgg    >  1:1512748/1‑71 (MQ=255)
                                                              cAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGaa                               >  1:3562163/1‑43 (MQ=255)
                                                              cAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGAt             >  1:1059566/1‑61 (MQ=255)
                                                              cAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGc   >  1:425496/1‑71 (MQ=255)
                                                               aaTTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGaa                               >  1:4151392/1‑42 (MQ=255)
                                                               aaTTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCt  <  1:3891586/71‑1 (MQ=255)
                                                               aaTTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCt  <  1:2646962/71‑1 (MQ=255)
                                                               |                                                                      
AATGTTGAAAGCGGCAAGAGCGGAAGGTTTTTTCTGTATTCATCGCCTGTTTATTGTTGATTCAATTTCGTTTCACAACATTGATAAGCAAGTTGCGCAATCGAATCCGGATTGTATTGAGATCCTGCCAGGCT  >  W3110S.gb/2893127‑2893260

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