Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,546,906 G→A 7.5% L74L (CTG→TTG)  cysW ← sulfate/thiosulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,546,9060GA 7.5% 239.8 / 5.6 80L74L (CTG→TTG) cysWsulfate/thiosulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (39/35);  new base A (3/3);  total (42/38)
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

AAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATC  >  W3110S.gb/2546839‑2546976
                                                                   |                                                                      
aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAgcgc                                                                     >  1:124116/1‑71 (MQ=255)
 aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGcgcg                                                                    <  1:15275/71‑1 (MQ=255)
 aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACAGGTACGGCAATCAGcgcg                                                                    <  1:3567024/71‑1 (MQ=255)
   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  <  1:3669467/71‑1 (MQ=255)
   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGa                                                                   <  1:2629025/70‑1 (MQ=255)
   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGa                                                                   <  1:1039991/70‑1 (MQ=255)
   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGa                                                                   <  1:761598/70‑1 (MQ=255)
   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATACCGAACACCAGGTTTACCGGTACGGCAATCAACGCGa                                                                   <  1:109577/70‑1 (MQ=255)
   ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATACCGAACACCAGGTTTACCGGTACGGCAATCAACGCGa                                                                   <  1:2475900/70‑1 (MQ=255)
    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:3848421/1‑70 (MQ=255)
    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:2570307/1‑70 (MQ=255)
    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:1390742/1‑70 (MQ=255)
    tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:780045/1‑70 (MQ=255)
     aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatca                                                                >  1:77770/1‑71 (MQ=255)
       aGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcatc                                                              <  1:1733100/71‑1 (MQ=255)
       aGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcatc                                                              <  1:1560657/71‑1 (MQ=255)
           cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCg                                                          <  1:3124180/71‑1 (MQ=255)
               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                      >  1:1597718/1‑71 (MQ=255)
               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                      >  1:888332/1‑71 (MQ=255)
               aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                                                      >  1:1097660/1‑71 (MQ=255)
               aCCAGCCAGGCCAGCAGAATGCCGAACACAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGc                                                     >  1:3014459/1‑71 (MQ=255)
                ccagccagGCCAGCAGAATACCGAACACCAGGTTTACCGGTACGGCAATCAAcgcgt                                                                   <  1:3757762/57‑2 (MQ=255)
                ccagccagGCCAGCAGAATACCGAACACCAGGTTTACCGGTACGGCAATCAACGCGATCATCACCGTCAGc                                                     >  1:2454369/1‑71 (MQ=255)
                ccagccagGCCAGCAGAATACCGAACACCAGGTTTACCGGTACGGCAATCAACGCGATCATCACCGTCAGc                                                     >  1:1741100/1‑71 (MQ=255)
                ccagccagGCCAGCAGAATACCGAACACCAGGTTTACCGGTACGGCAATCAACGCGATCATCACCGTCAGc                                                     >  1:2458646/1‑71 (MQ=255)
                 cagccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGc                                                     <  1:221611/70‑1 (MQ=255)
                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:2196249/1‑56 (MQ=255)
                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                   <  1:3660935/71‑1 (MQ=255)
                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                   <  1:3543267/71‑1 (MQ=255)
                  agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                   <  1:2188521/71‑1 (MQ=255)
                   gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAg                                                  <  1:3583587/71‑1 (MQ=255)
                    ccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                 <  1:655220/71‑1 (MQ=255)
                         ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:2097196/1‑49 (MQ=255)
                         ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:3475879/1‑49 (MQ=255)
                         ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                                                  >  1:806121/1‑49 (MQ=255)
                          cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCg                                            >  1:1680972/1‑70 (MQ=255)
                          cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCg                                            >  1:1312601/1‑70 (MQ=255)
                          cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGt                                           >  1:442948/1‑71 (MQ=255)
                           agcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTg                                          <  1:3165363/71‑1 (MQ=255)
                           agcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTg                                          <  1:1641828/71‑1 (MQ=255)
                            gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCg                                                          <  1:2712367/54‑1 (MQ=255)
                            gcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCg                                                          <  1:2270951/54‑1 (MQ=255)
                             cagAATGCCGAACACAAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgca                                        <  1:992280/71‑1 (MQ=255)
                              agAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                       >  1:1062696/1‑71 (MQ=255)
                               gAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcagc                                      <  1:1125755/71‑1 (MQ=255)
                                  tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                   >  1:536696/1‑71 (MQ=255)
                                  tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                   <  1:2612335/71‑1 (MQ=255)
                                  tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                   >  1:1554308/1‑71 (MQ=255)
                                  tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                   >  1:2206729/1‑71 (MQ=255)
                                  tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg                                   <  1:1007031/71‑1 (MQ=255)
                                    ccGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTc                                 <  1:3372586/71‑1 (MQ=255)
                                     cGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTcc                                >  1:447794/1‑71 (MQ=255)
                                        acacCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                             >  1:969530/1‑71 (MQ=255)
                                        acacCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGa                             <  1:2627423/71‑1 (MQ=255)
                                             aGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                   <  1:254697/44‑1 (MQ=255)
                                             aGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                        >  1:289965/1‑71 (MQ=255)
                                             aGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                        >  1:89773/1‑71 (MQ=255)
                                             aGGTTTACCGGTACGGCAATCAGCGCCATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                        >  1:2914345/1‑71 (MQ=255)
                                              ggTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                        >  1:2558348/1‑70 (MQ=255)
                                               gTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                                                   <  1:2250643/42‑1 (MQ=255)
                                               gTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAg                                                  <  1:1646531/43‑1 (MQ=255)
                                               gTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCa                      <  1:3385845/71‑1 (MQ=255)
                                                 ttACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                 >  1:3623669/1‑42 (MQ=255)
                                                 ttACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                                                 >  1:3157198/1‑42 (MQ=255)
                                                 ttACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGc                        >  1:670471/1‑67 (MQ=255)
                                                   aCCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGAt                   <  1:1468105/70‑1 (MQ=255)
                                                       gTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                       <  1:2591047/46‑1 (MQ=255)
                                                       gTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                       <  1:20795/46‑1 (MQ=255)
                                                       gTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                       <  1:2682057/46‑1 (MQ=255)
                                                       gTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTgcag                                       <  1:2633597/46‑1 (MQ=255)
                                                          cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCtt               >  1:836174/1‑66 (MQ=255)
                                                          cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa            >  1:1031383/1‑70 (MQ=255)
                                                          cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa            >  1:134538/1‑70 (MQ=255)
                                                          cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTaa            >  1:1440986/1‑70 (MQ=255)
                                                          cGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTa             >  1:3804913/1‑69 (MQ=255)
                                                             cAATCAGCGCGATCATCACCGTCAGCCAGATGGCGAGCAGCATGTCCGGAACGGCCAGATTCTGTaaaa          >  1:2925497/1‑69 (MQ=255)
                                                                  aGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCa    <  1:1880751/70‑1 (MQ=255)
                                                                  aGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt   >  1:30288/1‑71 (MQ=255)
                                                                  aGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt   >  1:1980460/1‑71 (MQ=255)
                                                                  aGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt   >  1:1139816/1‑71 (MQ=255)
                                                                  aGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCAt   >  1:34846/1‑71 (MQ=255)
                                                                   gcgcGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATc  <  1:1651129/71‑1 (MQ=255)
                                                                   |                                                                      
AAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATGTCCGGATCGGCCAGATTCTGTAAAACCGGCATC  >  W3110S.gb/2546839‑2546976

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: