Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 209,541 T→C 5.2% R55R (CGA→CGG yafC ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported209,5410TC 5.2% 266.8 / 3.8 96R55R (CGA→CGGyafCpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (43/48);  new base C (2/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CATCTCCTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCTTACCGCTG  >  Exported/209475‑209611
                                                                  |                                                                      
cATCTCCTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTcgcg                                                                    <  1:1088095/71‑1 (MQ=255)
   ctcGTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTgg                                                                 >  1:2574399/1‑71 (MQ=255)
      cTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTcc                                                              <  1:759366/71‑1 (MQ=255)
      cTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTcc                                                              <  1:5459892/71‑1 (MQ=255)
      cTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTcc                                                              <  1:173158/71‑1 (MQ=255)
      cTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTc                                                               >  1:662960/1‑70 (MQ=255)
      cTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTc                                                               >  1:3694864/1‑70 (MQ=255)
       tGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTg                                                                  <  1:687314/66‑1 (MQ=255)
         cAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTgg                                                                 >  1:3327271/1‑65 (MQ=255)
          aaaaTTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCGGTCAGGCTAAGTTGTCGCGTGGTCCGAt                                                           >  1:513071/1‑70 (MQ=255)
          aaaaTTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGAt                                                           <  1:2023716/70‑1 (MQ=255)
             aTTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAg                                                       <  1:4588148/71‑1 (MQ=255)
             aTTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAg                                                       <  1:367208/71‑1 (MQ=255)
             aTTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAg                                                       <  1:4654808/71‑1 (MQ=255)
                gACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTa                                                         <  1:1570954/66‑1 (MQ=255)
                gACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTa                                                         <  1:4750227/66‑1 (MQ=255)
                gACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCa                                                     <  1:2988541/70‑1 (MQ=255)
                gACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCa                                                     <  1:4954075/70‑1 (MQ=255)
                   tGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGc                                                  >  1:2465467/1‑70 (MQ=255)
                     tACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTTGTCCGATTAAGCAGGCt                                                 <  1:4734939/69‑1 (MQ=255)
                     tACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCt                                                 <  1:3454290/69‑1 (MQ=255)
                      aCGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAAc                                              >  1:2312154/1‑71 (MQ=255)
                       cgcgACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACg                                             >  1:5042373/1‑71 (MQ=255)
                          gacgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTaa                                                        <  1:2340004/57‑1 (MQ=255)
                          gacgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTaa                                                        <  1:5100870/57‑1 (MQ=255)
                          gacgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGcc                                           <  1:5157920/70‑1 (MQ=255)
                           acgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTcgc                                                                     <  1:3747961/43‑1 (MQ=255)
                           acgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTcgc                                                                     <  1:4738373/43‑1 (MQ=255)
                           acgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTcgc                                                                     <  1:1508570/43‑1 (MQ=255)
                            cgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTaa                                                        <  1:2130333/55‑1 (MQ=255)
                            cgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGc                                                      >  1:2031126/1‑57 (MQ=255)
                            cgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGc                                                      >  1:2841373/1‑57 (MQ=255)
                            cgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGc                                                      >  1:5007591/1‑57 (MQ=255)
                            cgaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGc                                                      >  1:5106460/1‑57 (MQ=255)
                             gaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTcgc                                                                     <  1:420302/41‑1 (MQ=255)
                             gaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGt                                                                   >  1:5116137/1‑43 (MQ=255)
                             gaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGt                                                                   >  1:763167/1‑43 (MQ=255)
                             gaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGt                                       >  1:1047471/1‑71 (MQ=255)
                             gaAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGt                                       >  1:896993/1‑71 (MQ=255)
                              aaaaTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGtt                                      <  1:5018679/71‑1 (MQ=255)
                               aaaTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGttt                                     >  1:1667703/1‑71 (MQ=255)
                               aaaTAACGCTCACCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGttt                                     >  1:5406325/1‑71 (MQ=255)
                                aaTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGttt                                     >  1:2135005/1‑70 (MQ=255)
                                aaTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGttt                                     >  1:4433701/1‑70 (MQ=255)
                                aaTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    <  1:1973749/71‑1 (MQ=255)
                                 aTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAg                                                    >  1:1252558/1‑54 (MQ=255)
                                 aTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAg                                                    >  1:3483075/1‑54 (MQ=255)
                                 aTAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAg                                                    >  1:4431107/1‑54 (MQ=255)
                                  tAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTg                                                                  >  1:4988832/1‑39 (MQ=255)
                                  tAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTg                                                                  >  1:2629077/1‑39 (MQ=255)
                                  tAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    >  1:103536/1‑69 (MQ=255)
                                  tAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCa                                   >  1:670310/1‑70 (MQ=255)
                                  tAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCAt                                  >  1:371287/1‑71 (MQ=255)
                                   gaCGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTa                                                         >  1:5220229/2‑47 (MQ=255)
                                   aaCGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCAtc                                 >  1:3691314/1‑71 (MQ=255)
                                    aCGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAg                                                       >  1:3382130/1‑48 (MQ=255)
                                    aCGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCAt                                  >  1:1565431/1‑69 (MQ=255)
                                     cGCTCGCCTTCTTCCGTCAGGCTAAGTTGCCGCGTGGTTCGATTAAGCAGGCTAACGCCAAGt                                       >  1:3700032/1‑63 (MQ=255)
                                     cGCTCGCCTTCTTCCGTCAGGCTAAGTTGCCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGt                                       >  1:3722674/1‑63 (MQ=255)
                                        tCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTaa                                                        <  1:5658430/43‑1 (MQ=255)
                                         cGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCt                                                 <  1:1718608/49‑1 (MQ=255)
                                         cGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCaa                                         >  1:2513183/1‑57 (MQ=255)
                                          gCCTTCTTCCGTCAGGCTAAGTTGTCGCGTTGTCCGATTAAGCAGGCTAACg                                             <  1:3596287/52‑1 (MQ=255)
                                          gCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGc                                                      >  1:3778821/1‑43 (MQ=255)
                                          gCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACg                                             <  1:2833671/52‑1 (MQ=255)
                                           ccTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTaa                                               >  1:1570765/1‑49 (MQ=255)
                                           ccTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTaa                                               >  1:303585/1‑49 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCTCGTGGTCCGATTAAGc                                                      >  1:1404287/1‑41 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGcc                                           >  1:5663420/1‑52 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGcc                                           >  1:20408/1‑52 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttt                        <  1:3896621/71‑1 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttt                        <  1:2983641/71‑1 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttt                        <  1:1633722/71‑1 (MQ=255)
                                            cttcttCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttt                        <  1:5683309/71‑1 (MQ=255)
                                             ttcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCt                                                 <  1:5633377/45‑1 (MQ=255)
                                             ttcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCt                                                 <  1:4876923/45‑1 (MQ=255)
                                             ttcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCt                                                 <  1:4497067/45‑1 (MQ=255)
                                             ttcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAg                                        <  1:3961342/54‑1 (MQ=255)
                                             ttcttcCGTCAGGCTAAGTTGCCGCGTGGTCCGATTAAGCAGGCTAACGCCAAg                                        <  1:2928746/54‑1 (MQ=255)
                                             ttcttcCGTCAGGCTAAGTTGCCGCGTGGTCCGATTAAGCAGGCTAACGCCAAg                                        <  1:3578392/54‑1 (MQ=255)
                                              tcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAAc                                              >  1:1377200/1‑47 (MQ=255)
                                              tcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAAc                                              >  1:3365149/1‑47 (MQ=255)
                                              tcttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttt                        <  1:3832184/69‑1 (MQ=255)
                                               cttcCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttttt                      <  1:3684117/70‑1 (MQ=255)
                                                     tCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    <  1:402011/50‑1 (MQ=255)
                                                     tCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    <  1:3858028/50‑1 (MQ=255)
                                                       aGGCTAAGTTTCCGCGTGGTCCGATTAAGCAGGCTAACGCCAAg                                        <  1:4909142/44‑1 (MQ=38)
                                                          cTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCCCCGGCtt         <  1:4486432/71‑1 (MQ=255)
                                                             aGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCtt         <  1:2082558/69‑1 (MQ=255)
                                                             aGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCtt         <  1:2663649/69‑1 (MQ=255)
                                                             aGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCTTAc       <  1:3341270/71‑1 (MQ=255)
                                                              gTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCttttt                      <  1:365602/55‑1 (MQ=255)
                                                                tGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCa                    >  1:1185671/1‑55 (MQ=255)
                                                                tGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCTTACCGc    >  1:1705760/1‑71 (MQ=255)
                                                                 gTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    <  1:3331801/38‑1 (MQ=255)
                                                                 gTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    <  1:1380838/38‑1 (MQ=255)
                                                                  tCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTc                                    <  1:3882937/37‑1 (MQ=255)
                                                                  tCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCTTACCGCTg  <  1:5452112/71‑1 (MQ=255)
                                                                  |                                                                      
CATCTCCTGCAAAATTGACTGTACGCGACGAAAATAACGCTCGCCTTCTTCCGTCAGGCTAAGTTGTCGCGTGGTCCGATTAAGCAGGCTAACGCCAAGTTTCATCTCCAGCTTTTTCACCGCCCGGCTTACCGCTG  >  Exported/209475‑209611

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: