Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,572,708 A→C 34.5% intergenic (+1476/‑303) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,572,7080AC34.5% 117.7 / 171.1 167intergenic (+1476/‑303)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (62/47);  new base C (35/22);  total (98/69)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

GGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTA  >  Exported/2572641‑2572777
                                                                   |                                                                     
 gTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCtt                                                                   >  1:2416482/1‑71 (MQ=255)
  tGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  <  1:1357980/71‑1 (MQ=255)
    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  >  1:1142815/1‑69 (MQ=255)
    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  >  1:4964614/1‑69 (MQ=255)
    gCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                <  1:1422235/71‑1 (MQ=255)
      tAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGc                                                              >  1:49578/1‑71 (MQ=255)
        aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGcc                                                             <  1:3723592/70‑1 (MQ=255)
         tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                           <  1:2969901/71‑1 (MQ=255)
          cccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                          >  1:3999477/1‑71 (MQ=255)
           ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACATACCTTCCAGCCCaa                                                          >  1:2296532/1‑70 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  <  1:1032664/61‑1 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  <  1:2106058/61‑1 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  <  1:3296740/61‑1 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                        >  1:73009/1‑71 (MQ=255)
            cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                        >  1:5319246/1‑71 (MQ=255)
            cAAGGGGATCGATTCCAGAAATAGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                         >  1:55964/1‑70 (MQ=255)
             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTc                                                                  <  1:917203/60‑1 (MQ=255)
             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:1465919/1‑71 (MQ=255)
             aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:4834487/1‑71 (MQ=255)
               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:2091529/1‑69 (MQ=255)
               ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:3756246/1‑69 (MQ=255)
                  gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                           <  1:2036989/62‑1 (MQ=255)
                  gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGcc                                                  >  1:1086747/1‑71 (MQ=255)
                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCa                                                                >  1:5590047/1‑49 (MQ=255)
                          ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                           <  1:1882051/70‑1 (MQ=255)
                           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         <  1:2017498/71‑1 (MQ=255)
                           cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                         >  1:28989/1‑71 (MQ=255)
                             gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                        <  1:1170903/67‑1 (MQ=255)
                             gAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTtata                                       <  1:1731037/71‑1 (MQ=255)
                              aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                >  1:2153154/1‑61 (MQ=255)
                                aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                             <  1:4897259/62‑1 (MQ=255)
                                 tCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGatc                                                   <  1:2424029/55‑3 (MQ=255)
                                 ccGGATTCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                           <  1:4494340/62‑2 (MQ=255)
                                 ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                   >  1:1660448/2‑71 (MQ=255)
                                 ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                   >  1:1231848/2‑71 (MQ=255)
                                  cGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAg                                                               <  1:171590/42‑1 (MQ=255)
                                    ggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGcccc                                <  1:2673676/71‑1 (MQ=255)
                                     attCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCgg                                                      <  1:5664344/47‑1 (MQ=255)
                                     attCAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGccc                                 >  1:3231085/2‑69 (MQ=255)
                                      ttCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCctg                              >  1:4336450/1‑71 (MQ=255)
                                      ttCAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATAcc                                     >  1:1843433/1‑64 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAgt                                             <  1:5182467/56‑1 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                              >  1:1136220/1‑69 (MQ=255)
                                      ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                          <  1:1049083/43‑1 (MQ=255)
                                       tCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCctg                              >  1:27694/1‑70 (MQ=255)
                                       tCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCct                               >  1:4019555/1‑69 (MQ=255)
                                        cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                             <  1:5249456/70‑1 (MQ=255)
                                        cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGatc                                                   >  1:813282/1‑46 (MQ=255)
                                        cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                              >  1:3173499/1‑67 (MQ=255)
                                        cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                            >  1:3994771/1‑71 (MQ=255)
                                        aaaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                           >  1:711242/2‑56 (MQ=255)
                                         aaCGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgcaa                           >  1:5397850/1‑69 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTg                                                >  1:3002682/1‑50 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGCCCctgcag                           >  1:1521892/1‑69 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAAt                                                         >  1:4998908/1‑41 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCctgcag                           <  1:2167400/71‑3 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:2630580/1‑43 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                       >  1:671114/1‑43 (MQ=255)
                                         aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgctg                           >  1:3513883/1‑71 (MQ=255)
                                          aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                           >  1:3815513/1‑68 (MQ=255)
                                          aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                      >  1:610003/1‑59 (MQ=255)
                                           cGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                           <  1:1588314/53‑1 (MQ=255)
                                            ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                        >  1:5171066/1‑71 (MQ=255)
                                            ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                           <  1:241324/52‑1 (MQ=255)
                                            ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                      >  1:1156018/1‑57 (MQ=255)
                                             gATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGccc                                 <  1:175360/61‑1 (MQ=255)
                                             gATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                        >  1:5212955/1‑70 (MQ=255)
                                             gATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                        >  1:428937/1‑38 (MQ=255)
                                             gATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                        >  1:4633215/1‑38 (MQ=255)
                                              aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                       >  1:2729850/1‑70 (MQ=255)
                                               tACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGATATATCGCCCCTGCAGCgggtc                     <  1:2850444/71‑4 (MQ=255)
                                               tACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                           <  1:329436/65‑3 (MQ=255)
                                               tACTGATAAAAATTCCATCCCCACTTCAGCCCAATCGGATCCTGAg                                              >  1:4576623/1‑46 (MQ=255)
                                               tACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGa                                                     <  1:1103184/39‑1 (MQ=255)
                                                aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGa                                               >  1:240013/1‑44 (MQ=255)
                                                 cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGa                                               >  1:1606250/1‑43 (MQ=255)
                                                 cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                 >  1:4294827/1‑57 (MQ=255)
                                                  tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                 <  1:2534791/40‑1 (MQ=255)
                                                  tGATAAAAATTCCACCTACCTCCCAGCCCAATCGGAGCCt                                                 <  1:432037/40‑1 (MQ=255)
                                                   gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                        >  1:4979797/1‑64 (MQ=255)
                                                   gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                        >  1:1934447/1‑64 (MQ=255)
                                                    aTAAAAATTCCATCCCGCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                      >  1:3379299/1‑48 (MQ=255)
                                                    aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                      >  1:2503737/1‑48 (MQ=255)
                                                    aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                       >  1:3729263/1‑48 (MQ=255)
                                                    aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                         >  1:2750914/1‑46 (MQ=255)
                                                    aTAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                      >  1:4383217/1‑49 (MQ=255)
                                                      aaaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGatat                                        >  1:4746828/1‑45 (MQ=255)
                                                       aaaaTTCCATCACCCCTTCAGCCCAATCGGATCCTGAGTGAtata                                       >  1:5468939/1‑45 (MQ=255)
                                                       aaaaTTCCACCTACCTTCCAGCCCAATCGGAGCATGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTAcc             >  1:4425874/1‑71 (MQ=255)
                                                         aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAcctc                                   <  1:4132582/47‑3 (MQ=255)
                                                             ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa        <  1:5065269/70‑1 (MQ=255)
                                                             ccACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataa        >  1:4211728/1‑70 (MQ=255)
                                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat       <  1:2640986/70‑1 (MQ=255)
                                                              cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat       <  1:967637/70‑1 (MQ=255)
                                                               acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                          <  1:309814/50‑1 (MQ=255)
                                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                          <  1:399991/48‑1 (MQ=255)
                                                                 ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCgg                        >  1:2143794/1‑50 (MQ=255)
                                                                  tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                             <  1:3548074/44‑1 (MQ=255)
                                                                  tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                             <  1:5145287/44‑1 (MQ=255)
                                                                  tacctTCCAGCCCAATCGGAGCCTGAGAGTTATACCGCCCCAGCTGCGGTTAAAGTTACCGAAAATAGAta  >  1:2432791/1‑71 (MQ=255)
                                                                   |                                                                     
GGTGGCTAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTA  >  Exported/2572641‑2572777

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: