Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 218,182 G→A 19.3% A57T (GCG→ACG)  yafF → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported218,1820GA19.3% 193.9 / 37.4 93A57T (GCG→ACG) yafFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (46/29);  new base A (10/8);  total (56/37)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACA  >  Exported/218113‑218250
                                                                     |                                                                    
ttCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATTTGGCGTCAGTCCTTGc                                                                     <  1:1664669/71‑1 (MQ=255)
 tCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                                                    >  1:5050505/1‑71 (MQ=255)
  cAAGGCATGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                                                   >  1:2166309/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCgg                                                                   >  1:2815696/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCgg                                                                   >  1:4845463/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCgg                                                                   >  1:792673/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCg                                                                    <  1:1015791/70‑1 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:1346359/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:3109827/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:1789837/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:2598466/1‑71 (MQ=255)
  cAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:380270/1‑71 (MQ=255)
   aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCgg                                                                   <  1:5237182/70‑1 (MQ=255)
   aaGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgg                                                                   >  1:2833266/1‑70 (MQ=255)
     ggCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGa                                                                <  1:4380476/71‑1 (MQ=255)
      gCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAg                                                               >  1:548341/1‑71 (MQ=255)
        aGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACg                                                                    >  1:1944452/1‑64 (MQ=255)
         gggTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCgg                                                            >  1:3965276/1‑71 (MQ=255)
         gggTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTAGCGTCAGTCCTTGCGGGGAGCgg                                                            >  1:3351741/1‑71 (MQ=255)
           gTTAAGACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCgggg                                                                 <  1:5216980/64‑1 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCggg                                                                  <  1:710175/62‑1 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGc                                                          >  1:3009295/1‑70 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                         >  1:5076748/1‑71 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                         >  1:3987943/1‑71 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCt                                                         >  1:1181451/1‑71 (MQ=255)
            ttAAGACGTAAGATGCGAAAAGAAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCggg                                                                  <  1:3304238/62‑1 (MQ=255)
              aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCtt                                                        <  1:2727779/70‑1 (MQ=255)
              aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                                      >  1:852468/1‑56 (MQ=255)
              aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTg                                                                      >  1:1125942/1‑56 (MQ=255)
              aaGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                                                       >  1:899542/1‑71 (MQ=255)
                gACGTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                                     <  1:1388404/71‑1 (MQ=255)
                gACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                                     >  1:3774724/1‑71 (MQ=255)
                 aCGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTGGCGGGGAGCGGGCTTTCgt                                                    <  1:3962684/71‑1 (MQ=255)
                 aCGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgt                                                    <  1:5643780/71‑1 (MQ=255)
                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGc                                                                     >  1:745403/1‑53 (MQ=255)
                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                                   <  1:2202383/71‑1 (MQ=255)
                  cGTAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                                   <  1:1786494/71‑1 (MQ=255)
                   gTAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaa                                                  >  1:4585127/1‑70 (MQ=255)
                    tAAGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCg                                                                    <  1:664308/52‑1 (MQ=255)
                    tAAGATGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaa                                                  <  1:4927959/70‑2 (MQ=255)
                    tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtag                                                  >  1:1978372/1‑70 (MQ=255)
                    tAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACgggg                                                                 <  1:3460643/55‑1 (MQ=255)
                     aaGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                      >  1:2716047/1‑65 (MQ=255)
                     aaGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                      >  1:4112011/1‑65 (MQ=255)
                     aaGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGAGGGGAGCGGGCTTTc                                                      >  1:4938711/1‑65 (MQ=255)
                     aaGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgta                                                   >  1:3841719/1‑68 (MQ=255)
                      aGATGCGAAAAGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCggg                                                                  >  1:3489218/1‑52 (MQ=255)
                       gATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTc                                                      >  1:2241696/1‑63 (MQ=255)
                        aTGCGAAAAGCAGCCATGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCg                                                     >  1:5043737/1‑63 (MQ=255)
                        aTGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACg                                                                    <  1:312583/48‑1 (MQ=255)
                        aTGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACg                                                                    <  1:2024620/48‑1 (MQ=255)
                        aTGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACg                                                                    <  1:2830334/48‑1 (MQ=255)
                            gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtagta                                                >  1:654207/1‑64 (MQ=255)
                            gAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATAcccc                                         <  1:22903/71‑1 (MQ=255)
                             aaaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:4200298/1‑44 (MQ=255)
                             aaaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:2010563/1‑44 (MQ=255)
                             aaaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:3709748/1‑44 (MQ=255)
                             aaaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:3215063/1‑44 (MQ=255)
                             aaaaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:3008033/1‑44 (MQ=255)
                               aaGCAGCCATGGATAGAAACTATCTCGCGTCAGTCCTTGCggg                                                                  <  1:4944522/43‑1 (MQ=255)
                               aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCgg                                                                   >  1:640366/1‑42 (MQ=255)
                               aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCgg                                                            >  1:2270936/1‑49 (MQ=255)
                               aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCg                                                             >  1:1994317/1‑48 (MQ=255)
                               aaGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCtt                                                        <  1:418260/53‑1 (MQ=255)
                                 gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTGCg                                                                    <  1:4232697/39‑1 (MQ=255)
                                 gcagcCATGGACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTc                                                      >  1:454470/1‑53 (MQ=255)
                                     ccATGGACAGAAACTATCTGGCGTCAGTCCTTGCgg                                                                   <  1:1481279/36‑1 (MQ=255)
                                        tGGATAGAAACTATCTCGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaa                                                  <  1:1025200/50‑2 (MQ=255)
                                        tGGACAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                                                       <  1:5545904/45‑1 (MQ=255)
                                        tGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGt                             <  1:32265/71‑1 (MQ=255)
                                          gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCg                                                             >  1:5167082/1‑37 (MQ=255)
                                          gACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGa                                       <  1:269346/59‑1 (MQ=255)
                                          gACAGAAACTACCTGGCGTCAGTCCTTACGGGGAg                                                               >  1:5442247/1‑35 (MQ=255)
                                          gACAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaat                                                 >  1:4749352/1‑47 (MQ=39)
                                           aCAGAAACTATCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCgtaa                                                  >  1:588114/1‑46 (MQ=255)
                                           aCAGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTCgtaat                                                 <  1:1560499/48‑3 (MQ=39)
                                            taGAAACTAGCTGACGTAAGTCCTTGCGGGGAGCAGGGTTTCgta                                                   <  1:4048999/44‑1 (MQ=255)
                                            cAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCttt                                                       >  1:4446192/1‑41 (MQ=255)
                                            cAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGAct                                     >  1:1316173/1‑59 (MQ=255)
                                             aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATTCCCCGCCTCTcccc                               >  1:2861394/1‑64 (MQ=255)
                                             aGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTa                        <  1:1925957/71‑1 (MQ=255)
                                             aGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTc                                                      >  1:5542994/1‑41 (MQ=255)
                                             aGAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCTTTc                                                      >  1:4544381/1‑41 (MQ=255)
                                              gAAACTACCTGGCGTCAGTCCTTACGGGGAGCGGGCtt                                                        <  1:5284256/38‑1 (MQ=255)
                                                   tACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCgtaaa                                                 <  1:4986723/40‑3 (MQ=255)
                                                   tACCTGACGTCAGTCCTTGCGGGGAGCAGGCTTTCgtaaa                                                 <  1:1559698/40‑3 (MQ=255)
                                                    aCCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTaa                                               <  1:5487540/41‑1 (MQ=255)
                                                        ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGCCTCTcccc                               >  1:3545847/1‑53 (MQ=255)
                                                        ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc                            >  1:1927930/1‑56 (MQ=255)
                                                        ggCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACACCGa                                       >  1:4164319/1‑45 (MQ=255)
                                                            tCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAacac                   <  1:2487951/61‑1 (MQ=255)
                                                               gTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGAc                                      >  1:2934784/1‑39 (MQ=255)
                                                                  cTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCgg                             >  1:3844474/1‑44 (MQ=255)
                                                                  cTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTc                            >  1:3270575/1‑46 (MQ=255)
                                                                    tGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTcacc     <  1:3697795/67‑1 (MQ=255)
                                                                     gCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACcaca  <  1:3447389/69‑1 (MQ=255)
                                                                     gCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACcaca  <  1:4680160/69‑1 (MQ=255)
                                                                     |                                                                    
TTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCAGCCATGGACAGAAACTACCTGGCGTCAGTCCTTGCGGGGAGCGGGCTTTCGTAGTAATACCCCGACTCTCCCCGTCCTTAAACACAACCCCCACTCACCACA  >  Exported/218113‑218250

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: