Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,897,791	C→T	11.1%	intergenic (‑191/+43)	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,791	0	C	T	11.1%	120.3 / 4.9	45	intergenic (‑191/+43)	PP_5172/PP_5173	hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base T (2/3);  total (21/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 = 7.01e-01

GGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCAATCATGGCG  >  NC_002947/5897644‑5897939                                                                                                                                                    |                                                                                                                                                      gggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                                      >  5:130870/1‑149 (MQ=255)     tttGCGCCTGGCTTGAGGGCAATATTAGTGAGGAGGTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTc                                                                                                                                                   <  8:155423/148‑1 (MQ=255)          gcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                                              >  4:193934/1‑148 (MQ=255)            cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                                 <  8:41822/143‑1 (MQ=25)            cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                                 >  7:41822/1‑143 (MQ=25)               gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                                         <  1:19959/148‑1 (MQ=255)                  tGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAAGACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCAGGGTATTGGCGGGTAAACCCg                                                                                                                                      <  4:150442/148‑1 (MQ=25)                   gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGt                                                                                                                                             >  1:165627/1‑140 (MQ=255)                   gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGt                                                                                                                                             <  2:165627/140‑1 (MQ=255)                      ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                                 >  1:24742/1‑133 (MQ=25)                      ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                                 <  2:24742/133‑1 (MQ=25)                           tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAg                                                                                                                           >  5:253882/1‑149 (MQ=255)                                 tGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCa                                                                                                                     >  1:54440/1‑149 (MQ=255)                                       ggtattattCTTAAAAATACTTTTTAAGAATTTTTATATACCCATTTAAGAGAACACGCCTGTGGTAGGAGGGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATggg                                                                                                               <  6:130870/146‑1 (MQ=38)                                             attCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTc                                                                                                          <  3:193934/148‑1 (MQ=255)                                              ttCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCa                                                                                                         <  2:54440/148‑1 (MQ=255)                                              ttCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAg                                                                                                        <  6:253882/149‑1 (MQ=255)                                              ttCTAAAAAATCATGTTAACCAATTTTTATAAAACCATTTAAGAGACCCGGGCCCTGGTAGGAGGGGGTTTACCCGCTAATACGGGGGCCCTGGCCCCATCGCTTTCGGGGTTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCa                                                                                                         <  4:289459/148‑1 (MQ=255)                                                   aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATgg                                                                                                                >  8:187209/1‑136 (MQ=255)                                                   aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATgg                                                                                                                <  7:187209/136‑1 (MQ=255)                                                      aaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGttt                                                                                                <  5:104628/149‑1 (MQ=255)                                                                aaGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcg                                                                                      >  4:77776/1‑149 (MQ=255)                                                                              aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCA‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                                          >  7:177834/1‑147 (MQ=255)                                                                               aCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                                       >  5:90868/1‑149 (MQ=255)                                                                                   tttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCat                                                                                                                    >  7:72817/1‑100 (MQ=255)                                                                                   tttAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCat                                                                                                                    <  8:72817/100‑1 (MQ=255)                                                                                        gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGCTGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTGCCAGGCCACGTAc                                                              >  5:301465/1‑149 (MQ=255)                                                                                           aaCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc                                                           >  5:240733/1‑149 (MQ=255)                                                                                           aaCACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCACGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCg                                                                                           >  6:228228/1‑117 (MQ=255)                                                                                                       gTAGGAGCGGTTTAACCCGCGAAAACGGTGGACGATGCACCATCGCATTCGCGGTTAAACCCGCTCC‑TACCAGGGGCCATATGGCCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                                                <  3:96537/148‑1 (MQ=255)                                                                                                             gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                         <  6:240733/149‑1 (MQ=255)                                                                                                                ggTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCCGGAACAGCAGCGTc                                      >  8:32412/1‑149 (MQ=255)                                                                                                                 gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGt                                       >  7:206976/1‑147 (MQ=255)                                                                                                                      cccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCaa                                 <  5:222156/148‑1 (MQ=255)                                                                                                                       ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg                               >  4:174087/1‑149 (MQ=255)                                                                                                                            aaTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGcc                          <  6:90868/149‑1 (MQ=255)                                                                                                                                     gtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                         <  6:232447/125‑1 (MQ=255)                                                                                                                                     gtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                                         >  5:232447/1‑125 (MQ=255)                                                                                                                                       cGATGCACCATCGCATTCGCGGGTAAACCCGCTCCGTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                                   >  1:207209/1‑114 (MQ=255)                                                                                                                                       cGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                                   <  2:207209/113‑1 (MQ=255)                                                                                                                                        gATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCACGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc               <  8:177834/148‑1 (MQ=255)                                                                                                                                         aTGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCa                         <  7:278617/137‑1 (MQ=255)                                                                                                                                               cATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgat                     <  7:50085/135‑1 (MQ=255)                                                                                                                                               cATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATAAGGTCCTTCATGCCTCTTTCTGCCGAACGCAGATCTAGAACCAGTCCACCTACAGCGCCGGCAGGAACACCAGCGTCACCAAGCTCGCCCCgaa                     >  8:50085/1‑134 (MQ=255)                                                                                                                                                   gCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCAATCATGGc   <  8:206976/149‑1 (MQ=255)                                                                                                                                                    cATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCAATCATGGCg  >  3:20069/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCCGCCAATCATGGCG  >  NC_002947/5897644‑5897939

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