Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,085,574 C→T 12.7% N71N (AAC→AAT metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,085,5740CT12.7% 139.6 / 10.2 55N71N (AAC→AATmetKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (16/32);  new base T (2/5);  total (18/37)
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.53e-01

TGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTCTG  >  W3110S.gb/3085510‑3085640
                                                                |                                                                  
tGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGtt                                                              >  1:107103/1‑71 (MQ=255)
  gTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTc                                                             <  1:363274/70‑1 (MQ=255)
     ttAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTcgcg                                                          >  1:71154/1‑70 (MQ=255)
       aGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGaa                                                        <  1:2634821/70‑1 (MQ=255)
         ttGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAAt                                                      <  1:1044153/70‑1 (MQ=255)
            gcggcgAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGc                                                  <  1:3319858/71‑1 (MQ=255)
              ggcgAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTa                                                <  1:1492418/71‑1 (MQ=255)
              ggcgAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTa                                                <  1:3557444/71‑1 (MQ=255)
               gcgAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTAt                                               <  1:1031229/71‑1 (MQ=255)
               gcgAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTAt                                               <  1:2545851/71‑1 (MQ=255)
                  aaaTCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTAtgtg                                            <  1:543389/71‑1 (MQ=255)
                    aTCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAATACCGTTCGCGAAATTGGCTATGTGCa                                          >  1:1668301/1‑71 (MQ=255)
                    aTCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAATACCGTTCGCGAAATTGGCTATGTGCa                                          >  1:1767542/1‑71 (MQ=255)
                    aTCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCa                                          >  1:641507/1‑71 (MQ=255)
                    aTCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCa                                          >  1:2366869/1‑71 (MQ=255)
                       accaccAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTc                                       >  1:3680630/1‑71 (MQ=255)
                         caccaGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTcc                                      <  1:2783220/70‑1 (MQ=255)
                           ccaGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGAc                                   <  1:588386/71‑1 (MQ=255)
                               cgcCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTgg                                                   <  1:559640/51‑1 (MQ=255)
                               cgcCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATgg                               <  1:2864776/71‑1 (MQ=255)
                                     ggTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCttt                          <  1:2750445/70‑1 (MQ=255)
                                     ggTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCtt                           <  1:3038601/69‑1 (MQ=255)
                                       tAGACATCGAAGATATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGAc                                   <  1:2583708/59‑1 (MQ=255)
                                        aGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGATTTGACg                      >  1:1145339/1‑71 (MQ=255)
                                         gACATCGAAGAGATCACCGGTAACACCGTTCGCGAAATTGGCTa                                                <  1:2314056/44‑1 (MQ=255)
                                         gACATCGAAGAGATCACCCGTAACACCGTTCGCGaaa                                                       >  1:985366/1‑37 (MQ=255)
                                         gACATCGAAGAGATCACCCGTAACACCGTTCGCGaaa                                                       >  1:2377879/1‑37 (MQ=255)
                                         gACATCGAAGAGATCACCCGTAACACCGTTCGCGa                                                         <  1:1269237/35‑1 (MQ=255)
                                         gACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTa                                                <  1:2221237/44‑1 (MQ=255)
                                         gACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTa                                                <  1:3839778/44‑1 (MQ=255)
                                           cATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTa                   <  1:1254209/71‑1 (MQ=255)
                                           cATCGAAGAGATCACCCGTAACAACGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTa                   <  1:529701/71‑1 (MQ=255)
                                             tCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTAtgtg                                            <  1:2212075/44‑1 (MQ=255)
                                             tCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTAtgtg                                            <  1:3283588/44‑1 (MQ=255)
                                              cGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATgg                               >  1:2235251/1‑56 (MQ=255)
                                               gAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACt                <  1:286033/70‑1 (MQ=255)
                                               gAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACt                <  1:2959297/70‑1 (MQ=255)
                                               gAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTc               <  1:1184937/71‑1 (MQ=255)
                                               gAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTc               >  1:2494529/1‑71 (MQ=255)
                                                aaGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTcc              <  1:3628799/71‑1 (MQ=255)
                                                aaGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTcc              <  1:3150228/71‑1 (MQ=255)
                                                 agagATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGCCGCTAACTCCt             >  1:2111560/1‑71 (MQ=255)
                                                 agagATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCt             >  1:3101173/1‑71 (MQ=255)
                                                          ccGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCttt                          <  1:3252773/49‑1 (MQ=255)
                                                          ccGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    >  1:2708140/1‑71 (MQ=255)
                                                          ccGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    >  1:1809825/1‑71 (MQ=255)
                                                           cGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACt                <  1:1318124/58‑1 (MQ=255)
                                                           cGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTCt   >  1:2963983/1‑71 (MQ=255)
                                                            ttAATACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    <  1:2923890/68‑1 (MQ=255)
                                                            gTAATACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    <  1:810474/69‑1 (MQ=255)
                                                            gTAATACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    <  1:3615424/69‑1 (MQ=255)
                                                            gTAATACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    <  1:422832/69‑1 (MQ=255)
                                                            gTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTCt   >  1:1675122/1‑70 (MQ=255)
                                                            gTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTCTg  <  1:624306/71‑1 (MQ=255)
                                                             tAATACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTc    <  1:2024643/68‑1 (MQ=255)
                                                                |                                                                  
TGGTTTTAGTTGGCGGCGAAATCACCACCAGCGCCTGGGTAGACATCGAAGAGATCACCCGTAACACCGTTCGCGAAATTGGCTATGTGCATTCCGACATGGGCTTTGACGCTAACTCCTGTGCGGTTCTG  >  W3110S.gb/3085510‑3085640

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: