Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,849,373 A→T 5.5% V166E (GTG→GAG)  ygcH ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,849,3730AT 5.5% 194.7 / 2.1 70V166E (GTG→GAG) ygcHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/20);  new base T (2/2);  total (48/22)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGTTCC  >  Exported/1849305‑1849441
                                                                    |                                                                    
tCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGcac                                                                     <  1:1086084/70‑1 (MQ=255)
 cGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCacac                                                                   >  1:1780171/1‑71 (MQ=255)
  gATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGcaca                                                                    >  1:3192397/1‑69 (MQ=255)
  gATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACAcc                                                                  >  1:3080663/1‑71 (MQ=255)
  gATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACAcc                                                                  >  1:3759596/1‑71 (MQ=255)
   aTTTAGCTGGCCAAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCt                                                                 <  1:1069886/71‑1 (MQ=255)
    tttAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                <  1:284162/71‑1 (MQ=255)
    tttAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                <  1:1518786/71‑1 (MQ=255)
    tttAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGAAGGTGTGCTCAcca                                                                 <  1:2642265/70‑2 (MQ=255)
    tttAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCCTTGAAGGGGTGCTCAcca                                                                 <  1:1004711/70‑2 (MQ=255)
     ttAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                >  1:2361716/1‑70 (MQ=255)
      tAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTc                                                               >  1:1644271/1‑70 (MQ=255)
       aGCTGGCCCAATACCTTGCTTTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCa                                                              <  1:2592577/70‑1 (MQ=255)
            gCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCa                                                         >  1:1559057/1‑70 (MQ=255)
               cAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:489782/1‑54 (MQ=255)
               cAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:1333212/1‑54 (MQ=255)
               cAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:3694209/1‑54 (MQ=255)
               cAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:3424106/1‑54 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                >  1:2794074/1‑58 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                >  1:27900/1‑58 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                >  1:2224525/1‑58 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCtt                                                                >  1:1160943/1‑58 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGt                                                   <  1:972033/71‑1 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGt                                                   <  1:3867535/71‑1 (MQ=255)
                 aTACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGAACACCTTCAAAGCAAACCGt                                                   >  1:9546/1‑71 (MQ=255)
                  tACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCaaa                                                            >  1:3268810/1‑61 (MQ=255)
                  tACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCaaa                                                            >  1:1560961/1‑61 (MQ=255)
                  tACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGtt                                                  <  1:2361811/71‑1 (MQ=255)
                  tACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTGCaaa                                                            >  1:4367359/1‑61 (MQ=255)
                    ccTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:285831/1‑49 (MQ=255)
                    ccTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:468933/1‑49 (MQ=255)
                    ccTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGca                                                                      >  1:632926/1‑49 (MQ=255)
                    ccTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGttt                                                 <  1:4104419/70‑1 (MQ=255)
                    ccTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGttt                                                 <  1:4356264/70‑1 (MQ=255)
                    ccTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGttt                                                 <  1:2818805/70‑1 (MQ=255)
                     cTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTgg                                               <  1:974271/71‑1 (MQ=255)
                      ttGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGa                                              >  1:1575291/1‑71 (MQ=255)
                       tGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGAt                                             >  1:967282/1‑71 (MQ=255)
                        gCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGAt                                             <  1:4264411/70‑1 (MQ=255)
                         cTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGAt                                             <  1:653697/69‑1 (MQ=255)
                         cTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAAACGTTTGGATCt                                           <  1:2833660/71‑1 (MQ=255)
                            tACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCaaa                                                       >  1:693519/1‑56 (MQ=255)
                            tACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTc                                        <  1:1227752/71‑1 (MQ=255)
                            tACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGAACACCTTCAAAGCaaa                                                       >  1:3864433/1‑56 (MQ=255)
                              cAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCt                                           >  1:517891/1‑66 (MQ=255)
                                aGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACt                                    >  1:1952357/1‑71 (MQ=255)
                                aGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACt                                    >  1:3761875/1‑71 (MQ=255)
                                 gATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACtt                                   >  1:3189020/1‑71 (MQ=255)
                                   tCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACttt                                  >  1:3608947/1‑70 (MQ=255)
                                    cTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGtt                                                  >  1:4412586/1‑53 (MQ=255)
                                    cTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGtt                                                  >  1:3523673/1‑53 (MQ=255)
                                    cTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGtt                                                  >  1:4183559/1‑53 (MQ=255)
                                      aTTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTAc                               >  1:3647898/1‑70 (MQ=255)
                                       ttAACGCTGGCGCTTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTAc                               <  1:1094292/69‑1 (MQ=255)
                                       ttAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCa                             >  1:3595270/1‑71 (MQ=255)
                                        tAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCaaa                                                       >  1:1329659/1‑44 (MQ=255)
                                        tAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCAt                            >  1:2843573/1‑71 (MQ=255)
                                        tAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCAt                            >  1:1749014/1‑71 (MQ=255)
                                         aaCGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGa                                              >  1:1931986/1‑52 (MQ=255)
                                          aCGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCAt                            <  1:3991224/69‑1 (MQ=255)
                                                gcgcGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAg                      <  1:822429/69‑1 (MQ=255)
                                                       gtTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACt             >  1:1824472/1‑71 (MQ=255)
                                                       gtTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACt             >  1:2842633/1‑71 (MQ=255)
                                                       gtTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACt             >  1:915837/1‑71 (MQ=255)
                                                         tGAAGGTGTGCTCACCATCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgt           >  1:4255711/1‑71 (MQ=38)
                                                         tGAAGGTGTGCTCACCATCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgt           >  1:654557/1‑71 (MQ=38)
                                                          gATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgtg          >  1:751953/1‑71 (MQ=255)
                                                          gATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgtg          >  1:536176/1‑71 (MQ=255)
                                                          gATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgtg          >  1:1623969/1‑71 (MQ=255)
                                                          gATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtg            >  1:3219112/1‑69 (MQ=255)
                                                               tGAGCACACCTTCAAAGTAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGt     >  1:2116741/1‑71 (MQ=255)
                                                               tGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGt     >  1:3775857/1‑71 (MQ=255)
                                                               tGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGt     <  1:812309/71‑1 (MQ=255)
                                                                gAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGtt    >  1:2624371/1‑71 (MQ=255)
                                                                gAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGtt    >  1:2585569/1‑71 (MQ=255)
                                                                 aGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATa               >  1:135101/1‑59 (MQ=255)
                                                                 aGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgt           >  1:3250589/1‑63 (MQ=255)
                                                                 aGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACtgt           >  1:1176891/1‑63 (MQ=255)
                                                                  gCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGTTcc  >  1:505363/1‑71 (MQ=255)
                                                                  gCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGTTcc  <  1:4018454/71‑1 (MQ=255)
                                                                   cacaCCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGt     <  1:2961896/67‑1 (MQ=255)
                                                                    acacCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATAc              <  1:1480301/57‑1 (MQ=255)
                                                                    acacCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATAc              <  1:839483/57‑1 (MQ=255)
                                                                    acacCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATAc              <  1:4148996/57‑1 (MQ=255)
                                                                    |                                                                    
TCGATTTAGCTGGCCCAATACCTTGCTGTACAAGATCTATTAACGCTGGCGCGTCGTTGATGGTGAGCACACCTTCAAAGCAAACCGTTTGGATCTTTCCACTTTTACCATCACCAGAAAAATACTGTGGCCGTTCC  >  Exported/1849305‑1849441

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: ‑