Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 561,102 T→G 11.7% G37G (GGA→GGC ybjI ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE561,1020TG11.7% 144.4 / 11.0 68G37G (GGA→GGCybjIconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base G (5/3);  total (40/28)
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

AGAAAGAGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGTTTT  >  minE/561050‑561158
                                                    |                                                        
agaaagaGATCAACTGATAATATTGATTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCAttt                                          <  1:1419605/69‑1 (MQ=255)
agaaagaGATCAACTGATAATATTGATTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCAttt                                          <  1:535100/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       >  1:1180669/1‑69 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1222210/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:987380/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1278303/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:3082425/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:467667/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1599280/69‑1 (MQ=255)
   aagaGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1593330/69‑1 (MQ=255)
     gagaTCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGa                                     <  1:2411488/69‑1 (MQ=255)
         tCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCt                                                 >  1:839606/1‑53 (MQ=255)
          cAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTg                                <  1:2968938/69‑1 (MQ=255)
          cAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTg                                <  1:1899833/69‑1 (MQ=255)
               gATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCa                           <  1:724936/69‑1 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCAt                          >  1:1038079/1‑68 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:3156038/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:1777281/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:3089602/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:608046/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:1187064/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:1899866/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:2658106/1‑69 (MQ=255)
                 taataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:2615524/1‑69 (MQ=255)
                  aataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATaa                        >  1:333750/1‑69 (MQ=255)
                  aataTTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATaa                        >  1:3042470/1‑69 (MQ=255)
                   atatTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTg                                >  1:991144/1‑60 (MQ=255)
                   atatTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCa                           <  1:2301322/65‑1 (MQ=255)
                   atatTGATTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCAttt                                          <  1:1247393/50‑1 (MQ=39)
                      ttGATTCCCGCTGGCGACTACAAAGCGTATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCg                    >  1:2953193/1‑69 (MQ=255)
                      ttGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:2128406/50‑1 (MQ=255)
                      ttGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:2282175/50‑1 (MQ=255)
                      ttGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:715017/50‑1 (MQ=255)
                      ttGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1252068/50‑1 (MQ=255)
                      ttGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCg                    >  1:2393640/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:3028892/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:487390/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:1456582/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:2040785/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:2486461/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:1227762/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:2182620/1‑69 (MQ=255)
                        gATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCt                  >  1:1498156/1‑69 (MQ=255)
                         aTTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:192953/1‑61 (MQ=255)
                         aTTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:1883768/1‑61 (MQ=255)
                         aTTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:3219571/1‑61 (MQ=255)
                         aTTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:1063884/1‑61 (MQ=255)
                         aTTCCCGCTGGCGACCACAAAGCGAATGCCTTGTGCTTTCATTTGCTGATACTGAGCCATa                         >  1:222002/1‑61 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:927504/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:992311/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1551029/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:1879587/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:233382/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:2266229/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCt                                       <  1:2131965/45‑1 (MQ=255)
                           tCCCGCTGGCGACTACAAAGCGAATTACTTGTGCTTTCATTTGCt                                       <  1:3129043/45‑1 (MQ=255)
                            cccGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAg                              >  1:1685895/1‑53 (MQ=255)
                            cccGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTc                 >  1:375997/1‑66 (MQ=255)
                            cccGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACg              >  1:3140662/1‑69 (MQ=255)
                            cccGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACg              >  1:1941220/1‑69 (MQ=255)
                            cccGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACg              >  1:295652/1‑69 (MQ=255)
                                 tGGCGACTACAAAGCTAATTCCTTGTGCTTTCATTTGCt                                       <  1:258302/39‑1 (MQ=255)
                                        tACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGtttt  >  1:3417014/1‑69 (MQ=255)
                                        tACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGtttt  >  1:1285573/1‑69 (MQ=255)
                                        tACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGtttt  >  1:1962954/1‑69 (MQ=255)
                                        tACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGtttt  >  1:2105275/1‑69 (MQ=255)
                                        tACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGtttt  >  1:881731/1‑69 (MQ=255)
                                        tACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGtttt  >  1:1013532/1‑69 (MQ=255)
                                                    |                                                        
AGAAAGAGATCAACTGATAATATTGATTCCCGCTGGCGACTACAAAGCGAATTCCTTGTGCTTTCATTTGCTGATACTGAGCCATAAACCGCTCACGGTTATAGGTTTT  >  minE/561050‑561158

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: