Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,895,826	A→G	L262L (TTG→CTG)	adeP ←	adenine permease, high affinity; adenine:H+ symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,895,826	0	A	G	97.6%	158.8 / ‑6.4	42	L262L (TTG→CTG)	adeP	adenine permease, high affinity; adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/19);  minor base C (0/1);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAACATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATCCAGCCCAGCAGCG  >  NC_000913/3895680‑3895966                                                                                                                                                   |                                                                                                                                             agTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATg                                                                                                                                          <  1:244926/151‑1 (MQ=255)       gATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAagag                                                                                                                                    <  2:336259/151‑1 (MQ=255)         tAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAgaga                                                                                                                                   >  2:540580/1‑150 (MQ=255)          aaaCGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGaa                                                                                                                                  >  1:12154/1‑150 (MQ=255)              gAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATc                                                                                                                              >  2:482098/1‑150 (MQ=255)                  cGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCAcccc                                                                                                                         <  1:168453/151‑1 (MQ=255)                  aggTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCAcccc                                                                                                                         >  2:205679/2‑151 (MQ=255)                   ggTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCt                                                                                                                        <  1:337636/151‑1 (MQ=255)                       aCGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCa                                                                                                                    <  1:99308/151‑1 (MQ=255)                                   cTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGtt                                                                                                          <  2:593685/149‑1 (MQ=255)                                     gTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTaaa                                                                                                       <  2:151915/150‑1 (MQ=255)                                     gTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAAc                                                                                                      <  2:396844/151‑1 (MQ=255)                                        gACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGAc                                                                                                   >  1:237897/1‑151 (MQ=255)                                        gACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGAc                                                                                                   >  1:633526/1‑151 (MQ=255)                                           aTACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCg                                                                                                >  2:213744/1‑151 (MQ=255)                                           aTACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCATGCCAGCCCGAGGTTAAACGACCCg                                                                                                >  2:93695/1‑151 (MQ=255)                                              cAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCt                                                                                             <  1:31296/151‑1 (MQ=255)                                                         tCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTa                                                                                      >  2:18759/1‑147 (MQ=255)                                                              cgcgGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGAccc                                                                              >  1:596500/1‑150 (MQ=255)                                                                  ggAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAAc                                                                         <  1:678205/151‑1 (MQ=255)                                                                                 cATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTcgg                                                          <  2:353117/151‑1 (MQ=255)                                                                                      gCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGgcggcg                                                     >  1:240252/1‑151 (MQ=255)                                                                                      gCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGgcggcg                                                     >  2:101673/1‑151 (MQ=255)                                                                                       ccAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGgcgc                                                    <  2:422835/151‑1 (MQ=255)                                                                                       ccAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGgcgc                                                    >  1:85876/1‑151 (MQ=255)                                                                                           aCCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGa                                                 <  1:114333/150‑1 (MQ=255)                                                                                              tGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACg                                             >  2:59033/1‑151 (MQ=255)                                                                                                cTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGAt                                           >  1:482015/1‑151 (MQ=255)                                                                                                 tttATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAAGCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAAc                                              <  2:624399/147‑1 (MQ=255)                                                                                                          cACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAg                                 <  1:615665/151‑1 (MQ=255)                                                                                                           aCGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGt                                >  1:461929/1‑151 (MQ=255)                                                                                                                aaTCAGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGcaca                           <  2:605806/151‑1 (MQ=255)                                                                                                                    aGCGTACCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCAc                       >  1:229005/1‑151 (MQ=255)                                                                                                                         aCCGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGc                  >  1:365519/1‑151 (MQ=255)                                                                                                                           cGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCAt                <  2:426998/151‑1 (MQ=255)                                                                                                                           cGGAGGAGTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCAt                <  2:82723/151‑1 (MQ=255)                                                                                                                                gagTCAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATccagc           >  1:206633/1‑151 (MQ=255)                                                                                                                                    cAAACAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATccagccc         >  1:211417/1‑149 (MQ=255)                                                                                                                                      aaCAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATccagcccagc      >  2:570174/1‑150 (MQ=255)                                                                                                                                        cAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATCCAGCCCagcagc   <  2:551963/151‑1 (MQ=255)                                                                                                                                        cAAGTTGACCAGCATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATCCAGCCCagcagc   >  1:316546/1‑151 (MQ=255)                                                                                                                                         aaGCTGACCACCATGAAATAGAAACTCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATCCCGCCGAGCAGcg  <  2:477903/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AGTACCGATAAACGAACCGGTCACGGAAGAGATACTGTCGACATACAGCGCCTGCTTCATGCGCGGGAATTTCCCCTTCTCATCCGCCAGACCTGCTTTATCGGTCACGCCAATCAGCGTACCGGAGGAGTCAAACAAGTTGACCAACATGAAAGAGAAAATCACCCCTGCCAGCCCGAGGTTAAACGACCCGGCTAAATCTACATGACCCACAACTGTCATTACGCTCGGCGGCGCAGAAACGATGCCATTGTAGTGCACATCACCCAGCATCCAGCCCAGCAGCG  >  NC_000913/3895680‑3895966

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