Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,246,455 T→C 36.6% I56I (ATT→ATC thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,246,4550TC36.6% 39.3 / 43.5 71I56I (ATT→ATCthiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (14/31);  new base C (8/18);  total (22/49)
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.84e-01

TATATCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTC  >  minE/2246387‑2246517
                                                                    |                                                              
tataTCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAAtt                                                                >  1:1091433/1‑69 (MQ=255)
tataTCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAAtt                                                                >  1:3345218/1‑69 (MQ=255)
             cacacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg                                                   <  1:451683/69‑1 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:1747805/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:3096274/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:3344853/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:2897262/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:2594491/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:776149/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:2365402/1‑69 (MQ=255)
             cacaaCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATCGGCGGTAGCAAAg                                                   >  1:104328/1‑69 (MQ=255)
              acacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg                                                   <  1:1835611/68‑1 (MQ=255)
              acacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg                                                   <  1:1694545/68‑1 (MQ=255)
               cacCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAg                                                   <  1:184800/67‑1 (MQ=255)
                      gcgtgcgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGc                                          <  1:2246164/69‑1 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  >  1:2225698/1‑56 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  >  1:2038279/1‑56 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  <  1:2359927/56‑1 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  <  1:3135629/56‑1 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  >  1:3316878/1‑56 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  >  1:1252622/1‑56 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  <  1:1165627/56‑1 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa                                      >  1:3117512/1‑68 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa                                      >  1:3020654/1‑68 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa                                      >  1:1321960/1‑68 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa                                      >  1:2777417/1‑68 (MQ=255)
                           cgcgtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTa                                      >  1:1056653/1‑68 (MQ=255)
                              gtgcCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGa                                                  <  1:433743/53‑1 (MQ=255)
                                           agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:315064/69‑1 (MQ=255)
                                           agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:3171746/69‑1 (MQ=255)
                                           agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:2861871/69‑1 (MQ=255)
                                           agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:496848/69‑1 (MQ=255)
                                           agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:2596257/69‑1 (MQ=255)
                                           agaTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:1348332/69‑1 (MQ=255)
                                             aTCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGa                     <  1:2827732/67‑1 (MQ=255)
                                               ccAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTcc                 >  1:575560/1‑69 (MQ=255)
                                               ccAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTc                  >  1:2874530/1‑68 (MQ=255)
                                                     tAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCt            >  1:1804548/1‑68 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACTAAGAAAACGAAGCGATTCCGGTCTAc          <  1:1117983/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:1716686/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:3235438/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:1235182/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:1965124/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:2824673/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:2115288/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:665987/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:2651602/67‑1 (MQ=255)
                                                        cccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:2384376/67‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:660826/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:503640/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:445378/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:685838/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:151068/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:2276408/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:1747325/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:1004102/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:270696/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:1753570/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:2773711/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:2800052/45‑1 (MQ=255)
                                                         ccGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:1950466/45‑1 (MQ=255)
                                                         ccGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:141613/66‑1 (MQ=255)
                                                          cGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:139209/44‑1 (MQ=255)
                                                          cGACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:134344/65‑1 (MQ=255)
                                                           gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:784606/64‑1 (MQ=255)
                                                           gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:566420/64‑1 (MQ=255)
                                                           gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:1944916/64‑1 (MQ=255)
                                                           gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:136276/64‑1 (MQ=255)
                                                           gACGCTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:2328650/64‑1 (MQ=255)
                                                              gCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACgaaga                                >  1:3318277/1‑39 (MQ=255)
                                                              gCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTc  >  1:940574/1‑69 (MQ=255)
                                                               cTAATCGGCGGTAGCAAAGAACAACCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTAc          <  1:2889076/60‑1 (MQ=255)
                                                                  aTTGGCGGTAGCAAAGAACAGCCGCAGTACgaagaa                               <  1:3222601/36‑1 (MQ=255)
                                                                    |                                                              
TATATCACTGGCACACACCCCGGCGTGCGCGTGCCGATGCGTGAGATCCAGCTTAGCCCGACGCTAATTGGCGGTAGCAAAGAACAGCCGCAGTACGAAGAAAACGAAGCGATTCCGGTCTACGACACCTC  >  minE/2246387‑2246517

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: