Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,989,358 T→A 5.9% S99R (AGT→AGA yecI → predicted ferritin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,989,3580TA 5.9% 188.1 / 2.9 67S99R (AGT→AGAyecIpredicted ferritin‑like protein
Reads supporting (aligned to +/- strand):  ref base T (36/27);  new base A (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATC  >  W3110S.gb/1989289‑1989428
                                                                     |                                                                      
cGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTa                                                                       >  1:4070895/1‑71 (MQ=255)
  gTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACg                                                                     >  1:3106574/1‑71 (MQ=255)
  gTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACg                                                                     >  1:3101792/1‑71 (MQ=255)
      aaaaCTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTgg                                                                 >  1:1207154/1‑71 (MQ=255)
         aCTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTgg                                                                 >  1:1005486/1‑68 (MQ=255)
         aCTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcac                                                              <  1:4037170/71‑1 (MQ=255)
          cTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGca                                                               >  1:1724183/1‑69 (MQ=255)
                 ctctGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGcc                                                      >  1:1411295/1‑71 (MQ=255)
                    tGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGAt                                                   >  1:1484627/1‑71 (MQ=255)
                     ggAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcac                                                              <  1:2299111/59‑1 (MQ=255)
                      gaagaaCTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATg                                                  <  1:2397629/70‑1 (MQ=255)
                       aagaaCTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATg                                                  >  1:3668621/1‑69 (MQ=255)
                           aCTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGa                                            <  1:975676/71‑1 (MQ=255)
                            cTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGaa                                           <  1:3583339/71‑1 (MQ=255)
                                tCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGaa                                       >  1:2712975/1‑71 (MQ=255)
                                 ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACg                                                                     >  1:3045227/1‑40 (MQ=255)
                                 ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcaca                                                             >  1:2591661/1‑48 (MQ=255)
                                 ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGaa                                       <  1:481235/70‑1 (MQ=255)
                                 ccAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAAc                                      >  1:1975969/1‑71 (MQ=255)
                                   aaaaaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGa                                                 <  1:1924712/58‑1 (MQ=255)
                                   aaaaaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                                     >  1:2311721/1‑70 (MQ=255)
                                    aaaaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGcac                                                              <  1:602557/44‑1 (MQ=255)
                                      aaaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAat                                 <  1:287349/71‑1 (MQ=255)
                                       aaCGATGGAAGAATACGAGCAACGTTCTAGTAGGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAat                                 <  1:564947/70‑1 (MQ=255)
                                       aaCGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatg                                <  1:2378276/71‑1 (MQ=255)
                                         cGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgat                              <  1:788572/71‑1 (MQ=255)
                                          gATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                             <  1:3424932/71‑1 (MQ=255)
                                          gATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                             <  1:1526557/71‑1 (MQ=255)
                                            tGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                                     >  1:2580830/1‑61 (MQ=255)
                                            tGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                                     >  1:2575722/1‑61 (MQ=255)
                                            tGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                                     >  1:67219/1‑61 (MQ=255)
                                               aagaaTACGATCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGaa                                       >  1:2597295/1‑56 (MQ=255)
                                                agaaTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                             >  1:1640498/1‑65 (MQ=255)
                                                  aaTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAAcc                     <  1:2824126/71‑1 (MQ=255)
                                                   aTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAAcc                     <  1:841721/70‑1 (MQ=255)
                                                   aTACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCg                    <  1:1975789/71‑1 (MQ=255)
                                                    tACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAg                                               >  1:2863778/1‑43 (MQ=255)
                                                    tACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt                   >  1:638149/1‑71 (MQ=255)
                                                     aCGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCg                                                     >  1:2374921/1‑36 (MQ=255)
                                                     aCGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt                   >  1:1190295/1‑70 (MQ=255)
                                                      cGAGCAACGTTCTAGTACGTTGTCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTc                  >  1:582570/1‑70 (MQ=255)
                                                      cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAg                                         >  1:3425003/1‑47 (MQ=255)
                                                      cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAg                                         >  1:1128228/1‑47 (MQ=255)
                                                      cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAAc                      >  1:3879954/1‑66 (MQ=255)
                                                      cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTc                  >  1:2209999/1‑70 (MQ=255)
                                                      cGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCa                 <  1:1887058/71‑1 (MQ=255)
                                                         gCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAAtt              <  1:341793/71‑1 (MQ=255)
                                                         ccAACGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt                   <  1:1604204/65‑1 (MQ=39)
                                                         ccAACGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGt                   <  1:3826346/65‑1 (MQ=39)
                                                          cAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAAtgatga                            >  1:100485/1‑56 (MQ=255)
                                                           aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGa                                            <  1:1008349/39‑1 (MQ=255)
                                                           aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTg                                    <  1:3066032/47‑1 (MQ=255)
                                                           aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAAttt             <  1:4075592/70‑1 (MQ=255)
                                                           aaCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTc            <  1:806378/71‑1 (MQ=255)
                                                           aaCGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCa                 >  1:1578734/1‑66 (MQ=255)
                                                           aaCGTACTAGAACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCa                 >  1:3293269/1‑66 (MQ=255)
                                                            aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACt                                     >  1:3815623/1‑45 (MQ=255)
                                                            aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAatgatg                             <  1:2504854/53‑1 (MQ=255)
                                                            aCGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTc            <  1:451451/70‑1 (MQ=255)
                                                             cGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCt          >  1:3366901/1‑71 (MQ=255)
                                                              gTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTc                  <  1:1996702/62‑1 (MQ=255)
                                                               ttCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCaa                >  1:317425/1‑63 (MQ=255)
                                                                 cTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAAttt             >  1:3432893/1‑64 (MQ=255)
                                                                 cTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTgcgc      >  1:4208888/1‑71 (MQ=255)
                                                                 cTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTgcg       >  1:583074/1‑70 (MQ=255)
                                                                    gTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGAt   <  1:1199559/71‑1 (MQ=255)
                                                                    gTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGAt   >  1:3250660/1‑71 (MQ=255)
                                                                     tACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATc  <  1:45133/71‑1 (MQ=255)
                                                                     tACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATc  <  1:3540602/71‑1 (MQ=255)
                                                                     |                                                                      
CGGTGAAAAACTGAACTCTCTGGAAGAACTGTTCCAAAAAACGATGGAAGAATACGAGCAACGTTCTAGTACGTTGGCACAGTTAGCCGATGAAGCGAAAGAACTGAATGATGATTCAACCGTCAATTTCCTGCGCGATC  >  W3110S.gb/1989289‑1989428

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: