Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,897,782	G→C	9.4%	intergenic (‑182/+52)	PP_5172 ← / ← PP_5173	hypothetical protein/RND family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,897,782	0	G	C	9.4%	131.4 / 3.9	52	intergenic (‑182/+52)	PP_5172/PP_5173	hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base C (2/3);  total (26/26)
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.30e-01

AAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGC  >  NC_002947/5897641‑5897928                                                                                                                                              |                                                                                                                                                    aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                              <  6:221831/149‑1 (MQ=255) aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCA‑TTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCg                                                                                                                                             <  4:377476/149‑1 (MQ=255)  aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTCCCGCGAATACGGTGGTGGGTGGCCCATCg                                                                                                                                             >  8:152355/1‑149 (MQ=255)  aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGGATTAGGTGGTCGATGCACCATCg                                                                                                                                             >  3:20888/1‑149 (MQ=255)               cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTa                                                                                                                                 >  6:90641/1‑148 (MQ=255)               cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTAGCGGGTa                                                                                                                                 >  3:151700/1‑148 (MQ=255)                  gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTCGATGGACCATCGCATTCGCGGGGaaaac                                                                                                                             >  6:50440/1‑147 (MQ=255)                    ttGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                           >  2:270054/1‑149 (MQ=255)                       agtgCAATAGTAGTTATGCGGCATTATTATTAAAAATACTGTTTAAGAATTTTTTTAAAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCt                                                                                                                         <  6:216739/145‑1 (MQ=38)                       aGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                                        <  3:291268/149‑1 (MQ=255)                       aGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCATCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                                        >  5:268167/1‑149 (MQ=255)                           cAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                                    >  6:157749/1‑149 (MQ=255)                            aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGGAAACCCGCTCCTAc                                                                                                                    >  8:101793/1‑148 (MQ=255)                               agaagtCAGGCGGGATAATTCTTAAAAATACTGTTTAAGAATTTTTAAAAAACCATTAAAGACACCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg                                                                                                                <  7:152355/146‑1 (MQ=37)                                gtagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAgg‑g                                                                                                              <  7:101793/149‑1 (MQ=255)                                         cGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATgg                                                                                                     <  5:157749/149‑1 (MQ=255)                                          gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGAGGGTTTACCCGCGAAAACCGTGGTCGAAGCACAATCGCATTAGCCGGTAAACACCATCCTACCAGG‑GGCCATATggg                                                                                                    >  6:136174/1‑149 (MQ=255)                                                  tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAgg                                                                                            >  7:332924/1‑149 (MQ=255)                                                  tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAg                                                                                             <  5:90641/148‑1 (MQ=255)                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCg                                                                                          >  4:284234/1‑149 (MQ=255)                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGGCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCGACCAGG‑GGCCATAGGGGCCTGCAGGCg                                                                                          >  8:249084/1‑149 (MQ=255)                                                      aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCgtg                                                                                        >  2:56772/1‑149 (MQ=255)                                                              tGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCg                                                                                >  3:329940/1‑149 (MQ=255)                                                                  tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                                    <  4:255459/110‑1 (MQ=255)                                                                  tAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTAGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAc                                                                                                                    >  3:255459/1‑110 (MQ=255)                                                                           ttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGCGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATc                                                                     <  1:270054/148‑1 (MQ=255)                                                                            tatataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCtt                                                                   <  1:56772/148‑1 (MQ=255)                                                                                taACCATTTAAGTGAACACGCCTCTGGTAGCAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAcc                                                              <  5:136174/149‑1 (MQ=255)                                                                                  aCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                            >  1:322453/1‑149 (MQ=255)                                                                                      tttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCGGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                            <  6:147379/145‑1 (MQ=255)                                                                                      tttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                            >  5:147379/1‑145 (MQ=255)                                                                                        tAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCAc                                                       <  3:289317/148‑1 (MQ=255)                                                                                          agagAACACGCCGCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                                    <  4:329940/149‑1 (MQ=255)                                                                                               acccGCCTCTGGTAGGTGCGGGTTTACCCGCGATTAGGGTGGTCGATGCACCTTCGCATTCGCGGTTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgcg                                               <  4:151700/146‑1 (MQ=255)                                                                                                    ccTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCgg                                           >  2:238943/1‑148 (MQ=255)                                                                                                      tctGGTAGGAGCGGGTTTTACCGCGAATACGGTGGTCGATGGACCATCGCATTAGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGGCTTCAGGCGTGTTTCTGCGGGTCGCGGATCTTGTACCAGGCCACGTACAGAGCCGGGAg                                        >  8:207145/1‑149 (MQ=37)                                                                                                        tGGTAGGAGCGGTTTACCCCGCTAATACGGTGGTCGAGCCACCATCGAATTCGCGGGAAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                        <  5:39645/147‑1 (MQ=255)                                                                                                              gcgCGGGTTTACCCGCTAATACGGGGGTCCATCCCCCATCGCTTTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagca                                <  2:322453/147‑1 (MQ=38)                                                                                                              gAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                                  >  5:14386/1‑147 (MQ=255)                                                                                                               aGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcagcag                               >  8:101849/1‑149 (MQ=21)                                                                                                                gCGGTTTAACCCGCTAATAGGGTGGCGCATGACCCACGGCTTTGCCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                              <  3:284234/149‑1 (MQ=255)                                                                                                                gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCCTTCGCGGGTTAACCCGCTGCTACCAGG‑GGCCATATGGGGCGTAAGGCGTGTTTATGCGGCGCGCGGATCGTGTACCAGGCCACGTTCAGCGCCGGCAGGCACagcagc                              >  6:56212/1‑149 (MQ=255)                                                                                                                gCGGGTTTACCCGCGAATACGGTGGGCGATGAACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                              <  4:20888/149‑1 (MQ=255)                                                                                                                         cccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAg                     <  8:332924/149‑1 (MQ=255)                                                                                                                          ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg                    >  8:114280/1‑149 (MQ=255)                                                                                                                           cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                                  <  4:220214/134‑1 (MQ=25)                                                                                                                           cgcgAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAAcag                                  >  3:220214/1‑134 (MQ=25)                                                                                                                               aagaCGGTGGTCGACGCACCATCTCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGcc               <  6:268167/146‑1 (MQ=255)                                                                                                                               aaTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGcc               <  7:207145/149‑1 (MQ=255)                                                                                                                                  aCGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAc             >  2:365383/1‑148 (MQ=255)                                                                                                                                            aTGCACGATCGCCTTCTCGGGTAATCCGGCGCCTACCAGG‑GGCCATATGGGACTTCAGGCGTGTTTCTGCGGCTGGCGGATCTTGTACCAGGCCAGGTACAGCGCCGTCAGGATCAGCAGCGTCAGATAGGTCGCCACGATGATCCCGc  <  8:27466/149‑1 (MQ=255)                                                                                                                                            aTGCACAATCGCTTCCGCGGTTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGc  <  5:50440/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                    AAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGG‑GGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGC  >  NC_002947/5897641‑5897928

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