Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	1,744,473	C→T	A67T (GCT→ACT)	apt ←	adenine phosphoribosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,744,473	0	C	T	97.8%	143.0 / NA	46	A67T (GCT→ACT)	apt	adenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (15/30);  total (15/31)
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

TAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAAC  >  CP000730/1744327‑1744616                                                                                                                                                   |                                                                                                                                                tAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTa                                                                                                                                               <  1:691455/149‑1 (MQ=255)  aaCACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGCATAAGTTAc                                                                                                                                              >  1:80059/1‑149 (MQ=255)                 ggTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGAt                                                                                                                               >  1:157752/1‑149 (MQ=255)                         tGCATCTTTGTGCATTGTTACAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCAc                                                                                                                       <  1:1039392/149‑1 (MQ=255)                              ctctgtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGACTGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCt                                                                                                                  <  1:367942/146‑1 (MQ=255)                                    gCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAATCCACGCGCTTCAGGt                                                                                                            <  1:707251/149‑1 (MQ=255)                                    gCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGt                                                                                                            <  1:1020757/149‑1 (MQ=255)                                       ttGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAATAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCt                                                                                                         <  1:724359/149‑1 (MQ=255)                                        tGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTa                                                                                                        <  1:232487/149‑1 (MQ=255)                                          ttAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACa                                                                                                      <  1:8453/149‑1 (MQ=255)                                           tAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCACTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCATaaaa                                                                                                     >  1:968886/1‑149 (MQ=255)                                             aaaCATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACg                                                                                                   <  1:526346/149‑1 (MQ=255)                                                 aTTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGatat                                                                                               >  1:791858/1‑149 (MQ=255)                                                      tACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACa                                                                                          <  1:315955/149‑1 (MQ=255)                                                          atatTCTAGGTCATACTCATAACGAATGACTTTCTGGTGGGATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                                                      >  1:1108864/1‑149 (MQ=255)                                                          atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAACGCCAATCCCCATTGAATACGTTACAGGACAGCCAATGATAAAGCCCCGCGCTTCAGGTCCTCCCACGTTATCAACAtctt                                                                                      >  1:290659/1‑148 (MQ=255)                                                          atatTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                                                      >  1:740570/1‑149 (MQ=255)                                                          atatTCTAGGTCATACTCATAACGAATGACTTCACGAGAAAATTACCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACAtctc                                                                                      >  1:439746/1‑149 (MQ=255)                                                                    tCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATACAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGc                                                                            <  1:876945/149‑1 (MQ=255)                                                                     cATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACTATATCAACATCTCTGTCTTTTGCg                                                                           <  1:626603/149‑1 (MQ=255)                                                                          tCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc                                                                      <  1:458458/149‑1 (MQ=255)                                                                          tCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTc                                                                      <  1:458493/149‑1 (MQ=255)                                                                           cATAACGAATGACTTCACGAGGTCAGTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCt                                                                     <  1:482026/149‑1 (MQ=255)                                                                             tAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCGTCAGATCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTAc                                                                   >  1:762537/1‑149 (MQ=255)                                                                                    tGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGGGTATTCTACAAttttt                                                            >  1:243665/1‑148 (MQ=255)                                                                                        ttCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTg                                                        <  1:297680/149‑1 (MQ=255)                                                                                              aGGTAATTCCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCTACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCAt                                                  <  1:473115/149‑1 (MQ=255)                                                                                               ggTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATa                                                 <  1:49895/149‑1 (MQ=255)                                                                                                      tCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCatat                                          >  1:998471/1‑149 (MQ=255)                                                                                                      tCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAAAGTACTCAGCATCTCTGTCTTTTGCGTATTCTACGATTTTATCTGTTGCATAGCCatat                                          >  1:1049196/1‑149 (MQ=255)                                                                                                       cccTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACACTTTTATCTGTTGCATAGCCATATg                                         >  1:170444/1‑149 (MQ=255)                                                                                                              tctCTAACAGGTGCAAACGCAAGCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCCCGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCAcc                                  <  1:860329/147‑1 (MQ=255)                                                                                                                    aCAGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATc                            <  1:465534/149‑1 (MQ=255)                                                                                                                     cAGGTGCAAAGCCAATCCCCATTTAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATcc                           <  1:264342/149‑1 (MQ=255)                                                                                                                      aGGTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCa                          <  1:138816/149‑1 (MQ=255)                                                                                                                      aGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGAACAAACAACGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGGTGCATAGCCATATGCTTCACCATTATCCa                          <  1:585779/149‑1 (MQ=255)                                                                                                                       ggTGCAAAGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTTTTCTACAATTTTGTCTGTTGCGTAGCCATTTGAttgtcctttttcttt                         >  1:106639/1‑136 (MQ=255)                                                                                                                          gCAAACCCACTCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAAt                      <  1:235216/149‑1 (MQ=255)                                                                                                                            aaaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACCACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCTCCATTTTTCATACttt                    >  1:352427/1‑148 (MQ=255)                                                                                                                             aaGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTgt                   <  1:417504/149‑1 (MQ=255)                                                                                                                              aGCCAATCCCCATTGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTgta                  <  1:159606/149‑1 (MQ=255)                                                                                                                                          ttGAATAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTg      <  1:272739/149‑1 (MQ=255)                                                                                                                                             aaTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGaaa   <  1:17580/149‑1 (MQ=255)                                                                                                                                             aaTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGaaa   <  1:186735/149‑1 (MQ=255)                                                                                                                                             aaTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGaaa   <  1:789292/149‑1 (MQ=255)                                                                                                                                              aTAAGTTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAAc  <  1:845542/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACGCGCTTCAGGTCCTACAACGATATCAACATCTCTGTCTTTTGCGTATTCTACAATTTTATCTGTTGCATAGCCATATGCTTCACCATTATCCATAATTGTAGTAATATCCTTGAAAC  >  CP000730/1744327‑1744616

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