Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,185,304 T→C 6.6% T95T (ACT→ACC gspA → general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,185,3040TC 6.6% 226.2 / 5.5 76T95T (ACT→ACCgspAgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base C (3/2);  total (37/39)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

GTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTTGATG  >  W3110S.gb/4185243‑4185368
                                                             |                                                                
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACt                                                                  >  1:1064567/1‑62 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCg                                                          >  1:2575097/1‑70 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:3028951/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:1794089/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:1925885/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:2190158/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:2691581/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:2811358/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:2947047/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:3003260/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:943621/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:3280343/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:3329947/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:3340749/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:6034/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:655518/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:757666/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         >  1:1300117/1‑71 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:1181397/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:1158538/71‑1 (MQ=255)
gTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:1209058/71‑1 (MQ=255)
 tGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:1184003/70‑1 (MQ=255)
 tGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:2188894/70‑1 (MQ=255)
 tGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:3079716/70‑1 (MQ=255)
 tGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:399995/70‑1 (MQ=255)
     aGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACCGATAGTc                                                           <  1:1510205/64‑1 (MQ=255)
     aGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACCGATAGTc                                                           <  1:1040358/64‑1 (MQ=255)
            ggACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttat                                              >  1:2888046/1‑70 (MQ=255)
                     atgagatgGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGaa                                    <  1:1144589/71‑1 (MQ=255)
                      tgagatgGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGaa                                    <  1:2630893/70‑1 (MQ=255)
                        agatgGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGc                                                         <  1:666908/47‑1 (MQ=255)
                                    aaTTACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGaa                                    <  1:2969975/56‑1 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttatta                                            >  1:3146117/1‑45 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttatta                                            >  1:1146596/1‑45 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttatta                                            >  1:765875/1‑45 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttatta                                            >  1:1318116/1‑45 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttatta                                            >  1:1494479/1‑45 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCttatta                                            >  1:1628153/1‑45 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGc                               >  1:3159849/1‑58 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:2435447/1‑67 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:2432284/1‑67 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:1068894/1‑67 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:34076/1‑67 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:3100257/1‑67 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:2713309/1‑67 (MQ=255)
                                       tACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:214729/1‑67 (MQ=255)
                                       tACTGCATGCGATTCATCGAACTGATAGTCGCTGccc                                                    >  1:694839/1‑37 (MQ=255)
                                       tACAGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGaa                      >  1:1992497/1‑67 (MQ=255)
                                        aCTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCAcc                 >  1:2552874/1‑71 (MQ=255)
                                            cATGCGCTTCATCGAACCGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCAccgccg             >  1:2830555/1‑71 (MQ=255)
                                            cATGCGCTTCATCGAACCGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCAccgccg             >  1:1987863/1‑71 (MQ=255)
                                            cATGCGCTTCATCGAACCGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCAccgccg             >  1:1205741/1‑71 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCTCAt          <  1:109216/70‑1 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCAt          <  1:2158065/70‑1 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCAt          <  1:1511111/70‑1 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCAt          <  1:1240166/70‑1 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCAt          <  1:2587399/70‑1 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCAt          <  1:1105534/70‑1 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATa         >  1:732240/1‑71 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATa         >  1:2381697/1‑71 (MQ=255)
                                                cgcTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATa         >  1:134415/1‑71 (MQ=255)
                                                      aTCGTACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:1144196/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCTCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:3314849/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCTGCATACTtga    <  1:854763/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:1009426/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:141925/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:2600388/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:1079097/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:381944/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:1617360/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:2743870/70‑1 (MQ=255)
                                                      aTCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:1147963/70‑1 (MQ=255)
                                                       tCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtgatg  <  1:369442/71‑1 (MQ=255)
                                                       tCGAACTGATAGTCGCTGCCCTCTTATTATTACTGAAAATGCTCACCTGAATCACCGCCGCATACTtgatg  <  1:2976153/71‑1 (MQ=255)
                                                        cGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:547295/68‑1 (MQ=255)
                                                         gAACTGATAGTCGCTGCCCTCTTATTATTATTGaaaa                                  >  1:2641082/1‑37 (MQ=255)
                                                          aaCTGATAGTAGCTGCCCTATTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTtga    <  1:2625234/66‑1 (MQ=255)
                                                             |                                                                
GTGACAGCGCCTGGACGCGAGATGAGATGGCGGGGCAATTACTGCATGCGCTTCATCGAACTGATAGTCGCTGCCCTCTTATTATTATTGAAAATGCTCACCTGAATCACCGCCGCATACTTGATG  >  W3110S.gb/4185243‑4185368

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: