Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,519,252	G→A	intergenic (+47/+5)	yfeD → / ← gltX	putative DNA‑binding transcriptional regulator YfeD/glutamate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,519,252	0	G	A	98.0%	171.8 / NA	51	intergenic (+47/+5)	yfeD/gltX	putative DNA‑binding transcriptional regulator YfeD/glutamate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/19);  minor base T (0/1);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 3.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTACGGCGACACG  >  NC_000913/2519115‑2519388                                                                                                                                          |                                                                                                                                         agacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACt                                                                                                                                        >  2:156926/1‑140 (MQ=255) agacaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACt                                                                                                                                        >  1:223988/1‑140 (MQ=255)   acaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTaa                                                                                                                                      >  1:1014424/1‑140 (MQ=255)   acaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTaa                                                                                                                                      >  2:177848/1‑140 (MQ=255)   acaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTaa                                                                                                                                      >  2:26092/1‑140 (MQ=255)   acaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTaa                                                                                                                                      >  1:901348/1‑140 (MQ=255)   acaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTaa                                                                                                                                      >  2:788729/1‑140 (MQ=255)   acaTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTaa                                                                                                                                      >  1:453343/1‑140 (MQ=255)     aTTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAAtt                                                                                                                                    >  1:371149/1‑140 (MQ=255)              gAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGAt                                                                                                                           >  2:326367/1‑140 (MQ=255)                   aGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCAAAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTcg                                                                                                                      >  2:294551/1‑140 (MQ=255)                       gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGtt                                                                                                                  <  2:958935/140‑1 (MQ=255)                       gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCAAAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGtt                                                                                                                  <  1:294551/140‑1 (MQ=255)                        cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTc                                                                                                                 >  2:555086/1‑140 (MQ=255)                            aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCa                                                                                                             <  2:901348/140‑1 (MQ=255)                                    cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                                                                                     <  1:650857/140‑1 (MQ=255)                                       gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                                                                                  >  2:859671/1‑140 (MQ=255)                                         gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                                                                                                <  1:156926/140‑1 (MQ=255)                                           aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                                                                                              <  2:825930/140‑1 (MQ=255)                                           aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                                                                                              <  2:453343/140‑1 (MQ=255)                                                 tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGA‑TTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                                                                                       >  2:800650/1‑140 (MQ=255)                                                  cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                                                                                       >  1:119125/1‑140 (MQ=255)                                                   gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                                                                                      <  2:333136/140‑1 (MQ=255)                                                           cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                                                                              <  2:987141/140‑1 (MQ=255)                                                               cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACt                                                                          >  1:1010733/1‑140 (MQ=255)                                                                   tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                                      <  1:26092/140‑1 (MQ=255)                                                                   tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                                      >  1:459994/1‑140 (MQ=255)                                                                   tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg                                                                      <  1:788729/140‑1 (MQ=255)                                                                     tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGt                                                                    <  1:971195/140‑1 (MQ=255)                                                                          atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc                                                               <  1:326367/140‑1 (MQ=255)                                                                              aCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAAt                                                           >  2:268783/1‑140 (MQ=255)                                                                                 aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGc                                                        >  2:825560/1‑140 (MQ=255)                                                                                    cAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTg                                                     <  2:75452/140‑1 (MQ=255)                                                                                      gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACGCTCTAATTACTGCTGATTTTCGCGTGCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGaa                                                   <  2:477892/140‑1 (MQ=255)                                                                                      gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGaa                                                   <  2:1014424/140‑1 (MQ=255)                                                                                         taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGA‑TTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGt                                               <  1:800650/140‑1 (MQ=255)                                                                                            ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTaa                                             >  1:878573/1‑140 (MQ=255)                                                                                                 aaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATcc                                        >  1:254344/1‑140 (MQ=255)                                                                                                     ggCAGGATATTATCTCCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTg                                    <  1:825560/140‑1 (MQ=255)                                                                                                                 tctcCTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGcc                        <  2:254344/140‑1 (MQ=255)                                                                                                                    ccTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCg                     >  2:41472/1‑140 (MQ=255)                                                                                                                    ccTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCg                     >  2:532081/1‑140 (MQ=255)                                                                                                                    ccTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCg                     >  2:875988/1‑140 (MQ=255)                                                                                                                    ccTGCCGTTTATCTTTTTACAAACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCg                     >  2:674916/1‑140 (MQ=255)                                                                                                                     cTGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGc                    >  2:175959/1‑140 (MQ=255)                                                                                                                      tGCCGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCa                   <  2:1010733/140‑1 (MQ=255)                                                                                                                        ccGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCTCCGCAcc                 >  2:149573/1‑140 (MQ=255)                                                                                                                        ccGTTTATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAcc                 >  1:599594/1‑140 (MQ=255)                                                                                                                             tATCTTTTTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTa            >  2:1021963/1‑140 (MQ=255)                                                                                                                                     tACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTACGGCGaca    >  2:248066/1‑140 (MQ=255)                                                                                                                                     tACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTACGGCGaca    >  2:391657/1‑140 (MQ=255)                                                                                                                                       cacaCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTACGGCGACACg  <  2:241270/140‑1 (MQ=255)                                                                                                                                       cacaCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTACGGCGACACg  <  1:555950/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         AGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCACCGGTTACGGCGACACG  >  NC_000913/2519115‑2519388

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