Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,625,537	T→A	9.7%	intergenic (‑50/+27)	PP_3195 ← / ← PP_3196	hypothetical protein/nodulation protein d1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,625,537	0	T	A	9.7%	108.7 / 3.2	41	intergenic (‑50/+27)	PP_3195/PP_3196	hypothetical protein/nodulation protein d1
Reads supporting (aligned to +/- strand):  ref base T (22/15);  new base A (2/2);  total (24/17)
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.14e-01

GTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAATT  >  NC_002947/3625390‑3625685                                                                                                                                                     |                                                                                                                                                     gtCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCt                                                                                                                                                      <  8:350901/148‑1 (MQ=255) gtCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTg                                                                                                                                                     >  7:124518/1‑149 (MQ=255)             tCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAGACACCGCAGCGGTG‑CTgggg                                                                                                                                         <  1:132364/149‑1 (MQ=255)                            cccACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCt                                                                                                                           >  8:180433/1‑148 (MQ=255)                            cccACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTc                                                                                                                          <  1:230708/149‑1 (MQ=255)                            cccACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGAACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTc                                                                                                                          <  2:321004/149‑1 (MQ=255)                                  tGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCgg                                                                                                                    >  6:107672/1‑149 (MQ=255)                                  tGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCg                                                                                                                     <  8:124518/148‑1 (MQ=255)                                        aCAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgcgc                                                                                                              <  3:118213/149‑1 (MQ=255)                                        aCAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGgcgc                                                                                                              <  7:180433/149‑1 (MQ=255)                                                       gcgGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCAGGCt                                                                                               >  3:188416/1‑149 (MQ=255)                                                                       gTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACt                                                                                >  4:239452/1‑148 (MQ=255)                                                                           gcaggcagGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAgg                                                                                                                       <  8:262952/105‑1 (MQ=255)                                                                           gcaggcagGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAgg                                                                                                                       >  7:262952/1‑105 (MQ=255)                                                                           gcaggcagGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATg                                                                           <  5:107672/149‑1 (MQ=255)                                                                               gcagggcaggGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTCTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAt                                                                       >  7:216862/1‑149 (MQ=255)                                                                                         gTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGAc                                                                         <  4:218539/137‑1 (MQ=255)                                                                                         gTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGAc                                                                         >  3:218539/1‑137 (MQ=255)                                                                                              gCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCTTCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTcgc                                                        >  5:263145/1‑149 (MQ=255)                                                                                              gCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTcgc                                                        >  4:232894/1‑149 (MQ=255)                                                                                                  ttcttGATCTGTGATGCGTGGTATCGCGAAACGCGATGCCGCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGGCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAg                                                    <  4:188416/145‑1 (MQ=255)                                                                                                                  cGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCTTCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATggg                                     <  6:263145/148‑1 (MQ=255)                                                                                                                   gAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCCGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGt                                   >  4:294254/1‑149 (MQ=255)                                                                                                                   gAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCAGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGt                                   >  2:163487/1‑149 (MQ=255)                                                                                                                      gTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGc                                >  4:344078/1‑149 (MQ=255)                                                                                                                        aTCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACAACGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCg                              >  5:325976/1‑149 (MQ=255)                                                                                                                         tCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGt                             <  3:344078/149‑1 (MQ=255)                                                                                                                         tCGCGAAACGCGATGCCCCAG‑CACCGCTGCCGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGt                             >  1:30670/1‑149 (MQ=255)                                                                                                                           gcgAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTcc                           <  8:231364/149‑1 (MQ=255)                                                                                                                           gcgAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTc                            >  7:231364/1‑148 (MQ=255)                                                                                                                            cgAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCCCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTcc                           >  4:309054/1‑148 (MQ=255)                                                                                                                                 cgcgATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCAtgctg                     >  7:276256/1‑149 (MQ=255)                                                                                                                                    atcgCCCCAGACACCGCAGCGGTG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCgg                  >  1:63577/4‑149 (MQ=255)                                                                                                                                        ccccAGACACCGCAGCGGCG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGc               >  8:153295/1‑148 (MQ=255)                                                                                                                                         cccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCca             >  5:52926/1‑149 (MQ=255)                                                                                                                                          ccAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCca             >  1:304593/1‑148 (MQ=255)                                                                                                                                              acACCGCAGCGGCG‑CTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCa        <  7:153295/148‑1 (MQ=255)                                                                                                                                               cACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAAt                                                                            >  6:137777/1‑81 (MQ=255)                                                                                                                                               cACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAAt                                                                            <  5:137777/81‑1 (MQ=255)                                                                                                                                                aCCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCCGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCAt       <  3:294254/149‑1 (MQ=255)                                                                                                                                                  cGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGa     <  2:30670/149‑1 (MQ=255)                                                                                                                                                    cTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAAt   >  7:310719/1‑149 (MQ=255)                                                                                                                                                     tGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAAtt  <  8:276256/149‑1 (MQ=255)                                                                                                                                                     tGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAAt   >  3:31600/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GTCCTTCTGCGGTCGTGGAGTCCGGCGCCCACTTGAGCGACAGGCAGGCAGGGCGCGGTTTGAAGGCAAGGTACGCAGGCAGGGCAGGGTGGGGCCAATCGTGATCAGGGATGCGAGGTATCGCGAAACGCGATGCCCCAG‑CACCGCTGCGGTGTCTGGGGCGATGGCCGCCGCTCAGGCGGTGGCGCGCAGAGCATCTGGCTGCGTGCTGGCGCACTGAATGACATCGCTGATGCGCTCGCGAAGCCATCTGTGCATGGGGTCGCCGTCCATGCTGCGGTGCCACAGCATGAATT  >  NC_002947/3625390‑3625685

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