Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,430,426	C→T	9.6%	intergenic (‑108/+53)	PP_4769 ← / ← PP_4770	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,430,426	0	C	T	9.6%	143.2 / 4.7	52	intergenic (‑108/+53)	PP_4769/PP_4770	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/22);  new base T (3/2);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCACTGCCCCAGCCACCGTT  >  NC_002947/5430288‑5430574                                                                                                                                            |                                                                                                                                                     gCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                                            >  3:32824/1‑149 (MQ=255)  cacaTGGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCGGCTCCCACAGGTAT‑CGCGCGTACACCGGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                                              <  2:62909/146‑1 (MQ=255)  cacaTCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAg                                                                                                                                              <  2:58523/146‑1 (MQ=255)                 cgcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCg                                                                                                                            >  1:99546/1‑149 (MQ=255)                  gcgATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCgg                                                                                                                           <  2:99546/149‑1 (MQ=255)                    gATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCCCAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAg                                                                                                                                    >  4:261504/1‑138 (MQ=255)                    gATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAg                                                                                                                                    <  3:261504/138‑1 (MQ=255)                        tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                                                                       <  5:126831/147‑1 (MQ=255)                        tttttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAAt                                                                                                                       >  6:126831/1‑147 (MQ=255)                          tttCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGaaa                                                                                                                        <  3:131941/144‑1 (MQ=255)                                  cgagctcTCCTTTTGGGCCTGGAGGGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCCCAGGTCT‑CGCGCGGACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATg                                                                                                           <  4:103367/143‑1 (MQ=38)                                   ggAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGc                                                                                                          <  4:32824/149‑1 (MQ=255)                                        cTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGCCACCCCCGCAACCCCACGTAT‑AGCCCGGACACCCCTGGGCCCCGGCGGTCCCGCCACCGGGCGCAAACCGCCCCCCGCAATGCCCGCAGAATTcccgcc                                                                                                     >  6:19223/1‑145 (MQ=255)                                         tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAcc                                                                                                    >  2:170708/1‑149 (MQ=255)                                         tCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAc                                                                                                     >  8:113565/1‑148 (MQ=255)                                                  gCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACg                                                                                           >  7:245446/1‑149 (MQ=255)                                                   ccTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACg                                                                                           >  7:46635/1‑148 (MQ=255)                                                   ccTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCCCGCCCCCAGGTAT‑AGCGCGTACACCCGTGGGAGCGGGGGTGCCCGCCCCGGGGCGCAAAGCGCCCCCGGCAATGCCAGCAGTATGCCa                                                                                                        >  8:158295/1‑135 (MQ=255)                                                    cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCGGTATGCCAGACCTGTCAGACGt                                                                                          >  5:135267/1‑148 (MQ=255)                                                    cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGt                                                                                          >  2:4456/1‑148 (MQ=255)                                                    cTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGt                                                                                          >  2:200836/1‑148 (MQ=255)                                                          tGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtc                                                                                   <  6:135267/149‑1 (MQ=255)                                                            ccGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcg                                                                                  <  1:200836/148‑1 (MQ=255)                                                             cGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                                                <  1:170708/149‑1 (MQ=255)                                                             cGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAtcgtc                                                                                <  1:4456/149‑1 (MQ=255)                                                                           gcAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                                                                     >  6:54595/1‑98 (MQ=255)                                                                           gcAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGc                                                                                                                     <  5:54595/98‑1 (MQ=255)                                                                                 cTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGt                                                            >  2:20308/1‑149 (MQ=255)                                                                                  ttCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGt                                                            >  7:242355/1‑148 (MQ=255)                                                                                    cgcgGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTc                                                           <  4:779/147‑1 (MQ=255)                                                                                         acacGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTg                                                    <  8:46635/149‑1 (MQ=255)                                                                                          cacGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGc                                                   >  1:31265/1‑149 (MQ=255)                                                                                             gCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAAGCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCa                                                 >  5:231435/1‑148 (MQ=255)                                                                                                 gCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                                                            >  7:109348/1‑117 (MQ=21)                                                                                                 gCTCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTcc                                                                            <  8:109348/117‑1 (MQ=21)                                                                                                   tCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAaccac                                          >  2:112533/1‑149 (MQ=37)                                                                                                   tCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCt                                                                                                   >  8:138819/1‑92 (MQ=18)                                                                                                   tCCCACGGGTGTACGCGCG‑ATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCt                                                                                                   <  7:138819/92‑1 (MQ=18)                                                                                                       acaGGTAT‑CGCGCGTACACCCGGGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGGCATCGCTGCCAAAACCACCGTc                                      <  5:19223/149‑1 (MQ=255)                                                                                                         aGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCgg                                    <  8:245446/149‑1 (MQ=255)                                                                                                         aGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCgg                                    >  2:191325/1‑149 (MQ=255)                                                                                                                  gcgcgatacCTGTGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCg                            <  1:112533/142‑1 (MQ=255)                                                                                                                     gTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                                        <  4:255809/86‑1 (MQ=255)                                                                                                                     gTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTaa                                                                                        >  3:255809/1‑86 (MQ=255)                                                                                                                     gTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGccc                         <  8:242355/149‑1 (MQ=255)                                                                                                                           ccGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCa                   <  7:113565/149‑1 (MQ=255)                                                                                                                           ccGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCa                   >  2:75424/1‑149 (MQ=255)                                                                                                                            cGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCAc                  >  3:17545/1‑149 (MQ=255)                                                                                                                              tGGGAGCGGGCGTGTCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGACCCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCACTg                >  3:158162/1‑149 (MQ=255)                                                                                                                              tGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCACTg                >  6:152194/1‑149 (MQ=255)                                                                                                                                 gAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCCGAAATGCCAGCAGTATGCCAGCGCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCCCCGTCATCGCTGCCAAACCCACCGTCGGCGTATCCGCCCTGGACACTGccc             >  7:66013/1‑149 (MQ=255)                                                                                                                                         gTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGAACTGTCAGACGTaa                                                                                        <  4:201383/66‑1 (MQ=255)                                                                                                                                            cccGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCACTGCCCCAGCCACCGtt  >  5:217257/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                     GCACATCGGCGAACAGCGCGATCTTTTTCACCTGGGAACCTCATTTTTGGCCTGGAGTGCCGGGGGGCTGCTGCGCAGACCTTCGCGGACACGCCCGCTCCCACAGGTAT‑CGCGCGTACACCCGTGGGAGCGGGCGTGCCCGCGAAGGGGCGCAAAGCGCCCCCGGAAATGCCAGCAGTATGCCAGACCTGTCAGACGTAAGAATCGTCGTCCGAGAAGAACCCACCGTCATCGCTGCCAAAACCACCGTCGGCGTATCCGCCCTGGTCACTGCCCCAGCCACCGTT  >  NC_002947/5430288‑5430574

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