Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,572,717 C→T 9.0% intergenic (+1485/‑294) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,572,7170CT 9.0% 500.4 / 26.8 178intergenic (+1485/‑294)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (87/75);  new base T (13/3);  total (100/78)
Fisher's exact test for biased strand distribution p-value = 3.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 187 total aligned reads displayed.

TAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCC  >  Exported/2572647‑2572787
                                                                      |                                                                      
  aTCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                                      <  1:5210802/71‑1 (MQ=255)
   tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGccc                                                                      >  1:137772/1‑70 (MQ=255)
   tCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                                     <  1:2969901/71‑1 (MQ=255)
    cccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                                    >  1:3999477/1‑71 (MQ=255)
    cccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                                    >  1:4792942/1‑71 (MQ=255)
     ccAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACATACCTTCCAGCCCaa                                                                    >  1:2296532/1‑70 (MQ=255)
      cAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                                  >  1:73009/1‑71 (MQ=255)
      cAAGGGGATCGATTCCAGAAATAGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   >  1:55964/1‑70 (MQ=255)
       aaGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                                 >  1:1465919/1‑71 (MQ=255)
         ggggATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                                 >  1:2091529/1‑69 (MQ=255)
            gATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                                     <  1:2036989/62‑1 (MQ=255)
             aTCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCt                                                           <  1:423499/71‑1 (MQ=255)
                gATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGa                                                         <  1:1099265/70‑1 (MQ=255)
                    ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     <  1:4341378/70‑1 (MQ=255)
                    ccAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     <  1:1882051/70‑1 (MQ=255)
                     cAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   <  1:4871295/71‑1 (MQ=255)
                       gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                                  <  1:1170903/67‑1 (MQ=255)
                       gagaTCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGatat                                                  <  1:3519381/67‑1 (MQ=255)
                        agaTCGGATTCAACGGATACTGATAAAAATTCCATCCTCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                                 >  1:4143937/3‑70 (MQ=255)
                        aaaTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTg                                                          >  1:2153154/1‑61 (MQ=255)
                          aTCTGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgt                                                       <  1:877913/62‑1 (MQ=255)
                          aTCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCg                                              <  1:1734648/71‑1 (MQ=255)
                           tCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCa                                                                     <  1:4050751/47‑1 (MQ=255)
                           tCGGATTCAACGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGatc                                                             <  1:2424029/55‑3 (MQ=255)
                           ccGGATTCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                     <  1:4494340/62‑2 (MQ=255)
                           ccGGATTCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGc                                             >  1:1660448/2‑71 (MQ=255)
                              ggTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGcccc                                          <  1:2673676/71‑1 (MQ=255)
                               attCAGCGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCgg                                                                <  1:5664344/47‑1 (MQ=255)
                               attCAGCGGATACGCATACAGACTCCATCCCCCCTCCAGCCCTATCGGGTCCTGAGTGATATACCGCCCct                                         <  1:2205840/70‑1 (MQ=255)
                               attCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGat                                                              >  1:320179/2‑49 (MQ=255)
                                ttCAACGTATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   >  1:2577749/1‑60 (MQ=255)
                                ttCAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATAcc                                               >  1:1843433/1‑64 (MQ=255)
                                ttCAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                        >  1:1136220/1‑69 (MQ=255)
                                ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                                    <  1:2808058/43‑1 (MQ=255)
                                ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCaa                                                                    <  1:1049083/43‑1 (MQ=255)
                                ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   <  1:4411233/60‑1 (MQ=255)
                                ttCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   >  1:2413257/1‑60 (MQ=255)
                                  cAACGGATACTGATACATATTCCATCCACCTCTCAACCCAATCGGGTCCGGGGTGATATACCGCCCctgc                                       <  1:5249456/70‑1 (MQ=255)
                                  cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGatc                                                             >  1:813282/1‑46 (MQ=255)
                                  cAACGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGccccag                                        >  1:3173499/1‑67 (MQ=255)
                                  cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   >  1:3854995/1‑42 (MQ=255)
                                  cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   >  1:4909787/1‑42 (MQ=255)
                                  cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTta                                                   >  1:1947126/1‑58 (MQ=255)
                                  cAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgct                                      >  1:3994771/1‑71 (MQ=255)
                                  aaaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     >  1:711242/2‑56 (MQ=255)
                                   aaCGGATACTGATAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTg                                                          >  1:3002682/1‑50 (MQ=255)
                                   aaCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAAt                                                                   >  1:4998908/1‑41 (MQ=255)
                                   aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                                 >  1:671114/1‑43 (MQ=255)
                                   aaCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCg                                                                 >  1:3729997/1‑43 (MQ=255)
                                    aCGGATACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCcagcag                                     >  1:554642/1‑68 (MQ=255)
                                    aCGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                >  1:610003/1‑59 (MQ=255)
                                      ggATACTGATAAAAATTCCATCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCTCCCCTGCAGCggg                                 >  1:2878568/1‑70 (MQ=255)
                                      ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAgtga                                                     <  1:5222836/52‑2 (MQ=255)
                                      ggATACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATACCGCCCCTGCAGCgg                                  >  1:5171066/1‑71 (MQ=255)
                                      ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAgtgt                                                     <  1:241324/52‑1 (MQ=255)
                                      ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                >  1:3432235/1‑57 (MQ=255)
                                      ggATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATAc                                                >  1:1156018/1‑57 (MQ=255)
                                       gATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                                  >  1:428937/1‑38 (MQ=255)
                                       gATACTGATAAAAATTCCACCTACCTTCCAGCCCAATc                                                                  >  1:4480579/1‑38 (MQ=255)
                                        aTACTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGa                                                               >  1:3232991/1‑40 (MQ=255)
                                        aTACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGt                                 >  1:2729850/1‑70 (MQ=255)
                                         tACTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGAg                                                        >  1:2754358/1‑46 (MQ=255)
                                         tACTGATAAAAATTCCATCCCCACTTCAGCCCAATCGGATCCTGAg                                                        >  1:4576623/1‑46 (MQ=255)
                                         tACTGATAAAAATTCCACCTACCTTCCAGCCCAAt                                                                   <  1:5058875/35‑1 (MQ=255)
                                          aCTGATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGa                                                         >  1:240013/1‑44 (MQ=255)
                                           cTGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGat                                                              <  1:1049524/38‑2 (MQ=16)
                                           cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGa                                                         >  1:1606250/1‑43 (MQ=255)
                                           cTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGccc                                           >  1:4294827/1‑57 (MQ=255)
                                            tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                           <  1:3630668/40‑1 (MQ=255)
                                            tGATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                           <  1:2519766/40‑1 (MQ=255)
                                            tGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGtt                                <  1:4968894/67‑1 (MQ=255)
                                             gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                  >  1:4979797/1‑64 (MQ=255)
                                             gATAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCAGCAGCgg                                  >  1:1934447/1‑64 (MQ=255)
                                             gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                           <  1:1263718/39‑1 (MQ=16)
                                             gATAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCt                                                           <  1:5332975/39‑1 (MQ=16)
                                              aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtatat                                                >  1:1130632/1‑48 (MQ=255)
                                              aTAAAAATTCCATCCCCCTTTCAGCCCAATCGGATCCTGAGTGAtata                                                 >  1:3729263/1‑48 (MQ=255)
                                              aTAAAAATTCCATCCCCCCTTCAGCCCAATCGGATCCTGa                                                         >  1:4320881/1‑40 (MQ=255)
                                               tAAAAATTCCATCCCCCCTTCAGTCCAATCGGATCCTGAGTGATATACCGCCCCagcagc                                    <  1:1807372/60‑1 (MQ=255)
                                                aaaaaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGatat                                                  >  1:4746828/1‑45 (MQ=255)
                                                 aaaaTTCCATCACCCCTTCAGCCCAATCGGATCCTGAGTGAtata                                                 >  1:5468939/1‑45 (MQ=255)
                                                   aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGAta                                                   >  1:1699164/1‑41 (MQ=255)
                                                   aaTTCCATCCCCCCTTCAGCCCAATCGGATCCTGAGTGATATAcctc                                             <  1:4132582/47‑3 (MQ=255)
                                                        cACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGataat                 <  1:967637/70‑1 (MQ=255)
                                                         acctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                    <  1:309814/50‑1 (MQ=255)
                                                          cctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt              <  1:2116581/71‑1 (MQ=255)
                                                           ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCtgctgc                                    <  1:399991/48‑1 (MQ=255)
                                                           ctacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGt              >  1:2631853/1‑70 (MQ=255)
                                                            tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                       <  1:5145287/44‑1 (MQ=255)
                                                            tacctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                       <  1:158107/44‑1 (MQ=255)
                                                             acctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta            >  1:2439551/1‑70 (MQ=255)
                                                              cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGa                     <  1:3796863/60‑1 (MQ=255)
                                                              cctTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTta            <  1:5310379/69‑1 (MQ=255)
                                                                  tcAGCCCAATCGGATCCTGAGTGATATACCGcccc                                          <  1:3332834/34‑1 (MQ=37)
                                                                  ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCctgc                                       >  1:274979/1‑38 (MQ=255)
                                                                  ccAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTTCTGCGGTTAATGTTACCGATAATAGTTATAGt        <  1:5113808/69‑1 (MQ=255)
                                                                   cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTa       >  1:2144679/1‑69 (MQ=255)
                                                                   cAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTAc      <  1:2124118/70‑1 (MQ=255)
                                                                      ccccATCGGGTCCTGAGTAATATACCGCCCCTGCAACGGATCGTAGTACCGGTGACGGTTATAGTACAGcc  <  1:1764035/71‑1 (MQ=255)
                                                                      cccAATCGGGTCCGGGGTGATATACCGCCCCTGCAACGGATCGTAGTGCCGGTGACGGTTATAGTACAGcc  >  1:64164/1‑71 (MQ=255)
                                                                      |                                                                      
TAATCCCAAGGGGATCGATTCCAGAAATCGGGTTCAACGGATACTGATAAAAATTCCACCTACCTTCCAGCCCAATCGGAGCCTGAGTGTTATACCGCCCCTGCTGCGGTTAATGTTACCGATAATAGTTATAGTACAGCC  >  Exported/2572647‑2572787

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: