Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,853,677 A→C 15.7% intergenic (+35/‑245) tisB → / → emrD toxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,853,6770AC15.7% 65.5 / 6.1 25intergenic (+35/‑245)tisB/emrDtoxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (9/12);  new base C (2/2);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

GCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTA  >  NC_000913/3853529‑3853822
                                                                                                                                                    |                                                                                                                                                 
gCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGTGGAGCGGGaaa                                                                                                                                                   <  3:918589/149‑1 (MQ=255)
gCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                   <  3:845828/149‑1 (MQ=255)
gCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGaaa                                                                                                                                                   <  2:32419/149‑1 (MQ=255)
    cTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGc                                                                                                                                               <  3:188117/149‑1 (MQ=255)
      gAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTc                                                                                                                                             <  1:623669/149‑1 (MQ=255)
      gAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTc                                                                                                                                             <  2:39611/149‑1 (MQ=255)
         ccACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTcccc                                                                                                                                          >  4:835304/1‑149 (MQ=255)
           acaGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCtt                                                                                                                                        <  2:387112/149‑1 (MQ=255)
            caGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGTGGGAAACCGCTCCCCttt                                                                                                                                       >  2:1102629/1‑149 (MQ=255)
                    gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGTGGGAAACCGCTCCCCTTTTATATTTa                                                                                                                               <  1:1102629/149‑1 (MQ=255)
                    gcgTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTa                                                                                                                               >  4:834413/1‑149 (MQ=255)
                              cTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCtttg                                                                                                                                      <  4:725755/132‑2 (MQ=255)
                              cTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCtttg                                                                                                                                      >  3:725755/1‑131 (MQ=255)
                               tggtggATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGTTTCCCGCTCCCCtt                                                                                                                                        <  2:355093/129‑1 (MQ=255)
                                             aTTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATTCCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcga                                                                                                      <  3:834413/149‑1 (MQ=255)
                                             aTTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcga                                                                                                      <  3:835304/149‑1 (MQ=255)
                                                            cTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTa                                                                                       >  2:332377/1‑149 (MQ=255)
                                                                       agcaCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTAcc                                                                            <  1:332377/149‑1 (MQ=255)
                                                                          aCTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCAtc                                                                         >  1:957800/1‑149 (MQ=255)
                                                                                  tGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcga                                                                                                      >  3:793446/1‑112 (MQ=255)
                                                                                  tGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTAGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGgatgcga                                                                                                      <  4:793446/112‑1 (MQ=255)
                                                                                         tGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACAtt                                                          <  2:957800/149‑1 (MQ=255)
                                                                                                                  ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAACCCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTAt                                                       >  4:469194/1‑127 (MQ=255)
                                                                                                                  ttCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTAt                                                       <  3:469194/127‑1 (MQ=255)
                                                                                                                                      aaGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCt                                                >  3:255148/1‑114 (MQ=255)
                                                                                                                                      aaGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCt                                                <  4:255148/114‑1 (MQ=255)
                                                                                                                                           gAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCt        >  1:817875/1‑149 (MQ=255)
                                                                                                                                                 gAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGta  >  3:1074008/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTCCTGAACCACAGGAGACGCGTATGAACCTGGTGGATATCGCCATTCTTATCCTCAAACTCATTGTTGCAGCACTGCAACTGCTTGATGCTGTTCTGAAATACCTGAAGTAATTCAGATTCAAGTCGCACCAAAGGGGAGCGGGAAACCGCTCCCCTTTTATATTTAGCGTGCGGGTTGGTGTCGGATGCGATGCTGACGCATCTTATCCGCCCTACCATCTCTCCCGGCAACATTTATTGCCGCTTTTGTTTACATATTCTGCCGCTAAACAATTCCCCATTCCTGGCGTA  >  NC_000913/3853529‑3853822

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: