Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,519,228	C→A	13.0%	intergenic (+23/+29)	yfeD → / ← gltX	DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,519,228	0	C	A	13.0%	150.0 / 13.1	46	intergenic (+23/+29)	yfeD/gltX	DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (21/19);  new base A (3/3);  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 = 2.31e-01

AATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAC  >  NC_000913/2519088‑2519372                                                                                                                                             |                                                                                                                                                 aaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                                                                                 >  1:26794/1‑142 (MQ=255) aaTCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                                                                                 <  2:26794/142‑1 (MQ=255)                      tatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                                             >  3:261857/1‑125 (MQ=255)                      tatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                                             >  3:524884/1‑125 (MQ=255)                      tatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                                             <  4:524884/125‑1 (MQ=255)                      tatGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTg                                                                                                                                             <  4:261857/125‑1 (MQ=255)                                 ttGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                                                                                  <  4:989140/109‑1 (MQ=255)                                 ttGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                                                                                  >  3:989140/1‑109 (MQ=255)                                 ttGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                                                                             <  3:113881/146‑1 (MQ=255)                                 ttGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                                                                             >  4:113881/1‑146 (MQ=255)                                                     cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                                                                 <  1:1017156/138‑1 (MQ=255)                                                     cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                                                                 >  2:1017156/1‑138 (MQ=255)                                                               cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttt                                                                                                        <  3:729758/121‑1 (MQ=255)                                                               cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttt                                                                                                        >  4:729758/1‑121 (MQ=255)                                                                     aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                                                                      <  4:1112791/133‑1 (MQ=255)                                                                     aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                                                                      >  1:640008/1‑133 (MQ=255)                                                                     aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                                                                      <  2:640008/133‑1 (MQ=255)                                                                     aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                                                                      >  3:1112791/1‑133 (MQ=255)                                                                      aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTcc                                                                                                                                               <  4:402432/75‑1 (MQ=255)                                                                      aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTcc                                                                                                                                               >  3:402432/1‑75 (MQ=255)                                                                        gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                                                                   <  2:365358/133‑1 (MQ=255)                                                                        gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                                                                   >  1:365358/1‑133 (MQ=255)                                                                                  tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                                                               >  1:21400/1‑143 (MQ=255)                                                                                  tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                                                               <  2:21400/143‑1 (MQ=255)                                                                                       tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttt                                                                                                        <  2:638430/97‑1 (MQ=255)                                                                                       tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttt                                                                                                        >  1:638430/1‑97 (MQ=255)                                                                                        cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                                                                    >  1:391142/1‑116 (MQ=255)                                                                                        cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                                                                    <  2:391142/116‑1 (MQ=255)                                                                                                                  cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                                                                 >  4:97721/1‑77 (MQ=255)                                                                                                                  cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                                                                 <  3:97721/77‑1 (MQ=255)                                                                                                                               cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                                                                 >  1:979006/1‑96 (MQ=39)                                                                                                                               cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                                                                 <  2:979006/96‑1 (MQ=39)                                                                                                                                 gCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCg                                       >  3:294411/1‑120 (MQ=255)                                                                                                                                 gCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCg                                       <  4:294411/120‑1 (MQ=255)                                                                                                                                  cAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                                                                 >  4:810489/1‑93 (MQ=255)                                                                                                                                  cAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCa                                             <  2:673096/113‑1 (MQ=255)                                                                                                                                    ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                                                                        <  2:349907/68‑1 (MQ=255)                                                                                                                                    ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                                                                        >  1:349907/1‑68 (MQ=255)                                                                                                                                    ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGcc        >  2:96860/1‑148 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                                                                          <  3:639641/80‑1 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                                                                          >  4:639641/1‑80 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc                                               >  2:4599/1‑107 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc                                               <  1:4599/107‑1 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTg                                         <  4:235813/113‑1 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTg                                         >  3:235813/1‑113 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACgg                                <  1:938391/122‑1 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACgg                                >  2:938391/1‑122 (MQ=255)                                                                                                                                      atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGcccc      >  3:490341/1‑148 (MQ=255)                                                                                                                                           atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAc  >  3:976188/1‑147 (MQ=255)                                                                                                                                           atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAc  <  4:976188/147‑1 (MQ=255)                                                                                                                                            tctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                                                                  <  3:751702/66‑1 (MQ=255)                                                                                                                                            tctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                                                                  >  4:751702/1‑66 (MQ=255)                                                                                                                                             |                                                                                                                                                 AATCATCAACGCGATAGATAATATGACAGACATTGAGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTACCAATTGCGTGAACGGTAACATCCAGTGCTGGAGACTGCCCCGCAC  >  NC_000913/2519088‑2519372

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