Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,090,126	C→T	9.7%	intergenic (+80/‑66)	hisL → / → hisG	his operon leader peptide/ATP phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,090,126	0	C	T	9.7%	94.1 / 2.3	42	intergenic (+80/‑66)	hisL/hisG	his operon leader peptide/ATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/22);  new base T (2/2);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

ACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGC  >  NC_000913/2089984‑2090268                                                                                                                                               |                                                                                                                                               acaGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCCGGCGATGTGTGCTGGGAGACATTCAGCTCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGGAGATCAc                                                                                                                                                >  2:46292/1‑143 (MQ=255)   agagACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGCTGCATGAACGCATGAGTCAGCCTCCGGCAGCTCAccg                                                                                                                                              >  1:839179/1‑142 (MQ=255)         ttttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCggg                                                                                                                                        <  2:497502/143‑1 (MQ=255)          tttATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGGTGATCTTCCgggg                                                                                                                                       <  2:42952/143‑1 (MQ=255)           ttATGACACGCGTTCTATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCggggg                                                                                                                                      <  2:139332/143‑1 (MQ=255)                acacGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCtttt                                                                                                                                 >  1:545194/1‑143 (MQ=255)                     cgTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTAtt                                                                                                                            <  2:367921/143‑1 (MQ=255)                       ttCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGACCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgc                                                                                                                          <  1:579628/143‑1 (MQ=255)                        tCAATCTAAACACCCCCATCATCACCATCAGCCTGACTAGTCTGTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTgcg                                                                                                                         <  2:764075/143‑1 (MQ=255)                           aTTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCgg                                                                                                                      <  1:574749/143‑1 (MQ=255)                            tttAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGACGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAACCCCCGAAAGATCACCTTCCGGGGGCTTTATTATTAcgcgat                                                                                                                     >  1:744253/1‑141 (MQ=255)                                   accatcATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAAc                                                                                                              <  1:53302/143‑1 (MQ=255)                                              caccatcatcCTGACTAGTCTTTCAGGCGATGTGTGCTTGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAg                                                                                                   <  2:418510/143‑1 (MQ=255)                                                 catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAGCATATAAATAATCGTTTGGGGGCAACCCAACCGGGGGCTTTTTTATTGCGCGGTCGATAACGGTTCAGACAGGtt                                                                                                >  2:188241/1‑143 (MQ=255)                                                 catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGCCAggat                                                                                                >  2:19051/1‑141 (MQ=255)                                                 catcatcCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGtt                                                                                                >  1:322828/1‑143 (MQ=255)                                                   gcagccTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATGCAGAGCTGCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTa                                                                                              <  1:556057/139‑1 (MQ=255)                                                       ccTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAaga                                                                                          <  2:147422/143‑1 (MQ=255)                                                         tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGGAAGCCCCCGGAAGATCACCTTCCGGGGGCTGTTTTATTGCGCGGTTGATGACGGTTCTGTGAGGttttcagggg                                                                                        >  2:667109/1‑136 (MQ=255)                                                         tGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAgg                                                                                        >  1:429801/1‑143 (MQ=255)                                                              aGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATaa                                                                                   <  2:465090/143‑1 (MQ=255)                                                               gTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATGGCGCGGGTGATCACGGTTCAGACGGGTTTGAGGAGGAATagc                                                                                  >  1:779094/1‑141 (MQ=255)                                                                tCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACa                                                                                 <  1:147570/143‑1 (MQ=255)                                                                  tttCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTATAAAGAGGAaaaaaaaa                                                                               >  1:213687/1‑143 (MQ=255)                                                                  ttgcAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaa                                                                               <  2:657435/140‑1 (MQ=255)                                                                   ttCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACaaaa                                                                              <  1:771167/143‑1 (MQ=255)                                                                      aGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATga                                                                           >  1:18294/1‑143 (MQ=255)                                                                                 tcTGGAAGACATTCAGATCTTCCCGTGGTGCAACAGCCTGCGTTTTGTCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCTGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACt                                                                <  2:867460/142‑1 (MQ=255)                                                                                     gAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGAGTCTCTTGTATTGCGCGGTTGATAACGGTTCAGACAGGTTCAACGAGCAATAACAAAATGACAGACACAACTCgct                                                            >  2:47451/1‑141 (MQ=255)                                                                                         aCATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCCTTTTCATTGCGCGGTTCATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACg                                                        <  1:227842/143‑1 (MQ=255)                                                                                            ttCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAt                                                     >  2:566434/1‑143 (MQ=255)                                                                                             tCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATa                                                    <  1:770581/143‑1 (MQ=255)                                                                                                    ttCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGccacac                                             >  1:196533/1‑138 (MQ=255)                                                                                                            gtgCTTGAATGCATGAGAAAGCACCCGGAAGATCAACTTCGGGGGGCGTTGTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAGTGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATcc                                     <  1:565634/143‑1 (MQ=255)                                                                                                             tgCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCg                                    <  2:70611/143‑1 (MQ=255)                                                                                                              gcttGAACGCATGAGAAAGCGCCCGGAAGATCACCTTCCGGGGGCTCTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAACGAGGAATAACAAAATGACAGTCAACCCTCGTTTACGCATAGCTACGCAGAAATCCgg                                   <  2:616844/140‑1 (MQ=255)                                                                                                                                 cccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCa                <  1:240806/143‑1 (MQ=255)                                                                                                                                  ccccGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCAGAGCTATGCAGAAATCCGGCCGATTAAGTGATGACTCAc               >  1:470384/1‑143 (MQ=255)                                                                                                                                   cccTGAAGATCACCTCCCGGGGGCGTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAAGAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCAcg              <  1:46292/143‑1 (MQ=255)                                                                                                                                   cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAACGAGGAAGAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGCCCGCGTAAGTGATGACTCAcg              >  2:2019/1‑143 (MQ=255)                                                                                                                                   cccGGAAGGTGATCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGACTAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATACCGCCGTTTAAGTGATGACTCAca              >  1:872791/1‑142 (MQ=255)                                                                                                                                       gAAGATCACCTTCCGGGGGCGTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGaa          <  1:681250/143‑1 (MQ=255)                                                                                                                                            tCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGCTTAAGTGATGACTCACGCGAATTGCt     >  2:354508/1‑143 (MQ=255)                                                                                                                                             taCCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGGTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACTCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTg    <  2:340493/142‑1 (MQ=255)                                                                                                                                               ccTTCCGGGGGCTTTTTTATTGCGCGGTCGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGgc  <  1:188241/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                               ACAGAGACTTTTATGACACGCGTTCAATTTAAACACCACCATCATCACCATCATCCTGACTAGTCTTTCAGGCGATGTGTGCTGGAAGACATTCAGATCTTCCAGTGGTGCATGAACGCATGAGAAAGCCCCCGGAAGATCACCTTCCGGGGGCTTTTTTATTGCGCGGTTGATAACGGTTCAGACAGGTTTAAAGAGGAATAACAAAATGACAGACAACACTCGTTTACGCATAGCTATGCAGAAATCCGGCCGTTTAAGTGATGACTCACGCGAATTGCTGGC  >  NC_000913/2089984‑2090268

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