Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 75,753 A→C 5.6% L516W (TTG→TGG)  sgrR ← DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP00904875,7530AC 5.6% 230.0 / 6.1 87L516W (TTG→TGG) sgrRDNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/52);  new base C (3/2);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGC  >  AP009048/75686‑75822
                                                                   |                                                                     
ccAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaa                                                                      >  1:6800218/1‑69 (MQ=255)
 cAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataa                                                                   <  1:6878711/71‑1 (MQ=255)
   gCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATg                                                                 >  1:7548464/1‑71 (MQ=255)
    ccGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATgg                                                                <  1:6732549/71‑1 (MQ=255)
        gggTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCAc                                                             >  1:7073934/1‑70 (MQ=255)
         ggTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCCATTATTGCACCa                                                           <  1:2615056/71‑1 (MQ=255)
         ggTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCCATTATTGCACCa                                                           <  1:1609628/71‑1 (MQ=255)
          gTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCAt                                                          <  1:5748295/71‑1 (MQ=255)
            aTTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCg                                                        >  1:4216873/1‑71 (MQ=255)
             ttCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataat                                                                  >  1:742786/1‑60 (MQ=255)
             ttCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCg                                                        <  1:1080724/70‑1 (MQ=255)
              tCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaata                                                                    <  1:3459127/57‑1 (MQ=255)
               cATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                     >  1:809611/1‑71 (MQ=255)
               cATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                     >  1:5578801/1‑71 (MQ=255)
               cATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                     >  1:4388715/1‑71 (MQ=255)
               cATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCt                                                      <  1:4460641/70‑1 (MQ=255)
               cATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCt                                                      <  1:2221445/70‑1 (MQ=255)
                aTGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                    >  1:3638427/1‑71 (MQ=255)
                aTGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                    >  1:2460333/1‑71 (MQ=255)
                aTGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                     <  1:6830172/70‑1 (MQ=255)
                aTGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCtt                                                     <  1:1666266/70‑1 (MQ=255)
                 tGCGCGGCCTGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATgg                                                                >  1:1077611/1‑58 (MQ=255)
                 tGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTg                                                   <  1:3875815/71‑1 (MQ=255)
                  gcgcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGc                                                               <  1:2491260/58‑1 (MQ=255)
                   cgcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaat                                                                     >  1:7709317/1‑51 (MQ=255)
                   cgcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                    <  1:7496157/68‑1 (MQ=255)
                   cgcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                    <  1:2760387/68‑1 (MQ=255)
                   cgcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCttt                                                    <  1:7486703/68‑1 (MQ=255)
                    gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                <  1:7251755/71‑1 (MQ=255)
                    gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                <  1:3870018/71‑1 (MQ=255)
                    gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                <  1:6202147/71‑1 (MQ=255)
                    gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                <  1:2187251/71‑1 (MQ=255)
                    gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                >  1:2762508/1‑71 (MQ=255)
                    gcAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                <  1:2630400/71‑1 (MQ=255)
                       ggCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCa                                                                       <  1:2942866/45‑1 (MQ=255)
                       ggCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCa                                                                       <  1:3050710/45‑1 (MQ=255)
                       ggCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCa                                                                       <  1:4170588/45‑1 (MQ=255)
                       ggCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCg                                             <  1:4500778/71‑1 (MQ=255)
                        gCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataat                                                                  <  1:4008370/49‑1 (MQ=255)
                        gCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaat                                                                     <  1:1824532/46‑1 (MQ=255)
                        gCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCAcc                                                            <  1:2779730/55‑1 (MQ=255)
                        gCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCAcc                                                            <  1:246369/55‑1 (MQ=255)
                          cgcgCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGAcc                                          <  1:6359796/71‑1 (MQ=255)
                           gcgcATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGAc                                           >  1:4856983/1‑69 (MQ=255)
                             gcATACTGCGTTGCCCCTGAATGCTCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGt                                       <  1:7203392/71‑1 (MQ=255)
                              cATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCt                                                 >  1:2731250/1‑60 (MQ=255)
                              cATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGAc                                           >  1:3176171/1‑66 (MQ=255)
                              cATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGAc                                           >  1:6628406/1‑66 (MQ=255)
                              cATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGtt                                      <  1:5069654/71‑1 (MQ=255)
                              cATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGtt                                      >  1:4511371/1‑71 (MQ=255)
                                  cTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCt                                   <  1:6697328/70‑1 (MQ=255)
                                  cTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTg                                  <  1:7547189/71‑1 (MQ=255)
                                  cTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTg                                  <  1:5487221/71‑1 (MQ=255)
                                   tGCGTTGCCCCTGAATGATCAGCCAGTGGTGCaataa                                                                   >  1:5310431/1‑37 (MQ=255)
                                   tGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAg                                        >  1:3104677/1‑64 (MQ=255)
                                       ttGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGAcc                                          <  1:5145577/58‑1 (MQ=255)
                                       ttGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCAcc                             <  1:785689/71‑1 (MQ=255)
                                        tGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGc                                                  <  1:5289485/49‑1 (MQ=255)
                                        tGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCAc                              >  1:1608540/1‑69 (MQ=255)
                                         gCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGc                                                  >  1:1905806/1‑48 (MQ=255)
                                          ccccTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGc                                                  <  1:1873544/47‑1 (MQ=255)
                                          ccccTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTg                                                >  1:7436635/1‑49 (MQ=255)
                                          ccccTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGt                          <  1:3189990/71‑1 (MQ=255)
                                          ccccTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGt                          <  1:4833361/71‑1 (MQ=255)
                                           cccTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGc                                <  1:7343092/64‑1 (MQ=255)
                                           cccTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGc                                <  1:6578629/64‑1 (MQ=255)
                                            acTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGc                                              >  1:658808/2‑49 (MQ=255)
                                              tGAATGATCAGCCAGTGGTGCAATAATGGGACCATCGCTTTGCTGGCGACCa                                         <  1:2925440/52‑1 (MQ=255)
                                              tGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGtt                                      <  1:2119949/55‑1 (MQ=255)
                                              tGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGtt                                      <  1:239497/55‑1 (MQ=255)
                                                aaTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGc                                    <  1:3804604/55‑1 (MQ=255)
                                                aaTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCa                               >  1:5694171/1‑60 (MQ=255)
                                                aaTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGc                      <  1:112404/69‑1 (MQ=255)
                                                aaTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACAAGTTGc                                    <  1:1017543/55‑1 (MQ=255)
                                                 aTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCATTTGCTGGCACCAGTTCGCCAg                   <  1:2598686/71‑1 (MQ=255)
                                                 aTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGc                                    <  1:3351573/54‑1 (MQ=255)
                                                 aTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGc                      >  1:286389/1‑68 (MQ=255)
                                                 aTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAg                   <  1:2838987/71‑1 (MQ=255)
                                                 aTGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAg                   <  1:5900631/71‑1 (MQ=255)
                                                  tGATCAGCCAGTGGTGCAATAATGGCACGATCGCTTTGCTGGCGACCAGtgg                                     >  1:6457474/1‑50 (MQ=39)
                                                  tGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTc                        >  1:1601805/1‑65 (MQ=255)
                                                  tGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGc                      >  1:885734/1‑67 (MQ=255)
                                                  tGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGc                      >  1:808503/1‑67 (MQ=255)
                                                   gATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCt                                   <  1:1823851/53‑1 (MQ=255)
                                                    aTCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGAt                 <  1:1133288/70‑1 (MQ=255)
                                                     tCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTc               <  1:4782687/71‑1 (MQ=255)
                                                     tCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTc               <  1:251237/71‑1 (MQ=255)
                                                       aGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCa              <  1:5730817/70‑1 (MQ=255)
                                                       aGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCAt             >  1:2980817/1‑71 (MQ=255)
                                                             tggtgCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGtt                         <  1:6363266/53‑1 (MQ=255)
                                                             tggtgCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGcc                     <  1:700390/57‑1 (MQ=255)
                                                              ggtgCCATTATTGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGcc       >  1:6797929/1‑70 (MQ=255)
                                                               gtgCCATTATTGCACCATCGCTTTGCTGGCGACCAGTTg                                     >  1:305866/1‑39 (MQ=25)
                                                               gtgCCATTATTGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCAt     >  1:2093876/1‑71 (MQ=255)
                                                               gtgCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTc                        >  1:6042433/1‑52 (MQ=255)
                                                                 gCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTg   >  1:6047578/1‑71 (MQ=255)
                                                                  cAATAATGGCACCATCGCTTTGCTGGCGACCAGtt                                      <  1:7636995/35‑1 (MQ=255)
                                                                  cAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTgc  >  1:511097/1‑71 (MQ=255)
                                                                   aataatGGCACCATCGCTTTGCTGGAGACCAGTTGCTGGCACCAg                           >  1:1583328/1‑45 (MQ=255)
                                                                   |                                                                     
CCAGCCGAGGGTATTCATGCGCAGGCCGCGCATACTGCGTTGCCCCTGAATGATCAGCCAGTGGTGCAATAATGGCACCATCGCTTTGCTGGCGACCAGTTGCTGGCACCAGTTCGCCAGATTCATCTCGCCATTGC  >  AP009048/75686‑75822

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: