Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,236 T→C 8.8% F482F (TTT→TTC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,2360TC 8.8% 216.4 / 9.4 79F482F (TTT→TTCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (33/39);  new base C (3/4);  total (36/43)
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.98e-01

ACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGT  >  minE/386166‑386306
                                                                      |                                                                      
aCGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGttt                                                                        <  1:1198125/71‑1 (MQ=255)
 cGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                                       >  1:1146638/1‑71 (MQ=255)
 cGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                                       <  1:950548/71‑1 (MQ=255)
   gcgcgcTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                                       >  1:2366779/1‑69 (MQ=255)
        cTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAc                                                                >  1:4307277/1‑71 (MQ=255)
         tGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAc                                                                >  1:5597332/1‑70 (MQ=255)
          gggCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcgg                                                              >  1:3454238/1‑71 (MQ=255)
                 atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                                       >  1:1404732/1‑71 (MQ=255)
                 atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                                       >  1:3274472/1‑71 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcggcgg                                                           >  1:1608238/1‑64 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                                     >  1:5324654/1‑70 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                                     <  1:1738957/70‑1 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                                     >  1:5325828/1‑70 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                    <  1:4526423/71‑1 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                    <  1:4406319/71‑1 (MQ=255)
                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                    <  1:3084162/71‑1 (MQ=255)
                       gatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                                       >  1:4132793/1‑65 (MQ=255)
                       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:806370/71‑1 (MQ=255)
                             gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                           <  1:2668341/71‑1 (MQ=255)
                             gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                           <  1:5476324/71‑1 (MQ=255)
                             gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgc                                            <  1:4521844/70‑1 (MQ=255)
                              gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTAcccg                                           >  1:3733837/3‑70 (MQ=255)
                              gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGt                                          <  1:3008848/69‑1 (MQ=255)
                              gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGt                                          <  1:305927/69‑1 (MQ=255)
                              gcgtCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGt                                                <  1:5138441/62‑1 (MQ=255)
                              cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                           <  1:3595086/70‑1 (MQ=255)
                              cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                          <  1:5282299/71‑1 (MQ=255)
                              cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                          <  1:4382262/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)
                                ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGttt                                                                        <  1:4219324/39‑1 (MQ=255)
                                ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTAt                                                                      <  1:3386721/41‑1 (MQ=255)
                                ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                    <  1:3724063/59‑1 (MQ=255)
                                ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgc                                        <  1:472205/71‑1 (MQ=255)
                                  cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTAt                                                                      >  1:656532/1‑39 (MQ=255)
                                  cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcgg                                                              >  1:558162/1‑47 (MQ=255)
                                  cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCa                                                   >  1:5628888/1‑58 (MQ=255)
                                  cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                       <  1:330016/70‑1 (MQ=255)
                                   tGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGc                                                    <  1:3908714/56‑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:3783402/1‑52 (MQ=255)
                                         cACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAg                                                  >  1:5262176/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:1514813/1‑57 (MQ=255)
                                               gCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                       >  1:2355952/1‑57 (MQ=255)
                                               gCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                       >  1:3073553/1‑57 (MQ=255)
                                                ccATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTg                        <  1:5128654/71‑1 (MQ=255)
                                                 cATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGc                       <  1:444402/71‑1 (MQ=255)
                                                 cATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGc                       <  1:1934363/71‑1 (MQ=255)
                                                   ttttGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAg                     <  1:3063722/71‑1 (MQ=255)
                                                    tttgtttgAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAg                     >  1:1372550/1‑70 (MQ=255)
                                                     ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGa                                                       >  1:1346262/1‑35 (MQ=255)
                                                     ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGa                                                       >  1:5313221/1‑35 (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:4928386/49‑1 (MQ=255)
                                                           gAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcgc                                   <  1:4128685/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:4421787/1‑38 (MQ=255)
                                                                   gTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgctg                                      >  1:1157827/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)
                                                                      |                                                                      
ACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGT  >  minE/386166‑386306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: