Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,962,363 T→C 43.8% E331E (GAA→GAG yjjI ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,962,3630TC43.8% 19.5 / 70.4 89E331E (GAA→GAGyjjIconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/25);  new base C (17/22);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAA  >  minE/2962298‑2962419
                                                                 |                                                        
cTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTca                                                       <  1:920902/69‑1 (MQ=255)
cTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTca                                                       <  1:2706501/69‑1 (MQ=255)
cTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTca                                                       <  1:2745754/69‑1 (MQ=255)
cTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTca                                                       <  1:1736066/69‑1 (MQ=255)
cTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTca                                                       <  1:3385725/69‑1 (MQ=255)
                cTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGTCAATCCCTTCTTTTTCACACAGCAAGTTAACCg                                       >  1:2150663/1‑69 (MQ=255)
                cTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGTCAATCCCTTCTTTTTCACACAGCAAGTTAACCg                                       >  1:1928964/1‑69 (MQ=255)
                cTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGTCAATCCCTTCTTTTTCACACAGCAAGTTAACCg                                       >  1:2217331/1‑69 (MQ=255)
                cTTCATTTGCGGAGGCTTCTTTACCGTAGCGCGCGTCAATCCCTTCTTTTTCACACAGCAAGTTAACCg                                       >  1:1112162/1‑69 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCa                                                >  1:943954/1‑58 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCa                                                >  1:439163/1‑58 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCa                                                >  1:1331713/1‑58 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCa                                                >  1:1578084/1‑58 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCt                                     >  1:349944/1‑69 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCt                                     >  1:1175460/1‑69 (MQ=255)
                  tCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCt                                     >  1:2907656/1‑69 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1244951/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:2193816/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1050773/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1826656/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:761160/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:2881236/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1750474/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1134063/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:3090890/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1655020/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:1347089/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:3407096/69‑1 (MQ=255)
                                ttcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       <  1:473218/69‑1 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATa                       >  1:985480/1‑68 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:3411346/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:2952226/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:2975286/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:680551/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:1023757/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:2553523/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:1701407/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:2428590/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:1873888/1‑69 (MQ=255)
                                 tcttTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      >  1:1229664/1‑69 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCtt                                    >  1:2502243/1‑53 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:3191440/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:3246491/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:3312579/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:3088551/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:3373674/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:351779/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:721331/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:1107963/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:2831965/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:1975759/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:965371/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:2765755/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:1970816/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:2198091/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:1826129/1‑54 (MQ=255)
                                   tttACCGTAGCGCGCGGCAATCCCTTCTTTCTCACACAGCAAGTTAACCGCTTc                                   >  1:184968/1‑54 (MQ=255)
                                    ttACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACa                     <  1:2121863/67‑1 (MQ=255)
                                            gcgcgcGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACAt           >  1:745388/1‑69 (MQ=255)
                                            gcgcgcGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACAt           >  1:3371212/1‑69 (MQ=255)
                                            gcgcgcGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACAt           >  1:3079938/1‑69 (MQ=255)
                                                    cAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACaa   <  1:1100160/69‑1 (MQ=255)
                                                     aaTCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACaaa  <  1:2328375/69‑1 (MQ=255)
                                                     aaTCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACaaa  <  1:2456396/69‑1 (MQ=255)
                                                     aaTCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACaaa  <  1:664171/69‑1 (MQ=255)
                                                     aaTCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACaaa  <  1:1373248/69‑1 (MQ=255)
                                                     aaTCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACaaa  <  1:1303700/69‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:484705/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:624075/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:1259340/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:860604/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:1013744/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:2977215/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:2098609/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:2120582/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:1815294/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:2763097/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:2829787/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:1761006/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:2886297/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:1736592/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:348212/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:3018785/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:3056841/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:3206452/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:3316527/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:3367277/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:3404506/42‑1 (MQ=255)
                                                            tcttTCTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATAc                      <  1:1307396/42‑1 (MQ=255)
                                                                 |                                                        
CTGATGCGATAACCTACTTCATTTGCGGCGGCTTCTTTACCGTAGCGCGCGGCAATCCCTTCTTTTTCACACAGCAAGTTAACCGCTTCCGCCAGCCCATACATGCCAAACATTGGCACAAA  >  minE/2962298‑2962419

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: