Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,897,779	G→C	12.4%	intergenic (‑179/+55)	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,779	0	G	C	12.4%	120.7 / 5.2	47	intergenic (‑179/+55)	PP_5172/PP_5173	hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (21/20);  new base C (2/4);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC  >  NC_002947/5897635‑5897926                                                                                                                                                 |                                                                                                                                                     ggTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGc                                                                                                                                                  <  1:32018/149‑1 (MQ=255)   ttggAAGGGGTTTTGCGCCTGTCTTTCGGGCAATAGTAGTGAGCCGGGATTATTCTTAAAAATACTGTTTAATAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                                <  8:103616/146‑1 (MQ=255)      gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCAGCGAATACGGTGGTCGATGAACCAt                                                                                                                                             >  7:122272/1‑149 (MQ=255)       aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                            <  6:14355/149‑1 (MQ=255)         ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCg                                                                                                                                           >  6:104628/1‑148 (MQ=255)          gggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAACCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                         >  5:130870/1‑149 (MQ=255)              tttGCGCCTGGCTTGAGGGCAATATTAGTGAGGAGGTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTc                                                                                                                                      <  8:155423/148‑1 (MQ=255)                   gcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                                 >  4:193934/1‑148 (MQ=255)                     cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                    >  7:41822/1‑143 (MQ=25)                     cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                    <  8:41822/143‑1 (MQ=25)                        gCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAAc                                                                                                                            <  1:19959/148‑1 (MQ=255)                           tGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAAGACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCAGGGTATTGGCGGGTAAACCCg                                                                                                                         <  4:150442/148‑1 (MQ=25)                            gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGt                                                                                                                                <  2:165627/140‑1 (MQ=255)                            gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGt                                                                                                                                >  1:165627/1‑140 (MQ=255)                               ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                    <  2:24742/133‑1 (MQ=25)                               ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgc                                                                                                                                    >  1:24742/1‑133 (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)                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGa                                                                                                                                                     >  5:276961/1‑89 (MQ=255)                                                          ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGa                                                                                                                                                     <  6:276961/89‑1 (MQ=255)                                                            aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATgg                                                                                                   <  7:187209/136‑1 (MQ=255)                                                            aaaaaTACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATgg                                                                                                   >  8:187209/1‑136 (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)                                                                                                                                                 |                                                                                                                                                     GGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCC‑TACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC  >  NC_002947/5897635‑5897926

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: ‑