Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	1,385,093	T→G	17.7%	intergenic (+71/‑44)	aspS → / → pmpR	Aspartate‑‑tRNA(Asp/Asn) ligase/Transcriptional regulatory protein PmpR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	1,385,093	0	T	G	17.7%	80.9 / 3.6	46	intergenic (+71/‑44)	aspS/pmpR	Aspartate‑‑tRNA(Asp/Asn) ligase/Transcriptional regulatory protein PmpR
Reads supporting (aligned to +/- strand):  ref base T (10/26);  new base G (7/2);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 1.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGT  >  ACB143/1384966‑1385226                                                                                                                                 |                                                                                                                                      gCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCCTCCCCAAGGATGCCCCGATGCGGGGGGACGCCCGGTTTTTACCCGCCCCCCCCCG‑CACggggc                                                                                                                               >  2:210699/1‑136 (MQ=255)   cAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCCTCCACATTGTTGGGCCGTTTGGGTTAGTCG‑ACGGGTTTTTACCGTGCCGCCCCGGGATggggcgg                                                                                                                             >  8:361618/1‑136 (MQ=255)        cccTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCCTCCACATTGCTGCGCCATTTAGTTTGGACG‑CCGGTTTTTTACCCTCCCCCCCCCTCAGGGGGGGGGGGCg                                                                                                                        >  7:371236/1‑136 (MQ=255)        cccTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCCCATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCg                                                                                                                        >  5:192654/1‑136 (MQ=255)         ccTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCCTCCACATGGTTGCGCCTTTGGGTTTCGGCG‑CAGGTTTTTTCCCGTCCCGCCCCGGCATTGGGGGGGGGCGt                                                                                                                       >  3:295314/1‑136 (MQ=255)              gcgAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCCTACCCAAGGATGCGCCTTTGGGGTTGGACG‑CCGGGGTTTTACCGGCCCGCCCCCGCCTggggcgggggcgggtggt                                                                                                                  >  8:217000/1‑130 (MQ=255)              gcgAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCCTACACATGGATGGGCCTTTTGGTTTAGGCG‑CAGGGTTTTTACCGTACCCCCCCGTAAATGGGCGCGGCCGgtttgg                                                                                                                  >  5:355990/1‑133 (MQ=255)              gcgAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGCTGCGCCTTTGCGTTTCGACG‑CCGGTTTTTTCCCGTCCCCCCCCGGCATGGGGGCGGGCCGGTtgtg                                                                                                                  >  1:97088/1‑136 (MQ=255)              gcgAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCCCCTCAACATGGATGCGCCTTTGCGGTTGGGCG‑CAGGTTTTTTCCCGTCCCGCCCCGGAATGGGGGGGGGCCGtgtggg                                                                                                                  >  7:113655/1‑134 (MQ=255)                gAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGAt                                                                                                                <  6:149955/136‑1 (MQ=255)                  aCTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTc                                                                                                              <  5:391210/136‑1 (MQ=255)                    tGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGGCATGGGGCGGGGGCGGTTGTGGttga                                                                                                             >  5:384970/1‑133 (MQ=255)                     gCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCTCGGGCGCCTCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCaa                                                                                                            >  2:233251/1‑135 (MQ=255)                           ccGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTc                                                                                                     <  1:272971/136‑1 (MQ=255)                            cGACTGCGCGTACAGACCAAGGGCGAGTAAGCAGGCCCCGGGCGGATCCACATGGATGCGCATTTTCGTTTCGACG‑CAGGTTTTTTCCCGTGCGGCCCCGTAATGGGGCTGGGCCGTTTGTTATTCAAAGGATTc                                                                                                     <  8:310995/135‑1 (MQ=255)                            cGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTc                                                                                                     <  1:84564/135‑1 (MQ=255)                            cGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTc                                                                                                     <  3:225153/135‑1 (MQ=255)                                              aaGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCGTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGGCCTGGGGGGGGGGCGTGTGTGGTTGAAAGGATTGGGGGTCGTGTTGGGTgggc                                                                                  >  6:161612/1‑134 (MQ=255)                                                 ggcgtgTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCAtt                                                                               <  7:331582/131‑1 (MQ=255)                                                          gCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTggg                                                                      <  2:396112/136‑1 (MQ=255)                                                              tCCCCGGGCGCATCGACTTGGATGCTCCTTTGCGTTTCGACG‑CTGGTTTTTTCCCGTCCCGCCCCGTGATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGa                                                                   <  6:150753/135‑1 (MQ=255)                                                               ccccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTGCAAGTGGGCgca                                                                  >  6:381324/1‑133 (MQ=255)                                                                 ccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACAt                                                               <  8:46788/136‑1 (MQ=255)                                                                  cGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATc                                                              <  5:126569/136‑1 (MQ=255)                                                                   gggCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCa                                                             <  5:362304/136‑1 (MQ=255)                                                                    ggCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCaa                                                            <  6:130199/136‑1 (MQ=255)                                                                          tCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCg                                                      <  7:66180/136‑1 (MQ=255)                                                                            cacaTGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGc                                                     <  7:203072/135‑1 (MQ=255)                                                                                  gATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTGAAAGGATTAGGAGTCGGTTTGGGGGGTGATTACAAGTGGGCGGACATGAAGCACCCCAAAGAgc                                              >  8:203072/1‑136 (MQ=255)                                                                                    tGCTCCTTTGCGTATCGACG‑CAGTTTTTTTCCCGTCCCGCGCCGTCATGGTGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAgcgc                                            <  7:86538/136‑1 (MQ=255)                                                                                    tGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAgcgc                                            <  4:187919/136‑1 (MQ=255)                                                                                      cgcCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTAAAAGGATTAGGAGTCGTTTTGGCTGGGCATTACAAAGTGGCGGACCTCAAGCCCCCCAAAGAGCGccc                                          >  2:2415/1‑135 (MQ=255)                                                                                         cTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGGTTGAAAGGGTTAGGATGCGTTTTGGCGTGTGATTCGAAGTTGGCGAACATACAGGACCCCCAAGAGGGCCCggg                                       >  1:199539/1‑135 (MQ=255)                                                                                             gCGTTTCCACG‑CTGTTTTTTTCCCGCCCCGCCCGGTTTTGGGGGGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGcc                                   <  4:15960/136‑1 (MQ=255)                                                                                                   cGACG‑CTGGTTTTTTCCGGTCCCGCCCCGGCTTGGGGGGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                             <  1:210699/136‑1 (MQ=255)                                                                                                    gACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagagag                            <  2:352723/136‑1 (MQ=255)                                                                                                    gACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagagag                            <  3:192408/136‑1 (MQ=255)                                                                                                    gACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagagag                            <  3:313396/136‑1 (MQ=255)                                                                                                    gACG‑CAGGTTTTTTCCCGTCCCGCCCCGGCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                             <  3:396650/135‑1 (MQ=255)                                                                                                    gACG‑AAGTTTTTTTCCCGTCCCGCCCCGTCAGGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAagaga                             <  2:407177/135‑1 (MQ=255)                                                                                                      cG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTGAAAGGATTCGGAGTCGTTTTGGCTGGTCATTGCCAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAgg                           >  1:24736/1‑135 (MQ=255)                                                                                                      cG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTGGTTTTTGGGGGGGATTGCAAGTGGGCGAAAAAGAAACACCGCAAAGAGCGCCCGGAAGCCAAAAGAgg                           >  4:66749/1‑135 (MQ=255)                                                                                                          aGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAg                       <  2:56605/136‑1 (MQ=255)                                                                                                          aGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAg                       <  1:331504/136‑1 (MQ=255)                                                                                                                         ccgCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTgg        <  8:170982/136‑1 (MQ=255)                                                                                                                               cGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGt  <  3:236683/136‑1 (MQ=255)                                                                                                                                 |                                                                                                                                      GCCAAGGCCCTGCGCGAACTGCACATCCGACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACG‑CAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCCAAGAGAGGCAAGGTCTTCACCAAGTGGATCCGT  >  ACB143/1384966‑1385226

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