Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,260,849	G→A	intergenic (‑46/+79)	dauA ← / ← prs	C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,260,849	0	G	A	98.0%	149.9 / ‑5.2	49	intergenic (‑46/+79)	dauA/prs	C4‑dicarboxylic acid transporter/phosphoribosylpyrophosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (26/22);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

CAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAC  >  NC_000913/1260701‑1260995                                                                                                                                                     |                                                                                                                                                   cAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGt                                                                                                                                                  <  1:649755/151‑1 (MQ=255)   aTCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTg                                                                                                                                                 >  1:606453/1‑150 (MQ=255)       ggTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAagg                                                                                                                                             >  1:512668/1‑150 (MQ=255)           aCGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCa                                                                                                                                        <  1:104709/151‑1 (MQ=255)             gggTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATa                                                                                                                                      <  1:443817/151‑1 (MQ=255)              ggTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAg                                                                                                                                     <  1:474249/151‑1 (MQ=255)                    ccGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCAt                                                                                                                                >  1:7751/1‑150 (MQ=255)                               tataTTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGatat                                                                                                                    >  2:514814/1‑151 (MQ=255)                                    ttttCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACa                                                                                                                <  2:678690/150‑1 (MQ=255)                                    ttttCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCACATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACa                                                                                                                <  2:256097/150‑1 (MQ=255)                                     tttctttcCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTAc                                                                                                                 >  1:625196/1‑148 (MQ=255)                                      ttctttcCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGAc                                                                                                             >  2:228135/1‑151 (MQ=255)                                         tttcCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACaa                                                                                                           >  2:145866/1‑150 (MQ=255)                                            ccAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCGTAGGTCATACAAATGGATATTACAGACaaaaa                                                                                                        <  2:462841/150‑1 (MQ=255)                                                 atgCGTTGATCAGAGCGCGGAAAGGCATCACATGTCAGAATAATATTTTGTTCACAATAATGTTTCAACCGTGAGCGCAACAGCTGTGAACTAAATGGCAAGTGAAGGGGGCATAGGTCATACAAATGTATACTTCAGGCAAGAAAACccg                                                                                                  <  1:38628/149‑1 (MQ=255)                                                 aaGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccg                                                                                                  >  1:278025/1‑151 (MQ=255)                                                  aGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACccgc                                                                                                 >  1:32776/1‑151 (MQ=255)                                                      tCGATCAGAGCGCGGAAAGGCATCACATGTGAGAAAAATATTTTGTTCAAAATAATGTTTCATCCGTGAGCGCATCATCTGTAAACTAAATGGAAAGTGAAGGAGGCATAGGTCATACAAAGGGATATTACAGAAAAAAAACCccgccg                                                                                               <  1:533944/149‑1 (MQ=255)                                                        gATCAGAGCGCGGAAAGGCATCACATGTGAGAAAAATATTTTGTTCAAAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATGGGTCATACAAATGGATATTACAGACAAAAAACCCCGCCGCAg                                                                                            <  1:140687/150‑1 (MQ=255)                                                                cgcgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCttt                                                                                   >  1:261507/1‑151 (MQ=255)                                                                cgcgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCttt                                                                                   >  1:178201/1‑151 (MQ=255)                                                                  cgGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                  >  2:607905/1‑150 (MQ=255)                                                                   ggAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTg                                                                                  >  1:580125/1‑149 (MQ=255)                                                                     aaaGGCATCACATGTGAGAAAAATATTTTGTTCAAAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACAAAAGGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGc                                                                               <  1:102928/150‑1 (MQ=255)                                                                              acaTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAAGGGCAAGTGAAGGAGGCATAGGTCATAAAAATGGATATTACAGACAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAt                                                                     <  2:42649/151‑1 (MQ=255)                                                                               caTGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGAtt                                                                    >  2:453898/1‑151 (MQ=255)                                                                                      ggAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTc                                                                        >  1:25727/1‑140 (MQ=255)                                                                                        aaaaTATTTTGTTCACAATAATGTTTCACCCGTGAGCGCATCATCTGTCAACTAAATGGAAAGGGAAGGAGGCATAGGTCATACAAAGGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTc                                                                        <  2:25758/138‑1 (MQ=255)                                                                                        aaaaTATTTTGTTCAAAATAAGGTTTCACCAGGGAGCCCATCATCTTTCAACAAAATGGAAAGTAAAGGAGGCAAGGGTCATACAAATGGATATTACAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCg                                                            <  1:276447/150‑1 (MQ=255)                                                                                              ttttGTACAAAAAAATGTTTCATCCGTGAGCAAATCATCTGTCAACTAAAGGAAAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAAAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATg                                                      <  1:521981/150‑1 (MQ=255)                                                                                                   ttCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCag                                                  >  1:455786/1‑149 (MQ=255)                                                                                                   ttCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATAAAAATGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCag                                                  <  2:512720/149‑1 (MQ=255)                                                                                                    tCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCag                                                  >  2:218991/1‑148 (MQ=255)                                                                                                      acaATAATGTTTCTTCGGTGAGCGGATCATCTGTCAACTAAATGGCAAGTGAAGGAGGAATGGGTCATAAAAATGGATATTAAAAAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCagag                                                <  2:429321/148‑1 (MQ=255)                                                                                                        aataatGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCg                                            >  2:534212/1‑150 (MQ=255)                                                                                                                ttCACCGGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAGGGAGGCATAGGTAATACAAAGGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCg                                    <  1:2631/150‑1 (MQ=255)                                                                                                                 tCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAATAAACCCGCCGCAGCGCTTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATGGATTCTTCGtt                                  >  1:314308/1‑151 (MQ=255)                                                                                                                    tCCGTGAGCGCATCATCTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                               >  2:541081/1‑151 (MQ=255)                                                                                                                    cccGTGAGCGAATCATCTGTCAACTAAAGAGAAGATGAGGGAGGCATAGGTCATACAAAGGGATATTAAAGACAAAAAAGCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCt                               <  1:16016/150‑1 (MQ=255)                                                                                                                                tcatcTGTCAACAAATTGGCAAGTGAGGGAGCCATAGGTCATACAAGTGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATacgac                       <  2:444581/147‑1 (MQ=255)                                                                                                                                  atcTGTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCg                 >  1:160109/1‑151 (MQ=255)                                                                                                                                      gTCAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCt               <  2:367446/149‑1 (MQ=255)                                                                                                                                        cAACTAAATGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCgg           >  2:47832/1‑151 (MQ=255)                                                                                                                                         aaCTAAATGGAAAGTGAGGGAGGCATAGGTCATACAAAGGGATATTACAGAAAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGc          <  1:181265/151‑1 (MQ=255)                                                                                                                                         aaCTAAAGGGAAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGAAAAAAAACCCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCg            <  1:607848/149‑1 (MQ=255)                                                                                                                                               aTGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCAt    >  2:244051/1‑151 (MQ=255)                                                                                                                                                tGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCAt    >  2:54030/1‑150 (MQ=255)                                                                                                                                                tGGCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCAt    >  2:710076/1‑150 (MQ=255)                                                                                                                                                 ggCAAGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAc  >  2:261825/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CAATCAGGTCACGGGTAAACCGTGCGGCAGTATATTTTTCTTTCCAGCAAGCGTCGATCAGAGCGCGGAAAGGCATCACATGTGAGGAAAATATTTTGTTCACAATAATGTTTCATCCGTGAGCGCATCATCTGTCAACTAAATGGCAGGTGAAGGAGGCATAGGTCATACAAATGGATATTACAGACAAAAAAACCCGCCGCAGCGGGTCTTTGAGCCGGGTTCGATTAGTGTTCGAACATGGCAGAGATCGATTCTTCGTTGCTGATACGACGAATCGCTTCGGCCAGCATAC  >  NC_000913/1260701‑1260995

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