Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,215 C→T 6.8% A475A (GCC→GCT) ‡ cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,2150CT 6.8% 301.8 / 7.8 102A475A (GCC→GCT) ‡cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base C (52/43);  new base T (3/4);  total (55/47)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 105 total aligned reads displayed.

GGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGC  >  minE/386146‑386285
                                                                     |                                                                      
ggggATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCa                                                                       >  1:1751017/1‑71 (MQ=255)
     tGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                  <  1:1016504/71‑1 (MQ=255)
     tGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                  <  1:4429443/71‑1 (MQ=255)
      ggCCTACATTCTGCACGGCGCGCTGGGCGGCTTGATGGATGTGGCGTTCTGGTATCACTTCGCCAtttt                                                                   >  1:4882563/1‑69 (MQ=255)
      ggCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCAtttt                                                                   >  1:2253432/1‑69 (MQ=255)
      ggCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                 >  1:1094609/1‑71 (MQ=255)
      ggCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                 >  1:309098/1‑71 (MQ=255)
      ggCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                 >  1:3278328/1‑71 (MQ=255)
        ccTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgttt                                                               >  1:1001011/1‑71 (MQ=255)
        ccTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgttt                                                               >  1:4796762/1‑71 (MQ=255)
         cTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                              <  1:240902/71‑1 (MQ=255)
         cTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgttt                                                               >  1:3469030/1‑70 (MQ=255)
          tACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                  <  1:3400997/66‑1 (MQ=255)
            cATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAg                                                           >  1:1375286/1‑71 (MQ=255)
            cATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAg                                                           >  1:4974094/1‑71 (MQ=255)
             aTTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGc                                                          <  1:5648531/71‑1 (MQ=255)
              ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                              <  1:664898/66‑1 (MQ=255)
              ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                              <  1:1762238/66‑1 (MQ=255)
              ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                         <  1:2568366/71‑1 (MQ=255)
              ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                         >  1:5432269/1‑71 (MQ=255)
               tCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                              <  1:5065524/65‑1 (MQ=255)
               tCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                 >  1:3112648/1‑62 (MQ=255)
               tCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                 >  1:728931/1‑62 (MQ=255)
                cTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgttt                                                               >  1:4358461/1‑63 (MQ=255)
                cTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACt                                                       >  1:4523335/1‑71 (MQ=255)
                   cACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACt                                                       <  1:463774/68‑1 (MQ=255)
                   cACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGtt                                                    >  1:3220209/1‑71 (MQ=255)
                    aCGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGttt                                                   <  1:1198125/71‑1 (MQ=255)
                     cGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                  <  1:950548/71‑1 (MQ=255)
                     cGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                  >  1:1146638/1‑71 (MQ=255)
                       gcgcgcTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCAt                                                                      >  1:2373590/1‑49 (MQ=255)
                       gcgcgcTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                  >  1:2366779/1‑69 (MQ=255)
                          cgcTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                  <  1:2168707/50‑1 (MQ=255)
                          cgcTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                  <  1:2303913/50‑1 (MQ=255)
                            cTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAc                                           >  1:4307277/1‑71 (MQ=255)
                             tGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAc                                           >  1:5597332/1‑70 (MQ=255)
                              gggCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcgg                                         >  1:3454238/1‑71 (MQ=255)
                                    cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                  <  1:4123589/40‑1 (MQ=255)
                                    cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                         >  1:5668714/1‑49 (MQ=255)
                                    cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                         >  1:1252992/1‑49 (MQ=255)
                                    cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                         >  1:1570879/1‑49 (MQ=255)
                                     atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                  >  1:3274472/1‑71 (MQ=255)
                                     atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                  >  1:1404732/1‑71 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGa                                                             >  1:1562947/1‑41 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGa                                                             >  1:2387878/1‑41 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCAc                                                        >  1:3907023/1‑46 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcggcgg                                      >  1:1608238/1‑64 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                <  1:1738957/70‑1 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                >  1:5324654/1‑70 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                >  1:5325828/1‑70 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                               <  1:4526423/71‑1 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                               <  1:3084162/71‑1 (MQ=255)
                                        atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                               <  1:4406319/71‑1 (MQ=255)
                                           gatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGaa                                                            >  1:1079282/1‑39 (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:2668341/71‑1 (MQ=255)
                                                 gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                      <  1:293709/71‑1 (MQ=255)
                                                 gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                      <  1:3166406/71‑1 (MQ=255)
                                                 gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                      <  1:806370/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: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)
                                                    ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACt                                                       >  1:182072/1‑35 (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:5136372/1‑68 (MQ=255)
                                                           aTCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTgcg               >  1:1395470/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: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:1514813/1‑57 (MQ=255)
                                                                   gCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                  >  1:2355952/1‑57 (MQ=255)
                                                                    ccATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTg   <  1:5128654/71‑1 (MQ=255)
                                                                     cATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGc  <  1:444402/71‑1 (MQ=255)
                                                                     |                                                                      
GGGGATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGC  >  minE/386146‑386285

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: