Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,134,811	Δ1 bp	coding (681/693 nt)	yghT →	putative ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,134,811	0	T	.	98.0%	202.2 / ‑2.8	51	coding (681/693 nt)	yghT	putative ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (23/27);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAA  >  NC_000913/3134689‑3134936                                                                                                                           |                                                                                                                              aTTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCaaa                                                                                                              <  2:1062784/139‑1 (MQ=255)               acaGTTAACCTTTGGTGGTGCACAACTGGTTGTTATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCCCTAATCAAAAATGGAAATGCCCg                                                                                                >  1:824859/1‑139 (MQ=255)               acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                <  1:1059086/139‑1 (MQ=255)               acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                <  1:632468/139‑1 (MQ=255)               acaGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                <  1:552220/139‑1 (MQ=255)                      aCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                <  1:1632003/132‑1 (MQ=255)                      aCCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                >  2:1632003/1‑132 (MQ=255)                          ttGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                                     <  2:671611/139‑1 (MQ=255)                           tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTAt                                                                                                                         >  1:1574781/1‑102 (MQ=255)                           tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTAt                                                                                                                         <  2:1574781/102‑1 (MQ=255)                           tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAcc                                                                                    >  1:1591573/1‑139 (MQ=255)                           tggtggTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAcc                                                                                    >  1:149464/1‑139 (MQ=255)                            ggtggtGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                >  2:453129/1‑126 (MQ=255)                            ggtggtGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                <  1:453129/126‑1 (MQ=255)                              tggtgCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                 <  2:1591573/139‑1 (MQ=255)                                        cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                       >  2:1579158/1‑139 (MQ=255)                                        cTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                       <  2:1382596/139‑1 (MQ=255)                                         tGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGC‑TGGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGa                                                                      <  1:1579158/139‑1 (MQ=255)                                           gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGaa                                                                                                       <  2:571261/104‑1 (MQ=255)                                           gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGaa                                                                                                       >  1:571261/1‑104 (MQ=255)                                           gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                 >  2:701536/1‑126 (MQ=255)                                           gTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCggg                                                                                 <  1:701536/126‑1 (MQ=255)                                                 atCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTaaa                                                              >  1:1230158/1‑139 (MQ=255)                                                   cGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATc                                                            <  2:149464/139‑1 (MQ=255)                                                    gATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCa                                                           >  1:887900/1‑139 (MQ=255)                                                     aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                                     <  1:1027699/112‑1 (MQ=255)                                                     aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGAc                                                                                     >  2:1027699/1‑112 (MQ=255)                                                     aTGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCaa                                                          >  2:1011449/1‑139 (MQ=255)                                                             cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  <  2:1005259/139‑1 (MQ=255)                                                             cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  <  1:1011449/139‑1 (MQ=255)                                                             cAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  <  2:1230158/139‑1 (MQ=255)                                                                  aCTGGAACAGGCGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATc                                             <  2:1189111/139‑1 (MQ=255)                                                                     ggAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                >  2:477642/1‑85 (MQ=255)                                                                     ggAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCg                                                                                                <  1:477642/85‑1 (MQ=255)                                                                            gtgtTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTAt                                              <  2:1506526/128‑1 (MQ=255)                                                                            gtgtTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTAt                                              >  1:1506526/1‑128 (MQ=255)                                                                              gttggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                       >  2:124383/1‑117 (MQ=255)                                                                              gttggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTg                                                       <  1:124383/117‑1 (MQ=255)                                                                                 ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGa                                                                      <  2:80239/99‑1 (MQ=255)                                                                                 ggttgATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGa                                                                      >  1:80239/1‑99 (MQ=255)                                                                                      aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCa                                                                                                                       <  1:268891/45‑1 (MQ=255)                                                                                      aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCa                                                                                                                       >  2:268891/1‑45 (MQ=255)                                                                                      aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  >  1:634355/1‑114 (MQ=255)                                                                                      aTGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  <  2:634355/114‑1 (MQ=255)                                                                                       tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                       >  1:1202782/1‑92 (MQ=255)                                                                                       tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAgg                                                                       <  2:1202782/92‑1 (MQ=255)                                                                                       tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  >  1:336588/1‑113 (MQ=255)                                                                                       tGCAGAAAAAGCGATTACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATg                                                  <  2:336588/113‑1 (MQ=255)                                                                                                      tACGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgt         >  2:61134/1‑139 (MQ=255)                                                                                                        cGGATTTTATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaa       <  1:61134/139‑1 (MQ=255)                                                                                                             tttATGACCGCGCG‑GGTTATCACTAATCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGgtaaagtaa  >  1:458023/1‑139 (MQ=255)                                                                                                                           |                                                                                                                              ATTGCCATAACGCCACAGTTAACCTTTGGTGGTGCACAACTGGTTGATATCGATGCCAATCAGCCACTGGAACAGGTGTTGGTTGATGCAGAAAAAGCGATTACGGATTTTATGACCGCGCGTGGTTATCACTAGTCAAAAATGGAAATGCCCGATCGCCAGGACCGGGCATTTTCAGGAAGGTTAAATCAACTGCAATGCTATCCAGTACAGCCCACCAGAAAGAAAAATTGCCGCCGGTAAAGTAA  >  NC_000913/3134689‑3134936

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