Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 217,903 A→G 8.5% S2S (AGT→AGC yafU ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported217,9030AG 8.5% 206.6 / 8.3 70S2S (AGT→AGCyafUpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (19/45);  new base G (3/3);  total (22/48)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTG  >  Exported/217836‑217973
                                                                   |                                                                      
cGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTc                                                                     <  1:4502323/71‑1 (MQ=255)
 ggCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCa                                                                    <  1:1938761/71‑1 (MQ=255)
    tataACTGCTTTGGCCACATCCATCGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTc                                                                     <  1:3817260/67‑1 (MQ=255)
      taACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCataat                                                                <  1:1634210/70‑1 (MQ=255)
      taACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCataat                                                                <  1:2817154/70‑1 (MQ=255)
            cTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTc                                                          <  1:2241739/70‑1 (MQ=255)
               tGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTa                                                        <  1:873863/69‑1 (MQ=255)
               tGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTa                                                        <  1:1989103/69‑1 (MQ=255)
                       tCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCa                                                                    <  1:2148059/49‑1 (MQ=255)
                       tCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTctgct                                              <  1:3210212/71‑1 (MQ=255)
                       tCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTctgc                                               <  1:3420465/70‑1 (MQ=255)
                         cATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCata                                                                  >  1:2801384/1‑49 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCataat                                                                <  1:3508506/49‑1 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCataat                                                                <  1:2838617/49‑1 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCataat                                                                <  1:2145853/49‑1 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTc                                                  >  1:1795562/1‑63 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTc                                                  >  1:4109501/1‑63 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTc                                                  >  1:4269031/1‑63 (MQ=255)
                           ttGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTc                                                  >  1:1022936/1‑63 (MQ=255)
                            tGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTc                                                  >  1:954280/1‑62 (MQ=255)
                            tGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGAtt                                          >  1:2119080/1‑70 (MQ=255)
                             gAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCAtaataa                                                              >  1:3630099/1‑49 (MQ=255)
                             gAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCAtaataa                                                              >  1:3187885/1‑49 (MQ=255)
                             gAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGAttt                                         >  1:790356/1‑70 (MQ=255)
                             gAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGAttt                                         >  1:4722093/1‑70 (MQ=255)
                              aaaaaTCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTaa                                       <  1:3725578/71‑1 (MQ=255)
                              aaaaaTCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTa                                        <  1:2556906/70‑1 (MQ=255)
                                  aTCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCg                                    >  1:4734730/1‑70 (MQ=255)
                                   tCGCCAAGAAAATTAACCAGATCTCGCTCTGAGCTCATAATAAATTCTa                                                        >  1:4246712/1‑49 (MQ=255)
                                   tCGCCAAGAAAATTAACCAGATCTCGCTCTGAGCTCATAATAAATTCTa                                                        >  1:1716179/1‑49 (MQ=255)
                                   tCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCt                                                         <  1:4163278/48‑1 (MQ=255)
                                   tCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCa                                  <  1:3892677/71‑1 (MQ=255)
                                   tCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCa                                  <  1:3409565/71‑1 (MQ=255)
                                     gCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCag                                 <  1:2375737/70‑1 (MQ=255)
                                     gCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCag                                 <  1:324417/70‑1 (MQ=255)
                                      ccAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:178353/71‑1 (MQ=255)
                                      ccAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:2429857/71‑1 (MQ=255)
                                        aaGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                             <  1:2768644/71‑1 (MQ=255)
                                        aaGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGaa                             <  1:298933/71‑1 (MQ=255)
                                           aaaaTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATg                                                      >  1:4071263/1‑43 (MQ=255)
                                              aTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGc                       <  1:3928627/71‑1 (MQ=255)
                                              aTTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGc                       <  1:3891925/71‑1 (MQ=255)
                                                  aCCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCcac                    <  1:3976307/70‑1 (MQ=255)
                                                  aCCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCcac                    <  1:869333/70‑1 (MQ=255)
                                                  aCCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCcac                    <  1:3022215/70‑1 (MQ=255)
                                                      gATCTCGCTCTGAGCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa               >  1:4423411/1‑71 (MQ=255)
                                                      gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCg                                    <  1:1650124/50‑1 (MQ=255)
                                                      gATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa               <  1:689928/71‑1 (MQ=255)
                                                       aTCTCGCTCTGAGCTCATAATAAATTCTATGGTTctgct                                              <  1:4746406/39‑1 (MQ=255)
                                                       aTCTCGCTCTGAGCTCATAATAAATTCTATGGTTctgct                                              <  1:588247/39‑1 (MQ=255)
                                                       aTCTCGCTCTGAGCTCATAATAAATTCTATGGTTctgct                                              <  1:4100523/39‑1 (MQ=255)
                                                       aTCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:4076849/54‑1 (MQ=255)
                                                       aTCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:2206816/54‑1 (MQ=255)
                                                       aTCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCag                                 <  1:482969/52‑1 (MQ=255)
                                                       aTCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCaa               <  1:2362654/70‑1 (MQ=255)
                                                         ctcGCTCTGAACTCATAATAAATTCTATGGTTctg                                                <  1:1922497/35‑1 (MQ=255)
                                                         ctcGCTCTGAACTCATAATAAATTCTATGGTTctg                                                <  1:555218/35‑1 (MQ=255)
                                                            gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAAcc                                     >  1:1468779/1‑43 (MQ=255)
                                                            gCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAAc                                      >  1:385278/1‑42 (MQ=255)
                                                               ctGAACTCATAATAAATTCTATGGTTCTGCTGATTTaa                                       >  1:393140/1‑38 (MQ=255)
                                                               ctGAACTCATAATAAATTCTATGGTTCTGCTGATTTaa                                       >  1:656978/1‑38 (MQ=255)
                                                               ctGAACTCATAATAAATTCTATGGTTCTGCTGATTTaa                                       >  1:3727446/1‑38 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:395733/45‑1 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:1805753/45‑1 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:2422136/45‑1 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCagag                               <  1:2641005/45‑1 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc       <  1:3571912/69‑1 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc       <  1:2346513/69‑1 (MQ=255)
                                                                tGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAc       <  1:1200839/69‑1 (MQ=255)
                                                                 gAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCAc    >  1:1951054/1‑71 (MQ=255)
                                                                  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGaaa        >  1:2584510/1‑66 (MQ=255)
                                                                  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt   >  1:3851588/1‑71 (MQ=255)
                                                                  aaCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACt   >  1:401005/1‑71 (MQ=255)
                                                                   aCTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTg  <  1:1679835/71‑1 (MQ=255)
                                                                   |                                                                      
CGGCTATAACTGCTTTGGCCACATCCATTGAAAAATCGCCAAGAAAATTAACCAGATCTCGCTCTGAACTCATAATAAATTCTATGGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTG  >  Exported/217836‑217973

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: