Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,479,644 T→C 25.2% intergenic (+103/+204) 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,6440TC25.2% 81.2 / 24.7 71intergenic (+103/+204)phoU/yieEDNA‑binding transcriptional regulator/predicted phosphopantetheinyl transferase
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base C (12/6);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATATAACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGT  >  minE/2479584‑2479712
                                                            |                                                                    
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:1751905/68‑1 (MQ=255)
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:225340/68‑1 (MQ=255)
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:1697484/68‑1 (MQ=255)
atatataACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:2309644/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:2601044/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:1028778/1‑68 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                        >  1:798157/1‑68 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                        >  1:1240783/1‑68 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                        <  1:253650/68‑1 (MQ=255)
       aCTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAgcg                                                        <  1:1130176/68‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:226505/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:1681085/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:853601/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:1292757/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:215078/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:2927669/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:794789/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:579602/58‑1 (MQ=255)
          ttATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:2860036/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:738028/68‑1 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                <  1:3312561/68‑1 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                <  1:3397268/68‑1 (MQ=255)
               aaTTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTa                                                <  1:1163806/68‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                              <  1:1817045/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                              <  1:2931164/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                              <  1:1707929/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                              <  1:13601/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                              <  1:2458256/67‑1 (MQ=255)
                 atCCTAAAATTACACAAAATTAATAACTGCGAGCATGGTCATACTTTTATAAATAGAGCATTGCTAtt                                              <  1:1102037/67‑1 (MQ=255)
                             cacaAAGTTAATAACTGCGAGCATGGTCATATTTTTATc                                                               <  1:2771019/39‑1 (MQ=255)
                                      aataaCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCaattaaatta                        <  1:3165677/69‑1 (MQ=255)
                                      aataaCTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCaattaaatta                        <  1:2655714/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:2825702/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:2146175/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:1981137/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:2939127/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:3055620/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:1606774/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:1417776/1‑69 (MQ=255)
                                             gCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCg                 >  1:3373437/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCt                                                 >  1:521086/1‑35 (MQ=25)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:779129/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:3428907/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:80761/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:2320033/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:3325816/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:2398709/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)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:1139486/1‑69 (MQ=255)
                                               gAGCATGGTCATACTTTTATAAATAGAGCATTGCTATTTTTTCTGCACGCAACTAAATTAATTTTCGaa               >  1:1103878/1‑69 (MQ=255)
                                                   aTGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTg           <  1:3030932/69‑1 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt        >  1:3145438/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt        >  1:300823/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt        >  1:2882058/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt        >  1:697953/1‑69 (MQ=255)
                                                      gTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGAt        >  1:1082227/1‑69 (MQ=255)
                                                        cATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATg       <  1:1308042/68‑1 (MQ=255)
                                                            tttttATCAATAGCGCATTGCTATTTTCTCTGCGCGCAATTAAATTAATTTCCGAACCTGGATGTTCGt  >  1:694315/1‑69 (MQ=255)
                                                            |                                                                    
ATATATAACTTTATAAATTCCTAAAATTACACAAAGTTAATAACTGCGAGCATGGTCATATTTTTATCAATAGCGCATTGCTATTTTCTCTGCACGCAATTAAATTAATTTCCGAACCTGGATGTTCGT  >  minE/2479584‑2479712

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: