Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 386,219 T→C 8.8% L477L (TTG→CTG)  cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE386,2190TC 8.8% 215.5 / 9.4 79L477L (TTG→CTG) cstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base C (3/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 110 total aligned reads displayed.

ATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAG  >  minE/386150‑386287
                                                                     |                                                                    
aTGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCAtttt                                                                     <  1:5212071/71‑1 (MQ=255)
 tGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                    <  1:1016504/71‑1 (MQ=255)
 tGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttg                                                                    <  1:3216677/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:309098/1‑71 (MQ=255)
  ggCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                   >  1:1094609/1‑71 (MQ=255)
  ggCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                   >  1:3278328/1‑71 (MQ=255)
    ccTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgttt                                                                 >  1:4796762/1‑71 (MQ=255)
    ccTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgttt                                                                 >  1:1001011/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:4974094/1‑71 (MQ=255)
        cATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAg                                                             >  1:1375286/1‑71 (MQ=255)
          ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                                <  1:1762238/66‑1 (MQ=255)
          ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                                <  1:664898/66‑1 (MQ=255)
          ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                           >  1:5432269/1‑71 (MQ=255)
          ttCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                           <  1:2568366/71‑1 (MQ=255)
           tCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgtttg                                                                <  1:5065524/65‑1 (MQ=255)
           tCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                   >  1:728931/1‑62 (MQ=255)
           tCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATtttgt                                                                   >  1:3112648/1‑62 (MQ=255)
            cTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACt                                                         >  1:4523335/1‑71 (MQ=255)
               cACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACt                                                         <  1:463774/68‑1 (MQ=255)
                aCGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGttt                                                     <  1:1198125/71‑1 (MQ=255)
                 cGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                    <  1:950548/71‑1 (MQ=255)
                 cGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTa                                                    >  1:1146638/1‑71 (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:1404732/1‑71 (MQ=255)
                                 atgatgGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGa                                    >  1:3274472/1‑71 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGa                                                               >  1:2387878/1‑41 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGa                                                               >  1:1562947/1‑41 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCAc                                                          >  1:3907023/1‑46 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGAcggcgg                                        >  1:1608238/1‑64 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                  >  1:5324654/1‑70 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                  >  1:5325828/1‑70 (MQ=255)
                                    atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATg                                  <  1:1738957/70‑1 (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:806370/71‑1 (MQ=255)
                                             gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                        <  1:5476324/71‑1 (MQ=255)
                                             gCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                        <  1:3166406/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:3870259/71‑1 (MQ=255)
                                              cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                       <  1:4382262/71‑1 (MQ=255)
                                               gTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtg                      <  1:3025731/71‑1 (MQ=255)
                                               gTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtg                      <  1:731110/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)
                                                  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:4394724/1‑70 (MQ=255)
                                                   tGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGttctgg                  >  1:3426748/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: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:3073553/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)
                                                                   ttttGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAg  <  1:3063722/71‑1 (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)
                                                                     |                                                                    
ATGGCCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAG  >  minE/386150‑386287

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: