Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,982,587 T→C 11.8% G97G (GGT→GGC creB → DNA‑binding response regulator in two‑component regulatory system with CreC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,982,5870TC11.8% 134.1 / 10.4 51G97G (GGT→GGCcreBDNA‑binding response regulator in two‑component regulatory system with CreC
Reads supporting (aligned to +/- strand):  ref base T (23/22);  new base C (3/3);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

CGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTAC  >  Exported/2982517‑2982651
                                                                      |                                                                
cGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGt                                                                  >  1:2898647/1‑71 (MQ=255)
cGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGt                                                                  <  1:2048723/71‑1 (MQ=255)
 ggCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTg                                                                 >  1:1165829/1‑71 (MQ=255)
   cGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCt                                                               <  1:2129594/71‑1 (MQ=255)
   cGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCt                                                               <  1:293632/71‑1 (MQ=255)
       aCCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTgac                                                            <  1:3837497/70‑1 (MQ=255)
           gTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTa                                                       >  1:533287/1‑71 (MQ=255)
           gTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTa                                                       >  1:394565/1‑71 (MQ=255)
              cTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGt                                                    >  1:580358/1‑71 (MQ=255)
                   ccTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCt                                                >  1:2158314/1‑70 (MQ=255)
                    cTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTa                                               <  1:1618463/70‑1 (MQ=255)
                     tGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTaaa                                             >  1:1618438/1‑71 (MQ=255)
                       aCGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAAcc                                           >  1:1353644/1‑71 (MQ=255)
                       aCGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAAcc                                           >  1:3190401/1‑71 (MQ=255)
                        cGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                          <  1:3675534/71‑1 (MQ=255)
                        cGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                          <  1:3293980/71‑1 (MQ=255)
                        cGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                          <  1:762847/71‑1 (MQ=255)
                            ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                      <  1:3583020/71‑1 (MQ=255)
                            ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                      <  1:899790/71‑1 (MQ=255)
                            ccGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGtttt                                      <  1:1206369/71‑1 (MQ=255)
                             cGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTc                                     <  1:1414884/71‑1 (MQ=255)
                              gAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCa                                    >  1:2839562/1‑71 (MQ=255)
                              gAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCa                                    >  1:1596080/1‑71 (MQ=255)
                              gAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCa                                    >  1:3594626/1‑71 (MQ=255)
                                    aaGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCcgc                              >  1:3118424/1‑71 (MQ=255)
                                    aaGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCcgc                              <  1:3198395/71‑1 (MQ=255)
                                      gagGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCcgcg                             >  1:3387739/1‑70 (MQ=255)
                                         gTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCg                                          <  1:1830333/54‑1 (MQ=255)
                                               cGCCTGCTTGGGCTGGAAATTGGTGATGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTgcgc                   >  1:328130/1‑71 (MQ=255)
                                                    gCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAggg              >  1:2887198/1‑71 (MQ=255)
                                                    gCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAggg              >  1:209167/1‑71 (MQ=255)
                                                     cTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAggg              >  1:1033054/1‑70 (MQ=255)
                                                     cTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAggg              >  1:182160/1‑70 (MQ=255)
                                                     cTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGt             <  1:707747/71‑1 (MQ=255)
                                                            tGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAgtg                        <  1:2341993/53‑1 (MQ=255)
                                                                aaTTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCa        >  1:300652/1‑65 (MQ=255)
                                                                aaTTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCCCCt     >  1:2911913/1‑68 (MQ=255)
                                                                aaTTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTAc  >  1:1190828/1‑71 (MQ=255)
                                                                 aTTGGCGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAgtg                        <  1:393676/48‑1 (MQ=38)
                                                                 aTTGGCGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAgtg                        <  1:3751695/48‑1 (MQ=38)
                                                                 aTTGGCGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAgtg                        <  1:3723703/48‑1 (MQ=38)
                                                                  ttGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGaa                           <  1:557895/44‑1 (MQ=255)
                                                                  ttGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGaa                           <  1:613414/44‑1 (MQ=255)
                                                                  ttGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGaa                           <  1:863050/44‑1 (MQ=255)
                                                                   tGGCGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGtgtg                      >  1:433756/1‑48 (MQ=38)
                                                                   tGGCGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGtgtg                      >  1:1267850/1‑48 (MQ=38)
                                                                   tGGCGCTGACGACTATGTGGCTAAACCGTTTTCACCCCGCGAAGtgtg                      >  1:334513/1‑48 (MQ=38)
                                                                    ggTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAg                          >  1:3819829/1‑43 (MQ=255)
                                                                    ggTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAg                          >  1:2833020/1‑43 (MQ=255)
                                                                     gTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTa   <  1:1342427/65‑1 (MQ=255)
                                                                     gTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTa   <  1:668286/65‑1 (MQ=255)
                                                                     gTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTa   <  1:2965863/65‑1 (MQ=255)
                                                                     gTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTa   <  1:3416482/65‑1 (MQ=255)
                                                                      |                                                                
CGGCGTTACCTGTACTGTTCCTGACGGCCCGAAGTGAAGAGGTCGATCGCCTGCTTGGGCTGGAAATTGGTGCTGACGACTACGTGGCTAAACCGTTTTCACCCCGCGAAGTGTGCGCCAGGGTGCGCACCTTAC  >  Exported/2982517‑2982651

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: