Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,479,657 C→A 36.7% intergenic (+116/+191) phoU → / ← yieE DNA‑binding transcriptional regulator/predicted phosphopantetheinyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,479,6570CA36.7% 43.4 / 33.8 49intergenic (+116/+191)phoU/yieEDNA‑binding transcriptional regulator/predicted phosphopantetheinyl transferase
Reads supporting (aligned to +/- strand):  ref base C (23/8);  new base A (12/6);  total (35/14)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGTT  >  minE/2479591‑2479713
                                                                  |                                                        
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAAGTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         <  1:611400/68‑1 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         <  1:253650/68‑1 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:2035650/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:1028778/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:1873511/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:1631754/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:2601044/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:1240783/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:798157/1‑68 (MQ=255)
aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         <  1:1130176/68‑1 (MQ=255)
    tataAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCAtt                                                     >  1:1185433/1‑68 (MQ=255)
        aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:1163806/68‑1 (MQ=255)
        aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:738028/68‑1 (MQ=255)
        aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:3312561/68‑1 (MQ=255)
        aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:3397268/68‑1 (MQ=255)
          atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:13601/67‑1 (MQ=255)
          atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:1707929/67‑1 (MQ=255)
          atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:2458256/67‑1 (MQ=255)
          atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:1817045/67‑1 (MQ=255)
          atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:1102037/67‑1 (MQ=255)
          atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:2931164/67‑1 (MQ=255)
                               aataaCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCaattaaatta                         <  1:2655714/69‑1 (MQ=255)
                               aataaCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCaattaaatta                         <  1:3165677/69‑1 (MQ=255)
                                   aCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAAttt                     >  1:980876/1‑69 (MQ=255)
                                   aCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAAttt                     >  1:1844729/1‑69 (MQ=255)
                                   aCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAAttt                     >  1:177130/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:2146175/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:1981137/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:3373437/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:1606774/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:1417776/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:2825702/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:2939127/1‑69 (MQ=255)
                                      gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:3055620/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCt                                                  >  1:521086/1‑35 (MQ=25)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:3428907/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:3325816/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:779129/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:80761/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2398709/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:1103878/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:1139486/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2320033/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2239381/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2191389/1‑69 (MQ=255)
                                        gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:1876941/1‑69 (MQ=255)
                                            aTGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTg            <  1:3030932/69‑1 (MQ=255)
                                               gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:1082227/1‑69 (MQ=255)
                                               gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:697953/1‑69 (MQ=255)
                                               gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:2882058/1‑69 (MQ=255)
                                               gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:3145438/1‑69 (MQ=255)
                                               gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:300823/1‑69 (MQ=255)
                                                 cATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATg        <  1:1308042/68‑1 (MQ=255)
                                                     tttttATCAATAGCGCATTGCTATTTTCTCTGCGCGCAATTAAATTAATTTCCGAACCTGGATGTTCGt   >  1:694315/1‑69 (MQ=255)
                                                       tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:2042137/1‑68 (MQ=255)
                                                       tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:685253/1‑68 (MQ=255)
                                                       tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:2431020/1‑68 (MQ=255)
                                                       tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:1647799/1‑68 (MQ=255)
                                                       tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGt   >  1:1813416/1‑67 (MQ=255)
                                                                  |                                                        
ACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGTT  >  minE/2479591‑2479713

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: