Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,860,448 C→T 8.9% G295E (GGG→GAG)  gllD ← 4‑oxalomesaconate tautomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,860,4480CT 8.9% 123.8 / 2.3 45G295E (GGG→GAG) gllD4‑oxalomesaconate tautomerase
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CAAACGGGCCGTCCTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACATTGCGCTGG  >  NC_002947/2860303‑2860594
                                                                                                                                                 |                                                                                                                                                   
cAACCGGGCCGTCCTCACCAACCCGCACCCGTTGATCACCCCATGTTCCAGGCTGATTTCCACCGTAATTTCCCCAGTCGGATGCGCCACCGCCAAACGCTGACGATCACCCCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCtt                                                                                                                                                  <  5:96010/149‑1 (MQ=255)
   aCGGGCCGTCCTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTtca                                                                                                                                                >  6:55043/1‑148 (MQ=255)
   aCGGGCCGTCCTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTtca                                                                                                                                                <  5:55043/148‑1 (MQ=255)
       gCCGTCCTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTtcaatcaa                                                                                                                                           <  4:254982/149‑1 (MQ=255)
           tCCTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACCTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAg                                                                                                                                       >  5:305544/1‑149 (MQ=255)
             cTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGc                                                                                                                                     <  8:332114/149‑1 (MQ=255)
              tCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGcc                                                                                                                                    >  3:96890/1‑149 (MQ=255)
                  cAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAg                                                                                                                                <  3:201815/149‑1 (MQ=255)
                         cACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGt                                                                                                                                  >  2:212002/1‑140 (MQ=255)
                         cACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCAGTGCAACCGACCCTTCAATCAAGCAGGCCGt                                                                                                                                  <  1:212002/140‑1 (MQ=255)
                            cccTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCAc                                                                                                                      >  1:244205/1‑149 (MQ=255)
                              cTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTtcaatcaa                                                                                                                                           >  3:324818/1‑126 (MQ=255)
                              cTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTtcaatcaa                                                                                                                                           <  4:324818/126‑1 (MQ=255)
                                     aCGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCaaa                                                                                                             <  6:305544/149‑1 (MQ=255)
                                       gCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAaca                                                                                                           >  4:262102/1‑149 (MQ=255)
                                                       aTTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCAt                                                                                           <  4:96890/148‑1 (MQ=255)
                                                       aTTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATg                                                                                          >  3:328993/1‑149 (MQ=255)
                                                           ccACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACAccc                                                                                                        <  5:136754/132‑1 (MQ=255)
                                                           ccACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACAccc                                                                                                        >  6:136754/1‑132 (MQ=255)
                                                            cACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGc                                                                                     >  3:175759/1‑149 (MQ=255)
                                                                        cccACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATg                                                                          >  5:27138/1‑148 (MQ=255)
                                                                                gATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAgcgc                                                                 <  2:244205/149‑1 (MQ=255)
                                                                                    cTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTa                                                              <  6:27138/148‑1 (MQ=255)
                                                                                           gCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGt                                                      >  2:96632/1‑149 (MQ=255)
                                                                                            ccAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATg                                                                          >  6:99526/1‑128 (MQ=255)
                                                                                            ccAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATg                                                                          <  5:99526/128‑1 (MQ=255)
                                                                                                  gCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACTCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAgcg                                                                  <  7:221181/130‑1 (MQ=255)
                                                                                                  gCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACTCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAgcg                                                                  <  7:219200/130‑1 (MQ=255)
                                                                                                  gCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACTCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAgcg                                                                  >  8:219200/1‑130 (MQ=255)
                                                                                                  gCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACTCTTCAATCAAGCAGGCCGTAGCAACACACACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAgcg                                                                  >  8:221181/1‑130 (MQ=255)
                                                                                                        gATCACCGCCGGACGTGCTTGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGGAACACTCACCGCACCAAACACCCCGATCGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTgcg                                          <  1:96632/149‑1 (MQ=255)
                                                                                                            aCCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGt                                      <  4:175759/148‑1 (MQ=255)
                                                                                                             ccgccgGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGa                                                          >  4:265850/1‑127 (MQ=255)
                                                                                                             ccgccgGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGa                                                          <  3:265850/127‑1 (MQ=255)
                                                                                                               gccgGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGAGGGATGACCGAGCGCGTATTGACCGGGCCGCCGTTGCGCGGTGCCg                                  >  8:219222/1‑149 (MQ=255)
                                                                                                                       tGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAgg                          <  3:262102/149‑1 (MQ=255)
                                                                                                                       tGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAgg                          >  3:342269/1‑149 (MQ=255)
                                                                                                                              aaGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCt                   >  5:44162/1‑149 (MQ=255)
                                                                                                                                        gCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACAtt         >  6:92609/1‑149 (MQ=255)
                                                                                                                                           aCCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACATTgcg      >  6:127923/1‑149 (MQ=255)
                                                                                                                                            ccGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACAtt         >  4:63488/1‑145 (MQ=255)
                                                                                                                                            ccGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACAtt         <  3:63488/145‑1 (MQ=255)
                                                                                                                                               aCCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACATTGCgctgg  >  7:236893/1‑149 (MQ=255)
                                                                                                                                               aCCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACATTGCgctgg  >  7:29495/1‑149 (MQ=255)
                                                                                                                                               aCCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGCACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACATTGCgctgg  <  7:219222/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAAACGGGCCGTCCTCACCAACCCGCACCCCTTGATCACGCCATGTTCCAGGCTGATTTCCACCGTAAATTCCCCACTCGGATGCTCCACCGCCAAACGCTGACGATCACCGCCGGACGTGCTGGCAAGGCCCTGTGCAACCGACCCTTCAATCAAGCAGGCCGTAGCAACACTCACCGCACCAAACACCCCGATGGACGCATGGCAGCGGTGCGGGATGAACGAGCGCGTATTGACCGTGCCGCCGTTGCGCGGTGCCGACAGCAGGCACATCTTCGGCACATTGCGCTGG  >  NC_002947/2860303‑2860594

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