Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,257 A→G 7.8% A489A (GCA→GCG cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,2570AG 7.8% 216.2 / 6.8 77A489A (GCA→GCGcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (30/41);  new base G (3/3);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

GATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACC  >  minE/386189‑386326
                                                                    |                                                                     
gatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAgg                                                                     <  1:1297316/71‑1 (MQ=255)
 atgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGTTGCAGGt                                                                    <  1:33828/71‑1 (MQ=255)
 atgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAgg                                                                     >  1:5211871/1‑70 (MQ=255)
      gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                               <  1:3166406/71‑1 (MQ=255)
      gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                               <  1:293709/71‑1 (MQ=255)
      gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                               <  1:5476324/71‑1 (MQ=255)
      gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                               <  1:2668341/71‑1 (MQ=255)
      gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                               <  1:806370/71‑1 (MQ=255)
      gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgc                                                                <  1:4521844/70‑1 (MQ=255)
       gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTAcccg                                                               >  1:3733837/3‑70 (MQ=255)
       gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGt                                                              <  1:305927/69‑1 (MQ=255)
       gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGt                                                              <  1:3008848/69‑1 (MQ=255)
       gcgtCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGt                                                                    <  1:5138441/62‑1 (MQ=255)
       cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                               <  1:3595086/70‑1 (MQ=255)
       cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                                              <  1:4382262/71‑1 (MQ=255)
       cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                                              <  1:5282299/71‑1 (MQ=255)
       cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                                              <  1:3870259/71‑1 (MQ=255)
        gTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtg                                                             <  1:731110/71‑1 (MQ=255)
        gTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtg                                                             <  1:3025731/71‑1 (MQ=255)
         ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgc                                                            <  1:472205/71‑1 (MQ=255)
           cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCa                                                                       >  1:5628888/1‑58 (MQ=255)
           cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                                           <  1:330016/70‑1 (MQ=255)
            tGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGttctgg                                                         >  1:3426748/1‑70 (MQ=255)
            tGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGttctgg                                                         >  1:4394724/1‑70 (MQ=255)
                aTCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcgcgc                                                     >  1:2687181/1‑71 (MQ=255)
                aTCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcg                                                        >  1:1395470/1‑68 (MQ=255)
                aTCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcg                                                        >  1:5136372/1‑68 (MQ=255)
                  cACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAg                                                                      >  1:5262176/1‑52 (MQ=255)
                  cACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAg                                                                      >  1:3783402/1‑52 (MQ=255)
                   aCTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCtt                                                   >  1:627540/1‑70 (MQ=255)
                     ttCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcg                                                                 <  1:4238556/54‑1 (MQ=255)
                     ttCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcg                                                                 <  1:714393/54‑1 (MQ=255)
                      tCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATg                                               >  1:1132488/1‑71 (MQ=255)
                       cGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTAt                                                <  1:5523042/69‑1 (MQ=255)
                        gCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                                           >  1:3073553/1‑57 (MQ=255)
                        gCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                                           >  1:2355952/1‑57 (MQ=255)
                        gCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                                           >  1:1514813/1‑57 (MQ=255)
                         ccATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTg                                            <  1:5128654/71‑1 (MQ=255)
                          cATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGc                                           <  1:1934363/71‑1 (MQ=255)
                          cATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGc                                           <  1:444402/71‑1 (MQ=255)
                            ttttGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAg                                         <  1:3063722/71‑1 (MQ=255)
                             tttgtttgAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAg                                         >  1:1372550/1‑70 (MQ=255)
                              ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTa                                                                   >  1:4652966/1‑43 (MQ=255)
                              ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgctgc                                                         >  1:5369582/1‑53 (MQ=255)
                              ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTAt                                                >  1:671389/1‑62 (MQ=255)
                                    gAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcgc                                                       <  1:4128685/49‑1 (MQ=255)
                                    gAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcgc                                                       <  1:4928386/49‑1 (MQ=255)
                                    gAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTa                                                 <  1:181016/55‑1 (MQ=255)
                                    gAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATctgc                                  <  1:2314553/70‑1 (MQ=255)
                                    gAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGATGCGCGCTTTa                                                 <  1:927550/55‑1 (MQ=255)
                                    gAAGCACTGTTTATTTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTa                                                 <  1:1710832/54‑1 (MQ=255)
                                     aaGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATctgctg                                <  1:4754531/71‑1 (MQ=255)
                                     aaGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATctgct                                 <  1:1764216/70‑1 (MQ=255)
                                           tGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtg                                                             <  1:5249229/36‑1 (MQ=255)
                                            gTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgctg                                                          >  1:1157827/1‑38 (MQ=255)
                                            gTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgctg                                                          >  1:4421787/1‑38 (MQ=255)
                                             tttATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATctgc                                  >  1:1737934/1‑61 (MQ=255)
                                             tttATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgtg                          <  1:3732176/69‑1 (MQ=255)
                                               tATTCTGACGGCGGTGGATGCAGGTACGCGtgctgc                                                         >  1:3637979/1‑36 (MQ=255)
                                               tATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgtggtg                       >  1:3963062/1‑70 (MQ=255)
                                               tATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGgtgt                      >  1:5088644/1‑71 (MQ=255)
                                                aTTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTAt                                                >  1:1715966/1‑44 (MQ=255)
                                                   cTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTcc                  >  1:2346000/1‑71 (MQ=255)
                                                      aCGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGgtgt                      >  1:456495/1‑64 (MQ=255)
                                                      aCGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGtc                     <  1:4305653/65‑1 (MQ=255)
                                                      aCGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTcc                  <  1:262788/68‑1 (MQ=255)
                                                      aCGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTgg               >  1:3769119/1‑71 (MQ=255)
                                                       cggcggTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgt                           <  1:3391405/58‑1 (MQ=255)
                                                       cggcggTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTcc                  <  1:269545/67‑1 (MQ=255)
                                                        ggcggTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGtct                    >  1:2776551/1‑64 (MQ=255)
                                                        ggcggTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGcc             >  1:1775749/1‑71 (MQ=255)
                                                             tGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGaa         <  1:4752994/70‑1 (MQ=255)
                                                               gATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCg      >  1:3893224/1‑71 (MQ=255)
                                                               gATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCg      <  1:2587683/71‑1 (MQ=255)
                                                                  gCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAc   >  1:3580561/1‑71 (MQ=255)
                                                                   cAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc  <  1:5419814/71‑1 (MQ=255)
                                                                   cAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc  <  1:4426125/71‑1 (MQ=255)
                                                                   cAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc  >  1:1694810/1‑71 (MQ=255)
                                                                    |                                                                     
GATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACC  >  minE/386189‑386326

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: