Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,581,839 C→T 57.7% M222I (ATG→ATA ydeM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,581,8390CT57.7% 21.1 / 65.3 52M222I (ATG→ATAydeMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (11/11);  new base T (12/18);  total (23/29)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACT  >  W3110S.gb/1581772‑1581901
                                                                   |                                                              
aGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGAc                                                                >  1:1865283/1‑68 (MQ=255)
 gCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                               >  1:432020/1‑68 (MQ=255)
  cTTTCAAACTTACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGAtata                                                             <  1:555766/69‑1 (MQ=255)
  cTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGAtata                                                             <  1:1638177/69‑1 (MQ=255)
  cTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGAtata                                                             <  1:3434872/69‑1 (MQ=255)
                gATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGc                                               <  1:411956/69‑1 (MQ=255)
                gATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGc                                               <  1:1980679/69‑1 (MQ=255)
                 aTGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGcc                                              <  1:1286962/69‑1 (MQ=255)
                  tGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCg                                             >  1:562161/1‑69 (MQ=255)
                                  aTCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaa                              >  1:536281/1‑68 (MQ=255)
                                   tggTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                             >  1:1125080/3‑68 (MQ=255)
                                   tCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGt                                            >  1:595082/1‑53 (MQ=255)
                                   tCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                             >  1:2788278/1‑68 (MQ=255)
                                   tCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                             >  1:868099/1‑68 (MQ=255)
                                   tCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                             >  1:569338/1‑68 (MQ=255)
                                   tCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                             >  1:3101012/1‑68 (MQ=255)
                                   tCGTTTTTAACCCATTGCATAAAAATGGTTGATATAAACTTGCAATAAGCCGt                                            >  1:3305833/1‑53 (MQ=255)
                                    cGTTTTTAACCCATTGCATAAAAATGGTTTATATaa                                                            <  1:478234/36‑1 (MQ=38)
                                     gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                          >  1:2864350/1‑69 (MQ=255)
                                     gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                          >  1:457968/1‑69 (MQ=255)
                                          tAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAg                                                <  1:2184642/42‑1 (MQ=255)
                                          tAACCCATTGCATAAAAATGGTTGATATAAACTTGCCATAAg                                                <  1:614053/42‑1 (MQ=255)
                                           aaCCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                     <  1:1857149/68‑1 (MQ=255)
                                                aTTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGTCACAGAAAAATCAAt                       <  1:2108742/61‑1 (MQ=255)
                                                aTTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                       >  1:1735405/1‑61 (MQ=255)
                                                aTTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                       >  1:865857/1‑61 (MQ=255)
                                                aTTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                       >  1:2639365/1‑61 (MQ=255)
                                                  tGCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAt             <  1:527989/69‑1 (MQ=255)
                                                   gCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACgg                <  1:812526/65‑1 (MQ=255)
                                                   gCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACgg                <  1:169943/65‑1 (MQ=255)
                                                   gCATAAAAATGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACgg                <  1:291117/65‑1 (MQ=255)
                                                   gCATAAAAATGGTTGACATAAACTTGCCATAAGGCGTGGGAGGCACAGAAAAATCAATGATACGGAAtg            >  1:117652/1‑69 (MQ=255)
                                                   gCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaa                              <  1:2900501/51‑1 (MQ=255)
                                                   gCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtg            >  1:1117481/1‑69 (MQ=255)
                                                     aTAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                          <  1:2313229/37‑1 (MQ=255)
                                                       aaaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        <  1:3133232/69‑1 (MQ=255)
                                                       aaaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        <  1:2234040/69‑1 (MQ=255)
                                                       aaaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        <  1:3443955/69‑1 (MQ=255)
                                                        aaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATAc                  >  1:972022/1‑58 (MQ=255)
                                                        aaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgtg          <  1:3006663/66‑1 (MQ=255)
                                                        aaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        >  1:2214719/1‑68 (MQ=255)
                                                        aaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        >  1:13279/1‑68 (MQ=255)
                                                        aaaaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        >  1:1636295/1‑68 (MQ=255)
                                                        aaaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        >  1:195240/1‑68 (MQ=255)
                                                        aaaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt        >  1:1880384/1‑68 (MQ=255)
                                                         aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct      <  1:1338741/69‑1 (MQ=255)
                                                          aaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc     <  1:2308845/69‑1 (MQ=255)
                                                          aaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct      <  1:3555877/68‑1 (MQ=255)
                                                          aaTGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct      <  1:3548556/68‑1 (MQ=255)
                                                          aaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct      <  1:2214483/68‑1 (MQ=255)
                                                          aaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct      <  1:2073548/68‑1 (MQ=255)
                                                          aaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct      <  1:3149745/68‑1 (MQ=255)
                                                            tGGTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc   <  1:2502479/69‑1 (MQ=255)
                                                             ggTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt  <  1:159295/69‑1 (MQ=255)
                                                             ggTTGATATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt  <  1:164058/69‑1 (MQ=255)
                                                                tGATATAAACTTGCCATAAGCCGTGGGAGGCACAGaa                               <  1:1559042/37‑1 (MQ=255)
                                                                   |                                                              
AGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACT  >  W3110S.gb/1581772‑1581901

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: