Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,521,502	A→T	10.9%	intergenic (+174/+91)	lysV → / ← xapR	tRNA‑Lys/transcriptional activator of xapAB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,521,502	0	A	T	10.9%	140.1 / 9.2	46	intergenic (+174/+91)	lysV/xapR	tRNA‑Lys/transcriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base T (3/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

TTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTC  >  NC_000913/2521356‑2521649                                                                                                                                                   |                                                                                                                                                    ttttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaa                                                                                                                                                    <  2:670275/148‑1 (MQ=32)  tttACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                   <  2:698655/148‑1 (MQ=32)    tACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAgg                                                                                                                                                 >  1:722113/1‑148 (MQ=31)                    cAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTg                                                                                                                                >  4:293691/1‑149 (MQ=32)                     aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               >  4:839615/1‑149 (MQ=32)                     aaaGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGc                                                                                                                               >  4:53265/1‑149 (MQ=32)                           tCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAttt                                                                                                                          >  2:636819/1‑148 (MQ=37)                                    cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGc                                                                                                                <  3:293691/149‑1 (MQ=255)                                               agttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCtt                                                                                                     <  3:179934/149‑1 (MQ=255)                                                 ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTc                                                                                                    <  2:722113/148‑1 (MQ=255)                                                                  acGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCAGATGCGCTTCACTTAGCAGACTACTAtt                                                                                  >  3:218413/1‑149 (MQ=255)                                                                      aGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTcc                                                                               >  1:122271/1‑148 (MQ=255)                                                                            tCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAAtt                                                                        <  3:53265/149‑1 (MQ=255)                                                                                 gcgACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCAGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGt                                                                   <  4:218413/149‑1 (MQ=255)                                                                                       gCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcct                                                              >  1:332738/1‑148 (MQ=255)                                                                                       gCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcct                                                              <  2:100211/148‑1 (MQ=255)                                                                                        cccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTctcctc                                                             >  1:218334/1‑148 (MQ=255)                                                                                            aaGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCt                                                         >  2:272613/1‑148 (MQ=255)                                                                                                   agcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTgtgt                                                  >  2:556640/1‑148 (MQ=255)                                                                                                          cgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCa                                           >  1:66590/1‑148 (MQ=255)                                                                                                                  ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTg                                                                    >  2:44979/1‑115 (MQ=255)                                                                                                                  ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTg                                                                    <  1:44979/115‑1 (MQ=255)                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGc                                                                                                                        >  3:282684/1‑62 (MQ=37)                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGc                                                                                                                        <  4:282684/62‑1 (MQ=37)                                                                                                                   cTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGc                                  <  1:636819/148‑1 (MQ=255)                                                                                                                     tctcGACGCACCAATGTAAAAAAGCGCACTAAAGGCGCTTTGTTGCTATCTGGGATTTGCGCAATTGCCGGACGCGCTTCACTTTGTAGACTACTATTTCCTGCAATTCCTGTCTCCTCACGTACTGTGTCTATGCAGCCAACTGCtt                                <  2:656893/145‑1 (MQ=255)                                                                                                                      cACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTaa                              <  3:839615/149‑1 (MQ=255)                                                                                                                       aCGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTaa                              <  2:66590/148‑1 (MQ=255)                                                                                                                         gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACtgt                                                    <  4:106523/124‑1 (MQ=255)                                                                                                                         gACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACtgt                                                    >  3:106523/1‑124 (MQ=255)                                                                                                                          aCCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCAGTTAGCAGACTAGTGTTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAATAGCTTaaaga                           >  1:656893/1‑145 (MQ=38)                                                                                                                                caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACa                                    <  3:508468/133‑1 (MQ=255)                                                                                                                                caATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACa                                    >  4:508468/1‑133 (MQ=255)                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACtgt                                                    >  3:845156/1‑111 (MQ=255)                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACtgt                                                    <  4:845156/111‑1 (MQ=255)                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCa               <  2:218334/148‑1 (MQ=255)                                                                                                                                       aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCAcc             >  3:568994/1‑149 (MQ=255)                                                                                                                                        aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCAcc             >  1:297179/1‑148 (MQ=255)                                                                                                                                          aaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACtgtg                                                   <  2:664659/108‑1 (MQ=255)                                                                                                                                          aaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACtgtg                                                   >  1:664659/1‑108 (MQ=255)                                                                                                                                             cgcCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACctgctg        <  1:675195/148‑1 (MQ=255)                                                                                                                                               cccTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCtgtg      <  2:332738/148‑1 (MQ=255)                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCa               <  2:25262/137‑1 (MQ=255)                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCa               >  1:25262/1‑137 (MQ=255)                                                                                                                                                 cTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGtt    <  2:120919/148‑1 (MQ=255)                                                                                                                                                  tAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTGCTATTTCCGGCAATTCCTGTCTCCTCACCTACTGCGTCAATGCAGCCAACAGATTAACCATCGCGGGCGTCACCTGCTGTGTTTc  >  3:762893/1‑149 (MQ=255)                                                                                                                                                  tAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGttt   <  2:122271/148‑1 (MQ=255)                                                                                                                                                  tAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGttt   >  1:852540/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TTTTACTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATTTGCGAAATTGCCTGATGCGCTTCACTTAGCAGACTACTATTTCCGGCAATTCCTGTCTCCTCACCTACTGTGTCAATGCAGCCAACAGCTTAACCATCGCGGGCGTCACCTGCTGTGTTTC  >  NC_000913/2521356‑2521649

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