Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,853,909	A→G	intergenic (+267/‑13)	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_000913	3,853,909	0	A	G	97.2%	246.1 / ‑2.2	72	intergenic (+267/‑13)	tisB/emrD	toxic membrane persister formation peptide, LexA‑regulated/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (37/33);  total (38/34)
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.56e-01

TCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAT  >  NC_000913/3853813‑3854004                                                                                                 |                                                                                                tCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCa                                                                                              <  2:691101/100‑1 (MQ=255)       cGTGTATCTGGCTATCATTCATCAATGTGATAGATTATTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTg                                                                                        <  2:960910/100‑1 (MQ=255)           tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                    >  1:781185/1‑100 (MQ=255)           tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                    <  1:842460/100‑1 (MQ=255)           tatCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATaa                                                                                    >  1:866515/1‑100 (MQ=255)            atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                   <  2:1496678/100‑1 (MQ=255)            atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                   <  1:1808193/100‑1 (MQ=255)            atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCACATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                   <  2:796219/100‑1 (MQ=255)                 gCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaaa                                                                              <  2:866512/100‑1 (MQ=255)                  cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                             >  2:1640565/1‑100 (MQ=255)                     aCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCa                                                                          >  2:261364/1‑100 (MQ=255)                       aTTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaaa                                                                        >  1:1757179/1‑100 (MQ=255)                       aTTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaaa                                                                        >  1:1980730/1‑100 (MQ=255)                        ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                       <  1:859495/100‑1 (MQ=255)                          catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                     >  2:1580238/1‑100 (MQ=255)                            tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                   >  2:946977/1‑100 (MQ=255)                            tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                   >  2:683770/1‑100 (MQ=255)                            tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                   <  1:946981/100‑1 (MQ=255)                            tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                   <  1:412679/100‑1 (MQ=255)                            tcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTGTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAAc                                                                   <  2:1254971/100‑1 (MQ=255)                             caATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACg                                                                  >  2:1522110/1‑100 (MQ=255)                              aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGt                                                                 <  2:569337/100‑1 (MQ=255)                               aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTc                                                                <  1:1162642/100‑1 (MQ=255)                               aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTc                                                                >  2:541966/1‑100 (MQ=255)                                tgtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCa                                                               <  2:1850250/100‑1 (MQ=255)                                tgtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCa                                                               >  1:1334737/1‑100 (MQ=255)                                 gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCaa                                                              >  2:1877224/1‑100 (MQ=255)                                 gtgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCaa                                                              >  1:1346800/1‑100 (MQ=255)                                  tgATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAt                                                             >  2:315109/1‑100 (MQ=255)                                      aGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                                                         <  2:175280/100‑1 (MQ=255)                                      aGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGt                                                         <  1:424313/100‑1 (MQ=255)                                         ttCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                      <  2:1857014/100‑1 (MQ=255)                                          tCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                                                     <  1:1263748/100‑1 (MQ=255)                                            ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGa                                                   >  2:1979765/1‑100 (MQ=255)                                               tcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGt                                                <  2:1544340/100‑1 (MQ=255)                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                               >  2:1622609/1‑100 (MQ=255)                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                               >  2:1222874/1‑100 (MQ=255)                                                cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                               <  1:2029559/100‑1 (MQ=255)                                                 ccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTg                                              >  2:131266/1‑100 (MQ=255)                                                 ccGCATTTATGGGAATGCGTAGTGACTTAGTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATg                                                 <  2:1759585/97‑1 (MQ=255)                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                             <  1:1111160/100‑1 (MQ=255)                                                  cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                             <  1:351508/100‑1 (MQ=255)                                                   gCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGt                                            >  2:1819923/1‑100 (MQ=255)                                                    cATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTa                                           <  1:1979783/100‑1 (MQ=255)                                                      tttATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGTTAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTAtt                                         <  2:1134582/100‑1 (MQ=255)                                                        tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                       <  1:1819939/100‑1 (MQ=255)                                                        tATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTAc                                       >  1:295427/1‑100 (MQ=255)                                                          tGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTc                                     >  2:1743165/1‑100 (MQ=255)                                                           gggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCg                                    <  1:1959175/100‑1 (MQ=255)                                                            ggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGt                                   >  2:1322421/1‑100 (MQ=255)                                                              aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                 <  1:459926/100‑1 (MQ=255)                                                                tGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACATCAATTTGTTATTGATGTTGGTATTACTCGTGGcc                               >  2:975392/1‑100 (MQ=255)                                                                  cGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                             >  2:248850/1‑100 (MQ=255)                                                                   gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                            <  1:315110/100‑1 (MQ=255)                                                                   gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                            <  1:541967/100‑1 (MQ=255)                                                                     aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                          >  1:370626/1‑100 (MQ=255)                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                       >  1:477912/1‑100 (MQ=255)                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                       <  1:1845685/100‑1 (MQ=255)                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                       >  2:1561415/1‑100 (MQ=255)                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                       >  2:411414/1‑100 (MQ=255)                                                                        gACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCa                       >  2:390844/1‑100 (MQ=255)                                                                            tATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATg                   <  2:781182/100‑1 (MQ=255)                                                                                cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc               <  1:1750231/100‑1 (MQ=255)                                                                                 tAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa              >  2:388267/1‑100 (MQ=255)                                                                                  aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGTCCGTCGGTCAGATGGCGCa              >  1:533151/1‑99 (MQ=255)                                                                                  aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaa             <  1:1681733/100‑1 (MQ=255)                                                                                    ttattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAc           >  1:1207113/1‑100 (MQ=255)                                                                                     tattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAcc          >  2:1114526/1‑100 (MQ=255)                                                                                            ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtata   >  1:1367569/1‑100 (MQ=255)                                                                                            ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTtata   >  2:1167410/1‑100 (MQ=255)                                                                                             taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat  >  1:526596/1‑100 (MQ=255)                                                                                             taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGGCAATTTGTTATTTATGTTGTTAGTACTCGTGGCCGTCGGTCAGATGGCTCAAACCATTTatat  >  2:678556/1‑100 (MQ=255)                                                                                                 |                                                                                                TCCTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAT  >  NC_000913/3853813‑3854004

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