Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	1,386,769	T→C	8.2%	H236R (CAT→CGT)	U712_07095 ←	Methylthioribose kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	1,386,769	0	T	C	8.2%	135.0 / 4.2	49	H236R (CAT→CGT)	U712_07095	Methylthioribose kinase
Reads supporting (aligned to +/- strand):  ref base T (25/20);  new base C (2/2);  total (27/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 = 9.44e-01

GCCCCGTCACGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGAAGCTCTTCTTCAAAATCATTTGTGTCAT  >  CP006881/1386634‑1386907                                                                                                                                        |                                                                                                                                           gCCCCGTCACGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGACCTACATCAAAGCCGATCTGGCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTAGTGCTCGCTTGCGAAAATGCTTACGGTAttc                                                                                                                                          >  5:286047/1‑136 (MQ=255)   cccGTCACGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCaaa                                                                                                                                       >  3:90928/1‑139 (MQ=255)   cccGTCACGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCaaa                                                                                                                                       <  8:194867/139‑1 (MQ=255)    ccGTCACGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAAt                                                                                                                                      <  4:311547/139‑1 (MQ=255)         tcGGATGAGTGCGTTAATAAATAAGTTGTCGATATATTGCTCTACATAAAAGCAGATCGGTACGTAAAATTCAATCTCTGTATCTATCAACTTTGTTTGGTGATCACTTGCGAAAATTCTTCCTGTATGCAAATctcg                                                                                                                                  <  6:286047/137‑4 (MQ=255)         aCGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCAcc                                                                                                                                  <  8:5091/138‑1 (MQ=255)         aCGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCAcc                                                                                                                                  <  2:318225/138‑1 (MQ=255)                aGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTACGCATATCACCATGAATc                                                                                                                           >  8:291808/1‑138 (MQ=255)                aGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTACGCATATCACCATGAATc                                                                                                                           >  8:287056/1‑138 (MQ=255)                 gTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAATCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCaa                                                                                                                         <  2:141165/139‑1 (MQ=255)                    cGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAACGCCGATCGGTCCGTAAAACGCAAACTCAGGAGCAATCCCCTTTGTTTCGTGCTCACTTGCGACAATGCTTCCTGTATGCAAATCACCAAGCATCCACGt                                                                                                                      >  7:4615/1‑139 (MQ=255)                       tCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAAt                                                                                                                                      <  6:297567/120‑1 (MQ=255)                       tCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAAt                                                                                                                                      >  5:297567/1‑120 (MQ=255)                           aaaTAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCa                                                                                                               <  7:345430/139‑1 (MQ=255)                                 gTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAATATGCTTCCTGTATGCAAATCACCATGAATAAACGTTTCTGCAGACGt                                                                                                          >  1:255050/1‑138 (MQ=255)                                           aaTTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGAc                                                                                                <  3:73234/138‑1 (MQ=255)                                           aaTTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTACGCATATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGAc                                                                                                <  7:287056/138‑1 (MQ=255)                                           aaTTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTACGCATATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGAc                                                                                                <  7:291808/138‑1 (MQ=255)                                                cccTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACttttt                                                                                           <  7:291528/138‑1 (MQ=255)                                                     cATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCtt                                                                                     >  1:339397/1‑139 (MQ=255)                                                          aaGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCgg                                                                                >  1:269986/1‑139 (MQ=255)                                                          aaGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCgg                                                                                >  2:97649/1‑139 (MQ=255)                                                           aGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATC‑ACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCt                                                                              <  4:278757/139‑1 (MQ=255)                                                                aTCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGa                                                                          >  8:258281/1‑139 (MQ=255)                                                                aTCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGa                                                                          >  5:341437/1‑139 (MQ=255)                                                                aTCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGa                                                                          >  2:19569/1‑139 (MQ=255)                                                                aTCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGa                                                                          >  1:270519/1‑139 (MQ=255)                                                                aTCGGTCCGTAAAAAGCAAAATCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCTAAAATGCTTCCTGAATGCAAATCACCATGAATCAACGTATCTGCAGACTTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGa                                                                          >  1:283617/1‑139 (MQ=255)                                                                      ccGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTaa                                                                    <  8:286183/139‑1 (MQ=255)                                                                       cGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAAc                                                                   >  3:208050/1‑139 (MQ=255)                                                                         tAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGc                                                                 <  1:19569/139‑1 (MQ=255)                                                                          aaaaaGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCt                                                                >  1:118235/1‑139 (MQ=255)                                                                           aaaaGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCt                                                                <  1:122596/138‑1 (MQ=255)                                                                           aaaaGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCt                                                                <  1:97649/138‑1 (MQ=255)                                                                               gCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCat                                                           <  8:111663/139‑1 (MQ=255)                                                                                                          tCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACAATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTAATTATTCCAAATCTTTTCGACAAAAGGa                                 <  7:72967/138‑1 (MQ=255)                                                                                                          tCGTGCTCAATTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGa                                 <  5:343341/138‑1 (MQ=255)                                                                                                            gTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACg                               >  3:3201/1‑138 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGc                                                                 <  4:300137/103‑1 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGc                                                                 >  3:300137/1‑103 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACg                               >  1:170123/1‑137 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACg                               >  3:176181/1‑137 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGaa                             >  4:39256/1‑139 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGaa                             >  5:215588/1‑139 (MQ=255)                                                                                                             tGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCAAGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACg                               >  3:349181/1‑137 (MQ=255)                                                                                                                  aCTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGAc                                                                                                >  7:311053/1‑67 (MQ=255)                                                                                                                  aCTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGAc                                                                                                <  8:311053/67‑1 (MQ=255)                                                                                                                       cGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGAAGCTcttctt                    >  3:262139/1‑138 (MQ=255)                                                                                                                            aTGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGAAGCTCTTCTTCaaaa               >  6:82134/1‑138 (MQ=255)                                                                                                                                 tCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGAAGCTCTTCTTCAAAATCAttt         >  3:53926/1‑139 (MQ=255)                                                                                                                                        tGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAACGCTTTTCGACAAAAGGACGAAGCTCTTCTTCAAAATCATTTGTGTCAt  >  2:180691/1‑139 (MQ=255)                                                                                                                                        tGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGTCAAAAGGACGAAGCTCTTCTTCAAAATCATTTGTGTCa   <  6:307092/138‑1 (MQ=255)                                                                                                                                        tGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGAAGCTCTTCTTCAAAATCATTTGTGTCa   <  2:66004/138‑1 (MQ=255)                                                                                                                                        |                                                                                                                                           GCCCCGTCACGGCTGAGTGCGTTCAAAAATAAGTTGGCGATAAATTGCCCTACATCAAAGCCGATCGGTCCGTAAAAAGCAAACTCAGGATCAATCACCTTTGTTTCGTGCTCACTTGCGAAAATGCTTCCTGTATGCAAATCACCATGAATCAACGTTTCTGCAGACGTTAAAAATGACTTTTTCAGCTTAGCGGCTTCGATTTTAACGCTGTCATTATTCCAAAGCTTTTCGACAAAAGGACGAAGCTCTTCTTCAAAATCATTTGTGTCAT  >  CP006881/1386634‑1386907

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