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	95.5%	211.2 / 0.2	67	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/2);  new base G (35/28);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853815‑3854005                                                                                               |                                                                                                 cTGTCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGAGTACTCTCTATTTCTAATTATTTTTATAAAGCCATc                                                                                             >  2:1373495/1‑100 (MQ=255) cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             >  2:891374/1‑100 (MQ=255) cTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATc                                                                                             >  2:1483338/1‑100 (MQ=255)     cGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTg                                                                                         <  2:560086/100‑1 (MQ=255)      gTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGa                                                                                        <  1:660824/100‑1 (MQ=255)       tataTCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGAt                                                                                       <  1:1020540/100‑1 (MQ=255)          atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    <  2:873629/100‑1 (MQ=255)          atCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAAt                                                                                    >  1:1156623/1‑100 (MQ=255)           tCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   <  2:634041/100‑1 (MQ=255)             tGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATg                                                                                   >  2:1429754/1‑98 (MQ=255)              ggCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGaaa                                                                                <  1:717572/100‑1 (MQ=255)                cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  1:395702/1‑100 (MQ=255)                cTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAg                                                                              >  2:140450/1‑100 (MQ=255)                 tAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAgg                                                                             <  2:574546/100‑1 (MQ=255)                  aaCATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATAAAAAGGc                                                                            <  1:377791/100‑1 (MQ=255)                    cATTCATCAATGTGATAGATTCCTCTTCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCaa                                                                          >  1:934442/1‑100 (MQ=255)                      ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        <  2:787635/100‑1 (MQ=255)                      ttCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAg                                                                        >  2:539433/1‑100 (MQ=255)                       tcatcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGa                                                                       <  1:414077/100‑1 (MQ=255)                        catcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaa                                                                      >  2:447241/1‑100 (MQ=255)                         atcaATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGaaa                                                                     >  2:1410896/1‑100 (MQ=255)                            aaTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGt                                                                  <  2:489063/100‑1 (MQ=255)                             aTGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTc                                                                 <  1:778540/100‑1 (MQ=255)                                 gATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAtt                                                             <  2:809791/100‑1 (MQ=255)                                   tAGATTCCTCTCCCGCAATTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAAttt                                                            >  1:803973/1‑99 (MQ=255)                                     gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         >  1:444158/1‑100 (MQ=255)                                     gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         >  2:1190982/1‑100 (MQ=255)                                     gATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGtt                                                         <  1:140450/100‑1 (MQ=255)                                       ttCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAt                                                       <  2:1243740/100‑1 (MQ=255)                                        tCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTAtt                                                      <  1:956553/100‑1 (MQ=255)                                          ctctCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGCTAATGAAAAGGGAAAGAAACGTCAATTTGTTATTg                                                     >  2:690885/1‑99 (MQ=255)                                            ctcCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATg                                                  >  2:493074/1‑100 (MQ=255)                                              cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGTCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:501272/1‑100 (MQ=255)                                              cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  1:1285565/1‑100 (MQ=255)                                              cccGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGtt                                                >  2:1141539/1‑100 (MQ=255)                                                cGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTgg                                              <  2:1031591/100‑1 (MQ=255)                                                   aTTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTAt                                           <  2:1324689/100‑1 (MQ=255)                                                       aTGGGAACGCGCAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAGGAAACGTCAATTTGTTATTGATGTTGGTATTACt                                       >  2:829382/1‑100 (MQ=255)                                                         gggAATGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCg                                     <  2:1258258/100‑1 (MQ=255)                                                            aaTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTgg                                  <  1:234483/100‑1 (MQ=255)                                                             aTGCGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGc                                 <  2:1164554/100‑1 (MQ=255)                                                                cGTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGt                              >  1:569368/1‑100 (MQ=255)                                                                 gTAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTc                             <  2:853416/100‑1 (MQ=255)                                                                  tAGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCg                            <  2:1195437/100‑1 (MQ=255)                                                                   aGTGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCgg                           <  1:1540672/100‑1 (MQ=255)                                                                     tGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTc                         <  2:1247250/100‑1 (MQ=255)                                                                     tGACTTATTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTc                         >  2:771119/1‑100 (MQ=255)                                                                           aTTCTAATTATTTTTATAAATGCATCAGGGATACTGAAAATGAGAGGAAAAAACAAATGGAGATTGATGGTAGTATGACTCGTGGCCCTCGGGCataag                    >  2:1577080/1‑95 (MQ=255)                                                                           aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTTGCCGTCGGTCAGATg                    >  1:798514/1‑99 (MQ=255)                                                                           aTTCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATgg                   >  2:785049/1‑100 (MQ=255)                                                                            ttCTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgc                  >  2:236672/1‑100 (MQ=255)                                                                              cTAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGgcgc                >  2:413199/1‑100 (MQ=255)                                                                               tAATTATTTTTATAAAGGCATCCGTGATAATGAAAAGACAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa               >  2:1348604/1‑100 (MQ=255)                                                                                aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGACACGTCAATCTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGAAGGCGCaa              >  2:1526190/1‑100 (MQ=255)                                                                                aaTTATTTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCa               >  2:1450430/1‑99 (MQ=255)                                                                                 ttccgatctataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaaa             <  1:433111/92‑1 (MQ=255)                                                                                 attattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGGCAATGTGTTATTGATGGTGGTATGACTCTTGGCCGTCGGTCAGATGGCGCaaa             >  2:684554/1‑100 (MQ=255)                                                                                   tattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAcc           <  1:1141544/100‑1 (MQ=255)                                                                                   tattTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAAcc           <  1:785054/100‑1 (MQ=255)                                                                                    attTTTATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCa          >  2:457865/1‑100 (MQ=255)                                                                                     tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaa              >  2:678490/1‑95 (MQ=255)                                                                                     tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAt         <  1:478927/100‑1 (MQ=255)                                                                                     tttttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCAt         >  1:1182240/1‑100 (MQ=255)                                                                                        ttATAAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTta      <  1:978516/100‑1 (MQ=255)                                                                                          ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTTGGTCAGATGGCGCAAACCATTtata    >  1:1415343/1‑100 (MQ=255)                                                                                          ataAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCaaaagatcggaa    >  2:433110/1‑91 (MQ=255)                                                                                           taAAGGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTatat   >  2:493510/1‑100 (MQ=255)                                                                                            aaagGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATAtt  >  2:1143938/1‑100 (MQ=255)                                                                                               |                                                                                                 CTGGCGTATATCTGGCTAACATTCATCAATGTGATAGATTCCTCTCCCGCATTTATGGGAATGCGTAGTGACTTATTCTAATTATTTTTATAAAAGCATCCGTGATAATGAAAAGGCAAAGAAACGTCAATTTGTTATTGATGTTGGTATTACTCGTGGCCGTCGGTCAGATGGCGCAAACCATTTATATT  >  NC_000913/3853815‑3854005

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