Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 386,263 G→C 5.6% T491T (ACG→ACC cstA → carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported386,2630GC 5.6% 218.8 / 2.4 71T491T (ACG→ACCcstAcarbon starvation protein
Reads supporting (aligned to +/- strand):  ref base G (32/35);  new base C (2/2);  total (34/37)
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.41e-01

GCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATT  >  Exported/386195‑386330
                                                                    |                                                                   
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:305927/69‑1 (MQ=255)
 gcttCTGGTATCACTTCGGTATTCTGTTCGAAGCCCTGTTCATCCTGACCGCGCTGGATGCGGGTACCCGt                                                                  <  1:3008848/69‑1 (MQ=255)
 cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcgcg                                                                   <  1:3595086/70‑1 (MQ=255)
 cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                                                  <  1:3870259/71‑1 (MQ=255)
 cGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGt                                                                  <  1:4382262/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)
   ttCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgc                                                                <  1:472205/71‑1 (MQ=255)
     cTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGtgct                                                               <  1:330016/70‑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)
             aCTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCtt                                                       >  1:627540/1‑70 (MQ=255)
               ttCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcg                                                                     <  1:714393/54‑1 (MQ=255)
               ttCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTAcg                                                                     <  1:4238556/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)
                        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)
                                                       tGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGaa             <  1:4752994/70‑1 (MQ=255)
                                                         gATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCg          >  1:3893224/1‑71 (MQ=255)
                                                         gATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCg          <  1:2587683/71‑1 (MQ=255)
                                                            gCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAc       >  1:3580561/1‑71 (MQ=255)
                                                             cAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc      <  1:5419814/71‑1 (MQ=255)
                                                             cAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc      >  1:1694810/1‑71 (MQ=255)
                                                             cAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc      <  1:4426125/71‑1 (MQ=255)
                                                               ggTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGGGTCTCCTGGCCTGAAGCGGACCGa    <  1:3930613/71‑1 (MQ=255)
                                                                gTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCg                                <  1:4928545/42‑1 (MQ=255)
                                                                 tACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgt                               <  1:4492896/42‑1 (MQ=255)
                                                                 tACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgt                               <  1:4861116/42‑1 (MQ=255)
                                                                 tACGCGGGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgt                               <  1:3472685/42‑1 (MQ=255)
                                                                  aCGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGgtgt                          >  1:777534/1‑46 (MQ=255)
                                                                  aCGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAg            >  1:1563152/1‑60 (MQ=255)
                                                                  aCGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGAcc      >  1:5655772/1‑66 (MQ=255)
                                                                  aCGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGAtt  >  1:2668712/1‑70 (MQ=255)
                                                                  aCGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGAtt  >  1:2691392/1‑70 (MQ=255)
                                                                  aCGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGAtt  >  1:1200959/1‑70 (MQ=255)
                                                                   cgcgTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCgtgg                             >  1:5437911/1‑42 (MQ=255)
                                                                   cgcgTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTTGTGTCTCCTGGCCTGa              <  1:1267999/57‑1 (MQ=255)
                                                                   cgcgTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGa              <  1:4541683/57‑1 (MQ=255)
                                                                   cgcgTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGa              <  1:1187200/57‑1 (MQ=255)
                                                                   cgcgTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGa              <  1:3940326/57‑1 (MQ=255)
                                                                    |                                                                   
GCGTTCTGGTATCACTTCGCCATTTTGTTTGAAGCACTGTTTATTCTGACGGCGGTGGATGCAGGTACGCGTGCTGCGCGCTTTATGTTGCAGGATCTGCTGGGCGTGGTGTCTCCTGGCCTGAAGCGGACCGATT  >  Exported/386195‑386330

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: