Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,479,651 C→A 30.0% intergenic (+110/+197) 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,6510CA30.0% 77.3 / 33.6 60intergenic (+110/+197)phoU/yieEDNA‑binding transcriptional regulator/predicted phosphopantetheinyl transferase
Reads supporting (aligned to +/- strand):  ref base C (31/11);  new base A (12/6);  total (43/17)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATATATAACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGTT  >  minE/2479584‑2479713
                                                                   |                                                              
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:225340/68‑1 (MQ=255)
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:1751905/68‑1 (MQ=255)
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:2309644/68‑1 (MQ=255)
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:1697484/68‑1 (MQ=255)
   tataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                >  1:3304777/1‑65 (MQ=255)
     taACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAg                                                           >  1:1096288/1‑68 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAAGTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         <  1:611400/68‑1 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:1028778/1‑68 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:2035650/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:253650/68‑1 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         >  1:798157/1‑68 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                         <  1:1130176/68‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:2860036/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:1292757/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:2927669/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:215078/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:853601/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:226505/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:1681085/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:579602/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:794789/58‑1 (MQ=255)
           tataAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCAtt                                                     >  1:1185433/1‑68 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                >  1:1135310/1‑53 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                >  1:1084871/1‑53 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:3397268/68‑1 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:3312561/68‑1 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:1163806/68‑1 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                 <  1:738028/68‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:2458256/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:2931164/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:1817045/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:1707929/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:13601/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                               <  1:1102037/67‑1 (MQ=255)
                             cacaAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                                <  1:2771019/39‑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:177130/1‑69 (MQ=255)
                                          aCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAAttt                     >  1:1844729/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:1606774/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:2825702/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:2146175/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:3373437/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:1981137/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:2939127/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:1417776/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                  >  1:3055620/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCt                                                  >  1:521086/1‑35 (MQ=25)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2398709/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:3428907/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:3325816/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:1876941/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:779129/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:1139486/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:80761/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:1103878/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2191389/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2239381/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa                >  1:2320033/1‑69 (MQ=255)
                                                   aTGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTg            <  1:3030932/69‑1 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:697953/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:1082227/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:3145438/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:300823/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt         >  1:2882058/1‑69 (MQ=255)
                                                        cATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATg        <  1:1308042/68‑1 (MQ=255)
                                                            tttttATCAATAGCGCATTGCTATTTTCTCTGCGCGCAATTAAATTAATTTCCGAACCTGGATGTTCGt   >  1:694315/1‑69 (MQ=255)
                                                              tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:685253/1‑68 (MQ=255)
                                                              tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:1647799/1‑68 (MQ=255)
                                                              tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:2042137/1‑68 (MQ=255)
                                                              tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGtt  >  1:2431020/1‑68 (MQ=255)
                                                              tttATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGt   >  1:1813416/1‑67 (MQ=255)
                                                                   |                                                              
ATATATAACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGTT  >  minE/2479584‑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: