Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,249 G→C 8.7% V487L (GTG→CTG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,2490GC 8.7% 234.3 / 8.4 80V487L (GTG→CTG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (33/40);  new base C (3/4);  total (36/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 = 1.00e+00

ATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCC  >  minE/386183‑386315
                                                                  |                                                                  
atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                                                >  1:3274472/1‑71 (MQ=255)
atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                                                >  1:1404732/1‑71 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcggcgg                                                                    >  1:1608238/1‑64 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                                              >  1:5325828/1‑70 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                                              >  1:5324654/1‑70 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                                              <  1:1738957/70‑1 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                             <  1:3084162/71‑1 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                             <  1:4406319/71‑1 (MQ=255)
   atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                             <  1:4526423/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:806370/71‑1 (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)
            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)
               ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                                             <  1:3724063/59‑1 (MQ=255)
               ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgc                                                 <  1:472205/71‑1 (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: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)
                                    ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGa                                                                >  1:5313221/1‑35 (MQ=255)
                                    ttgtttgAAGCACTGTTTATTCTGACGGCGGTGGa                                                                >  1:1346262/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: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)
                                                                  |                                                                  
ATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCC  >  minE/386183‑386315

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: