Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 386,224 T→C 8.7% F478F (TTT→TTC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported386,2240TC 8.7% 219.2 / 9.1 80F478F (TTT→TTCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base T (38/35);  new base C (3/4);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAG  >  Exported/386154‑386287
                                                                      |                                                               
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)
    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)
        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)
               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)
                            cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                           >  1:1570879/1‑49 (MQ=255)
                            cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                           >  1:5668714/1‑49 (MQ=255)
                            cATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCa                                                           >  1:1252992/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:4406319/71‑1 (MQ=255)
                                atggatgTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGc                                 <  1:3084162/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:293709/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)
                                         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:4382262/71‑1 (MQ=255)
                                          cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                       <  1:3870259/71‑1 (MQ=255)
                                          cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                       <  1:5282299/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: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)
                                                             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)
                                                                      |                                                               
CCTACATTCTGCACGGCGCGCTGGGCGGCATGATGGATGTGGCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAG  >  Exported/386154‑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: