Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,483,912 T→A 11.0% V53D (GTC→GAC)  evgA → DNA‑binding transcriptional activator EvgA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,483,9120TA11.0% 86.5 / 3.2 37V53D (GTC→GAC) evgADNA‑binding transcriptional activator EvgA
Reads supporting (aligned to +/- strand):  ref base T (18/12);  new base A (2/2);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

TTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGG  >  NC_000913/2483771‑2484047
                                                                                                                                             |                                                                                                                                       
ttGATGACCATCCTCTTGCTATCGGAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                                                        <  1:75112/143‑1 (MQ=255)
ttGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTc                                                                                                                                        <  1:666197/143‑1 (MQ=255)
     gaCCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGatat                                                                                                                                   >  1:413654/1‑143 (MQ=255)
         ttCCTCTTGCTAGCGCAGTAATTCGTAATTTATTGATCACAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATcccc                                                                                                                               <  2:712881/142‑1 (MQ=255)
         gtCCTGTTGCTATCGCAGCAATTCGTAATTTATTGATCAAGAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGACATCAATGATGACGATATcagg                                                                                                                               <  2:499937/142‑4 (MQ=255)
                   tATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                     >  2:511561/1‑143 (MQ=255)
                   ggtCGCAGCCATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACg                                                                                                                     <  1:69862/141‑1 (MQ=255)
                     tCGCAGCAATTCGTAATTTATTGATTAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGt                                                                                                                   <  2:218164/143‑1 (MQ=255)
                         agcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGAGATCCACGGCGTTAACGGTATcc                                                                                                               >  2:225824/1‑143 (MQ=255)
                          gcaATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCa                                                                                                              <  1:181336/143‑1 (MQ=255)
                               tcgggaTTTAGTGAGCAAAGAGGATAGTGGAATCTGAGCAGAGGTGAGTGAAGGCGGAAGGGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGAGATCCCCGGAGTTAACGGTATCCAGgtgt                                                                                                         <  2:673415/138‑1 (MQ=255)
                                  tAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAg                                                                                                      >  1:812279/1‑143 (MQ=255)
                                   aagtgaTTGATCAACAACGATAGTGAGATCGTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGa                                                                                                     <  1:12237/138‑1 (MQ=255)
                                             tCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGGATCGAGGTGTTAGAAACGCTGAgg                                                                                           >  1:86145/1‑143 (MQ=255)
                                              cAAAAACGATATTCAAATCTTAGCAGAGTTGACTGAAGGCTGAAGTGCCGTTCAGCGCGTGGTGACACTTAAGCCTGCTATCGTCAGCATTGATAGCGATATCACCGGAGGTCACGGTGTCCACGTGTTAGAAAAGCTGCGGa                                                                                          >  2:196819/1‑143 (MQ=255)
                                                       tatTGAAATCTTAGCAGAGTTGACTGAAGGCGGACGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAat                                                                                 >  1:244809/1‑143 (MQ=255)
                                                          tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCactcta                                                                              >  2:456957/1‑138 (MQ=255)
                                                          tGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAAtata                                                                              >  1:748983/1‑143 (MQ=255)
                                                           gAAATCTTGGCAGGGTGGACTGAAGGCGGAAGGGCCGGTCGGCGGGTGGAAACACTTAAGCGTGATATCGTCATCATTGATGGCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAg                                                                             <  2:280819/143‑1 (MQ=255)
                                                           gAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGCTAGAAACGCTGAGGAAGCGCCCATATAg                                                                             >  2:868555/1‑143 (MQ=255)
                                                                          ttGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATtattattat                                                              >  2:786937/1‑143 (MQ=255)
                                                                               tggaGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAACGGAATTGTTATTATCGtac                                                         >  2:285940/4‑141 (MQ=255)
                                                                               tGAAGGCGGAAGTGCCGTTCAGCGGGAGGAAACACTTAAGCCTGATAGCGGCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGCGGAAGAGCCAATATAGCGGAATAATTAGTATCGtct                                                         >  1:860777/1‑143 (MQ=255)
                                                                                ggaGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGtctc                                                        <  1:470806/141‑1 (MQ=255)
                                                                                          gTGCCGTGCAGCGGGTGGAAACCCTTAAGCCTTATGTCCGCATCATTGATGTCGATATCCTCGGGGTTAACTGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAAt                                              <  2:54723/143‑1 (MQ=255)
                                                                                           tGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTAGATATACCCGGAGTTAACGGTATCAAGGTGTTAGAAACGCTGAGGAAGCCCGAATATAGAGAAATTATGATTATCGTCTCCGCTcaaaggt                                             >  2:342890/1‑140 (MQ=255)
                                                                                            gcagTTCAGCGGGTGGAAGCACGTAAGCCGGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGa                                            <  1:729269/140‑1 (MQ=255)
                                                                                                    gCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCAttttt                                    >  1:556527/1‑143 (MQ=255)
                                                                                                        gTGGAAACACTTAAGCCTGATATCGTCATGATTGATGTCGATATCCCAGGAGTTAACGGTATCCAGTTGTTCGACACGCTGAGGCCGCGCCACTAAACCCGAACCATACTTCTCCTCTCACCTCCCAACGAACATTTTTACgg                                >  2:557583/1‑143 (MQ=255)
                                                                                                         gggAGACCCTTCAGCCTGATATCGACATCATTGAGGGCGATCACCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCGGAGGCAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACggg                               <  1:264202/142‑1 (MQ=255)
                                                                                                            aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATCTCCCCGGAGTTAACGGTATCCCGGTGTTAGACACTCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCATCTCCGCTAAAAATGACCATTTTTACGGGaaa                            >  1:329965/1‑143 (MQ=255)
                                                                                                            aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACTCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGaaa                            <  1:532295/143‑1 (MQ=255)
                                                                                                            aaaCACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCACGGAGTTAACGGTATCCAGGTGTTAGAAGCGCTGAGGAAGCGCCAATATAGCCGAATTATTATTGTCGTCTCCGCTAAAAAAGACCATTTTTACGGCaaa                            >  2:752111/1‑143 (MQ=255)
                                                                                                                            tatCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAGACATTGTGCTGATGCTg            >  1:516911/1‑143 (MQ=255)
                                                                                                                             atCGACATCAATGATGACGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTgg           >  1:432147/1‑143 (MQ=255)
                                                                                                                               cGTCATCATTGATGACGATATCCCCGGAGTTAACGGTATAAGGGTGTTATAGTCACTGAGGAAGCGCCAATATATCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGgcg         <  1:834934/143‑1 (MQ=255)
                                                                                                                                      atTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAAAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATAGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATgg  >  1:311938/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                       
TTGATGACCATCCTCTTGCTATCGCAGCAATTCGTAATTTATTGATCAAAAACGATATTGAAATCTTAGCAGAGTTGACTGAAGGCGGAAGTGCCGTTCAGCGGGTGGAAACACTTAAGCCTGATATCGTCATCATTGATGTCGATATCCCCGGAGTTAACGGTATCCAGGTGTTAGAAACGCTGAGGAAGCGCCAATATAGCGGAATTATTATTATCGTCTCCGCTAAAAATGACCATTTTTACGGGAAACATTGTGCTGATGCTGGCGCTAATGG  >  NC_000913/2483771‑2484047

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