Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,546,874 G→A 7.8% G84G (GGC→GGT cysW ← sulfate/thiosulfate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,546,8740GA 7.8% 226.8 / 6.7 77G84G (GGC→GGTcysWsulfate/thiosulfate transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (41/30);  new base A (3/3);  total (44/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATG  >  W3110S.gb/2546805‑2546943
                                                                     |                                                                     
gTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGc                                                                      >  1:2335265/1‑71 (MQ=255)
gTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGc                                                                      >  1:3663788/1‑71 (MQ=255)
  ccAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGcc                                                                     >  1:3573594/1‑70 (MQ=255)
     gTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGaa                                                                  >  1:1819887/1‑70 (MQ=255)
     gTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGaa                                                                  >  1:2404122/1‑70 (MQ=255)
             agcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGtt                                                         >  1:3761254/1‑71 (MQ=255)
             agcagTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGtt                                                         >  1:2713199/1‑71 (MQ=255)
                agTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTAc                                                      <  1:2024923/71‑1 (MQ=255)
                   aaCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAaca                                                                >  1:875751/1‑58 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   >  1:889416/1‑70 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   >  1:1240075/1‑70 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   >  1:3320065/1‑70 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   >  1:890166/1‑70 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   <  1:3391524/70‑1 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   >  1:3472623/1‑70 (MQ=255)
                    aCTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCgg                                                   <  1:25029/70‑1 (MQ=255)
                         cGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGc                                             <  1:1527678/71‑1 (MQ=255)
                              aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAgg                                                           <  1:3442071/52‑1 (MQ=255)
                              aGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAAt                                          >  1:679740/1‑69 (MQ=255)
                               gggAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAgg                                                           <  1:1086771/51‑1 (MQ=255)
                                  aaGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAgcgc                                    >  1:124116/1‑71 (MQ=255)
                                   aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGcgcg                                   <  1:15275/71‑1 (MQ=255)
                                   aGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACAGGTACGGCAATCAGcgcg                                   <  1:3567024/71‑1 (MQ=255)
                                    gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTAcc                                                     >  1:3553064/1‑52 (MQ=255)
                                    gTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTAcc                                                     >  1:1542768/1‑52 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCa                                                             <  1:3649665/43‑1 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 <  1:3669467/71‑1 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGa                                  <  1:1039991/70‑1 (MQ=255)
                                     ttAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGa                                  <  1:2629025/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)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGaa                                                                  >  1:3494911/1‑37 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:780045/1‑70 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:3848421/1‑70 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:1390742/1‑70 (MQ=255)
                                      tAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:2570307/1‑70 (MQ=255)
                                       aaaGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatca                               >  1:77770/1‑71 (MQ=255)
                                         aGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcatc                             <  1:1560657/71‑1 (MQ=255)
                                         aGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGatcatc                             <  1:1733100/71‑1 (MQ=255)
                                             cgTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCg                         <  1:3124180/71‑1 (MQ=255)
                                                 aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                     >  1:1097660/1‑71 (MQ=255)
                                                 aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                     >  1:888332/1‑71 (MQ=255)
                                                 aCCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAg                     >  1:1597718/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)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATc                                         >  1:1131437/1‑48 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATc                                         >  1:849851/1‑48 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:2196249/1‑56 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                  <  1:3543267/71‑1 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                  <  1:2188521/71‑1 (MQ=255)
                                                    agccagGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCa                  <  1:3660935/71‑1 (MQ=255)
                                                     gccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAg                 <  1:3583587/71‑1 (MQ=255)
                                                      ccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCa                                        >  1:2594532/1‑47 (MQ=255)
                                                      ccaggccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGa                <  1:655220/71‑1 (MQ=255)
                                                           ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:3475879/1‑49 (MQ=255)
                                                           ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:806121/1‑49 (MQ=255)
                                                           ccagCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGat                                 >  1:2097196/1‑49 (MQ=255)
                                                            cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCg           >  1:1312601/1‑70 (MQ=255)
                                                            cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCg           >  1:1680972/1‑70 (MQ=255)
                                                            cagcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGt          >  1:442948/1‑71 (MQ=255)
                                                             agcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTg         <  1:1641828/71‑1 (MQ=255)
                                                             agcagAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTg         <  1:3165363/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:2206729/1‑71 (MQ=255)
                                                                    tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg  >  1:1554308/1‑71 (MQ=255)
                                                                    tGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATg  <  1:1007031/71‑1 (MQ=255)
                                                                     |                                                                     
GTCCAGTAGCGTCAGCAGTAACTGGCGTCCAGGGAAGTTAAAGCGCGTCACCAGCCAGGCCAGCAGAATGCCGAACACCAGGTTTACCGGTACGGCAATCAGCGCGATCATCACCGTCAGCCAGATGGCGTGCAGCATG  >  W3110S.gb/2546805‑2546943

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: