Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,835,651 G→A 22.7% N82N (AAC→AAT yjeJ ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,835,6510GA22.7% 125.9 / 27.0 92N82N (AAC→AATyjeJhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/40);  new base A (8/13);  total (39/53)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

TCCAGCGTGTTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCATACTGAAGGTG  >  minE/2835586‑2835704
                                                                 |                                                     
tCCAGCGTGTTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGttt                                                    >  1:2672414/1‑69 (MQ=255)
       tgtTAAGCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCa                                             >  1:43160/1‑69 (MQ=255)
       tgtTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCa                                             >  1:2421653/1‑69 (MQ=255)
       tgtTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCa                                             >  1:357516/1‑69 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           >  1:2850976/1‑69 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           >  1:2403235/1‑69 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           <  1:1656586/69‑1 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           >  1:144532/1‑69 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           <  1:298252/69‑1 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAt                                            >  1:385157/1‑68 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAt                                            >  1:2972529/1‑68 (MQ=255)
         ttAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCAt                                            >  1:832826/1‑68 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:2321629/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:142251/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:1483120/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:2341772/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:3365575/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:1208510/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:2703790/69‑1 (MQ=255)
              ccGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCtt                                      <  1:200606/69‑1 (MQ=255)
                    aTTAATATCGGCATTTTTAAGCTCATCAACCAGACTTTCTGGGATGTTTTCTGCCATTTTCTTAATGAc                                <  1:3234700/69‑1 (MQ=255)
                    aTTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGAc                                <  1:2307319/69‑1 (MQ=255)
                      tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:1789476/1‑69 (MQ=255)
                      tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:228767/1‑69 (MQ=255)
                      tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:2371510/1‑69 (MQ=255)
                      tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:484164/1‑69 (MQ=255)
                      tAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACtt                              >  1:2890901/1‑69 (MQ=255)
                           tCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATc                          >  1:2021458/1‑68 (MQ=255)
                           tCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATc                          >  1:285813/1‑68 (MQ=255)
                           tCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATc                          >  1:1545291/1‑68 (MQ=255)
                               cATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAg                     <  1:508822/69‑1 (MQ=255)
                                aTTTTTAAGCTCATGAACCAGAATTTCTGGGATGTTTTCTGCCAtt                                           <  1:3051006/46‑1 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTg                                                >  1:1569559/1‑41 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:2365184/1‑69 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:3387066/1‑69 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:1901480/1‑69 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:341710/1‑69 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:1366767/1‑69 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:53663/1‑69 (MQ=255)
                                aTTTTTAAGCTCATCAACCAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGc                    >  1:1182404/1‑69 (MQ=255)
                                 tttttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCt                   <  1:2748227/69‑1 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:3104024/58‑1 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:245318/58‑1 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:592017/58‑1 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:2938744/58‑1 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATcgcg                       >  1:218738/1‑63 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:336702/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1022593/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:462575/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:597492/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:869678/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:2154491/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1004129/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1127895/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:2728221/1‑69 (MQ=255)
                                   tttAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTg                 >  1:1331043/1‑69 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:2990791/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:3348979/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:49112/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:2268788/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:1355350/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:432817/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:1556750/68‑1 (MQ=255)
                                      aaGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:156291/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCTTTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:1416975/66‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:30890/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:1965484/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:210621/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2856811/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2359659/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2624923/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:2861124/68‑1 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCAt            >  1:2807949/1‑69 (MQ=255)
                                        gCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCAt            >  1:2811911/1‑69 (MQ=255)
                                         cTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:3374591/67‑1 (MQ=255)
                                                aCCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTa                            <  1:2597152/45‑1 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:963050/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:1683346/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:2928939/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:733946/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:72005/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:699341/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:1231824/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:1271975/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:3008960/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:3135102/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:1897617/69‑2 (MQ=255)
                                                 ccAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:405248/69‑2 (MQ=255)
                                                  cAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCATACTGAAGgtg  <  1:1382226/69‑1 (MQ=255)
                                                  cAGAATTTCTGGGATATTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGTTCATACTGAAgg  <  1:2901062/68‑2 (MQ=255)
                                                      aTTTCTGGGATGTTTTCTGCAATTTTCTTAATGACTTTATCGCGAGCTTGCt               <  1:3328470/52‑1 (MQ=255)
                                                         tCTTGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCa             <  1:1424159/51‑1 (MQ=255)
                                                                 |                                                     
TCCAGCGTGTTAACCCGACGATTAATATCGGCATTTTTAAGCTCATCAACCAGAATTTCTGGGATGTTTTCTGCCATTTTCTTAATGACTTTATCGCGAGCTTGCTCATACTGAAGGTG  >  minE/2835586‑2835704

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