Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,809,905	G→T	9.7%	P129T (CCT→ACT)	PP_5086 ←	nuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,809,905	0	G	T	9.7%	117.6 / 2.8	41	P129T (CCT→ACT)	PP_5086	nuclease
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base T (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

GATACCTTCGATACGACCACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCACCAGCCCGACGCGATCGTCATTGGCCT  >  NC_002947/5809756‑5810054                                                                                                                                                      |                                                                                                                                                      gATACCTTCGATACGACCACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAgg                                                                                                                                                       <  2:98200/150‑1 (MQ=255) gATACCTTCGATACGACCACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAgg                                                                                                                                                       <  2:156158/150‑1 (MQ=255)    aCCTTCGATACGACCACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGc                                                                                                                                                    >  2:580497/1‑150 (MQ=255)        tCGATACGACCACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATa                                                                                                                                                <  2:674876/150‑1 (MQ=255)                 ccACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGCCAATAGCCATGCgg                                                                                                                                       <  2:593121/150‑1 (MQ=255)                   gcGCACAACAGCGAAAGCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACAGTAGGGGCAATAGCCACGCGGTa                                                                                                                                     <  2:41392/149‑1 (MQ=255)                   aCGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTa                                                                                                                                     >  2:520457/1‑150 (MQ=255)                         aacaGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCaa                                                                                                                               >  2:379004/1‑150 (MQ=255)                                    cTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAgg                                                                                                                                                       >  2:637696/1‑115 (MQ=255)                                    cTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAgg                                                                                                                                                       <  1:637696/115‑1 (MQ=255)                                    cTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTaa                                                                                                                    >  2:562128/1‑150 (MQ=255)                                           ccTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTACTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTg                                                                                                             >  1:462971/1‑150 (MQ=255)                                             tGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTgcg                                                                                                           >  1:170116/1‑150 (MQ=255)                                              gACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTgcgc                                                                                                          <  2:134347/150‑1 (MQ=255)                                               aCATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCt                                                                                                         >  1:121664/1‑150 (MQ=255)                                               aCATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCt                                                                                                         <  1:259919/150‑1 (MQ=255)                                                cATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCtt                                                                                                        >  1:335318/1‑150 (MQ=255)                                                          aGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTc                                                                                              <  2:468481/150‑1 (MQ=255)                                                            ccACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGc                                                                                            <  2:170116/150‑1 (MQ=255)                                                                tACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGTGGCaa                                                                                                                                                  <  2:520073/92‑1 (MQ=255)                                                                tACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGTGGCaa                                                                                                                                                  >  1:520073/1‑92 (MQ=255)                                                                                 acaGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTgcg                                                                       >  1:207196/1‑150 (MQ=255)                                                                                       cAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAgtg                                                                 >  2:655660/1‑150 (MQ=255)                                                                                                 cTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTa                                                       <  2:55097/150‑1 (MQ=255)                                                                                                  tGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCAtt                                                                                        <  1:407620/116‑1 (MQ=255)                                                                                                  tGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCAtt                                                                                        >  2:407620/1‑116 (MQ=255)                                                                                                         cAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCtt                                               >  1:598436/1‑150 (MQ=255)                                                                                                          aGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCttt                                              <  2:121664/150‑1 (MQ=255)                                                                                                               gTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGTGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGt                                                  >  2:330158/1‑141 (MQ=255)                                                                                                               gTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGTGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGt                                                  <  1:330158/141‑1 (MQ=255)                                                                                                               gTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGCTTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCtt                                         <  1:580497/150‑1 (MQ=255)                                                                                                               gTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGCAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCtt                                         <  1:562128/150‑1 (MQ=255)                                                                                                                   ggCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTa                                                                               <  2:663420/108‑1 (MQ=255)                                                                                                                   ggCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTa                                                                               >  1:663420/1‑108 (MQ=255)                                                                                                                   ggCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAac                                     <  2:335318/150‑1 (MQ=255)                                                                                                                       aGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCa                                 <  2:323043/150‑1 (MQ=255)                                                                                                                          tgctgACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGAt                                                                            >  2:677027/1‑104 (MQ=255)                                                                                                                          tgctgACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGAt                                                                            <  1:677027/104‑1 (MQ=255)                                                                                                                            ctgACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCAc                            >  1:481004/1‑150 (MQ=255)                                                                                                                                       tGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCACCAGCCCGAcgc                 >  1:417887/1‑150 (MQ=255)                                                                                                                                             gTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCAtt                                                                                        <  2:619123/73‑1 (MQ=255)                                                                                                                                             gTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCAtt                                                                                        >  1:619123/1‑73 (MQ=255)                                                                                                                                             gTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCACCAGCCCGACGCGATCGt           >  1:594614/1‑150 (MQ=255)                                                                                                                                                 aTTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCa                                 >  2:512466/1‑124 (MQ=255)                                                                                                                                                 aTTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCCACCGTACTTGTCTTTACCTTCAACACCa                                 <  1:512466/124‑1 (MQ=255)                                                                                                                                                      gggCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCCGATTGTCGCCATTGCGGCTGTAGCTGCGCGGCAGTGTCCCACCTTCCTTGTCTTTACCTTCACAGCCCGGCACCAGCACGCCGCGATCGTCATTGGCCt  >  1:317459/1‑150 (MQ=255)                                                                                                                                                      gggCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCACCAGcc                       >  2:155364/1‑129 (MQ=255)                                                                                                                                                      gggCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCACCAGcc                       <  1:155364/129‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GATACCTTCGATACGACCACGCACAACAGCGAAACCTGACTGCCTGACATCGCTTGCAGCCACTACGGGCGAATGCCGCCACAGGCCAACGCCGGCCTGGCGCGCAGCTTGTTCGGCAAGCTGCTGACAACCGCTGAGACGTACATTAGGGGCAATAGCCACGCGGTAGCCCAAGCCCTCGCTAAGCAATTGCGCTTCCAGATTGTCGCCATTGCGGCTGTAGATGTGCGCCAGTGTGCGACCGTACTTGTCTTTACCTTCAACACCAGGCACCAGCCCGACGCGATCGTCATTGGCCT  >  NC_002947/5809756‑5810054

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