Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,354,061 C→T 24.6% G218G (GGG→GGA dcuR ← DNA‑binding response regulator in two‑component regulatory system with DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,354,0610CT24.6% 114.0 / 35.2 61G218G (GGG→GGAdcuRDNA‑binding response regulator in two‑component regulatory system with DcuS
Reads supporting (aligned to +/- strand):  ref base C (20/26);  new base T (7/8);  total (27/34)
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

ATTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTA  >  W3110S.gb/4353997‑4354127
                                                                |                                                                  
aTTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                               <  1:2402788/70‑1 (MQ=255)
aTTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                               <  1:2804789/70‑1 (MQ=255)
  tGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGTCCCGTGa                                                             <  1:663008/70‑1 (MQ=255)
  tGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGTCCCGTGAc                                                            >  1:750925/1‑71 (MQ=255)
        tatTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                               <  1:4012506/62‑1 (MQ=255)
         atTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATaa                                                     >  1:1082426/1‑71 (MQ=255)
         atTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATaa                                                     >  1:1012847/1‑71 (MQ=255)
         atTGTTACAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAtga                                                     >  1:524997/1‑69 (MQ=255)
          ttGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAAt                                                    >  1:1587165/1‑71 (MQ=255)
            gTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATg                                                   >  1:3745922/1‑70 (MQ=255)
             tttCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGa                                                 <  1:3971928/71‑1 (MQ=255)
              ttCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGa                                                 <  1:3970384/70‑1 (MQ=255)
               tCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGAt                                                <  1:4052377/70‑1 (MQ=255)
                 agtagtgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACt                                             <  1:740287/71‑1 (MQ=255)
                     gtgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGt                                          >  1:2544896/1‑70 (MQ=255)
                     gtgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGt                                          >  1:2253625/1‑70 (MQ=255)
                      tgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                       <  1:2974124/56‑1 (MQ=255)
                      tgagtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                       <  1:1862307/56‑1 (MQ=255)
                        agtagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGaa                                       >  1:2072656/1‑70 (MQ=255)
                          tagtgTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACaa                                    <  1:2583058/71‑1 (MQ=255)
                            gtgtTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAg                                   <  1:66804/70‑1 (MQ=255)
                             tgtTCTGCCTGAATACGGTAAGGGTAAACTGGACGCCCCGTGACGCCAt                                                       <  1:2111961/49‑1 (MQ=255)
                               ttCTGCCTGAATACGGTAACTGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGAtgt                               >  1:1587494/1‑71 (MQ=255)
                               ttCTGCCTGAATACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGa                                  >  1:425130/1‑68 (MQ=255)
                               ttCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACaa                                    <  1:1304211/66‑1 (MQ=255)
                               ttCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGAtgt                               >  1:3030864/1‑71 (MQ=255)
                               atCTGCCTGAATACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGa                                  >  1:1119529/2‑68 (MQ=255)
                                 cTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATa                                               >  1:1773974/1‑53 (MQ=255)
                                 cTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATa                                               >  1:2823538/1‑53 (MQ=255)
                                 cTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATa                                               >  1:443453/1‑53 (MQ=255)
                                  tGCCTGAATACGGTAACGGTAAACTGGACGCCCCGt                                                               <  1:1870927/36‑1 (MQ=255)
                                  tGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                       <  1:599014/44‑1 (MQ=255)
                                   gCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGAAGCCATAATGGATACTGGTGAACAAGATGTgg                             >  1:3434368/1‑69 (MQ=255)
                                     cTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATg                                                   <  1:426098/45‑1 (MQ=255)
                                       taaTACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                       <  1:2271770/38‑1 (MQ=255)
                                       gAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGa                                                 >  1:543522/1‑45 (MQ=255)
                                        aaTACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGtt                        <  1:619085/69‑1 (MQ=255)
                                        aaTACGGTAACGGTAAACTGGACGCCCCGTGACGCCAt                                                       <  1:1179409/38‑1 (MQ=255)
                                         aTACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGa                      <  1:3737938/70‑1 (MQ=255)
                                         aTACGGTAACGGTAAACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAc                     <  1:1114119/71‑1 (MQ=255)
                                          tACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAg                                   <  1:618192/56‑1 (MQ=255)
                                          tACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAg                                   <  1:581810/56‑1 (MQ=255)
                                            cggtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccag                  >  1:385942/1‑71 (MQ=255)
                                             ggtaacggtaaACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc                 >  1:3174016/1‑71 (MQ=255)
                                             ggtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc                 <  1:3433783/71‑1 (MQ=255)
                                              gtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc                 <  1:1078931/70‑1 (MQ=255)
                                              gtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc                 >  1:1358491/1‑70 (MQ=255)
                                              gtaacggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagc                 >  1:1363297/1‑70 (MQ=255)
                                               taacggtaaACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAccagcca               <  1:3997180/71‑1 (MQ=255)
                                                aacggtaaACTGGACGCCCCGTGACGCCATAATGGa                                                 >  1:4077479/1‑36 (MQ=255)
                                                 acggtaaACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAc                     >  1:2536742/1‑63 (MQ=255)
                                                   ggtaaACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGAt            >  1:1514323/1‑70 (MQ=255)
                                                    gtaaACTGGACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGAc                     <  1:2570200/60‑1 (MQ=255)
                                                      aaaCTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGAtg                                <  1:1364988/47‑1 (MQ=255)
                                                        aCTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTa      <  1:2431071/71‑1 (MQ=255)
                                                            gACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTAtt    <  1:1178893/69‑1 (MQ=255)
                                                            gACGTCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTa  >  1:3163715/1‑71 (MQ=255)
                                                            gACGTCACGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTa  >  1:3417910/1‑71 (MQ=255)
                                                            gACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTAtt    <  1:414579/69‑1 (MQ=255)
                                                            gACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTAtt    <  1:4008872/69‑1 (MQ=255)
                                                            gACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTAtt    <  1:2102045/69‑1 (MQ=255)
                                                                |                                                                  
ATTGGCAATATTGTTTCAGTAGTGAGTAGTGTTCTGCCTGAATACGGTAACGGTAAACTGGACGCCCCGTGACGCCATAATGGATACTGGTGAACAAGATGTGGCAGTTGACCAGCCAGATGAGGTATTTA  >  W3110S.gb/4353997‑4354127

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: