Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,203,832	G→C	8.2%	intergenic (+96/‑364)	PP_4583 → / → yajD	putative Peptidase/endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,203,832	0	G	C	8.2%	149.3 / 4.8	49	intergenic (+96/‑364)	PP_4583/yajD	putative Peptidase/endonuclease
Reads supporting (aligned to +/- strand):  ref base G (30/15);  new base C (2/2);  total (32/17)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

GATGTGGCTCTGGAGTATGCGTTCGTGTTTGGCGAGCTGGGTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCACCAGGCACTGGC  >  NC_002947/5203686‑5203980                                                                                                                                                   |                                                                                                                                                         gATGTGGCTCTGGAGTATGCGTTCGTGTTTGGCGAGCTGGGTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGc                                                                                                                                                         >  3:352361/1‑148 (MQ=255)          ctGGAGTATGCGTTCGTGTTTGGCGAGCTGGGTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGAt                                                                                                                                               >  3:181365/1‑149 (MQ=255)                         gtgtTTGGCGAGCTGGGTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGcgggccgg                                                                                                                                >  8:117026/1‑149 (MQ=255)                                        ggTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATc‑ggcg                                                                                                                >  6:83690/1‑149 (MQ=255)                                         gTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATc‑ggcgg                                                                                                               <  7:117026/149‑1 (MQ=255)                                               tttGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑g                                                                                                        >  6:186646/1‑149 (MQ=255)                                                 tGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGGAATGGCGGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTg                                                                                                      >  3:67537/1‑149 (MQ=255)                                                             gggCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCggct                                                                                        >  6:126120/1‑149 (MQ=255)                                                               gcgcGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCggctgg                                                                                      >  6:46091/1‑149 (MQ=255)                                                               gcgcGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCggctgg                                                                                      <  5:83690/149‑1 (MQ=255)                                                               gcgcGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCggctgg                                                                                      <  4:352361/149‑1 (MQ=255)                                                                cgcgTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCggctggc                                                                                     >  5:161490/1‑149 (MQ=255)                                                                            cccAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAg                                                                            <  4:185703/146‑1 (MQ=255)                                                                            cccAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAg                                                                            >  3:185703/1‑146 (MQ=255)                                                                                 cccccTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGc                                                                    >  1:200296/1‑149 (MQ=255)                                                                                   cccTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCgg                                                                  >  2:229218/1‑149 (MQ=255)                                                                                    ccTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCCGGCACTGCCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCAGGCTGTCGTAAGCCTTGGGCGgt                                                                 >  7:242725/1‑149 (MQ=255)                                                                                           tCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCg                                                          >  4:108278/1‑149 (MQ=255)                                                                                             gCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTaa                                                                             <  2:321932/128‑1 (MQ=255)                                                                                               aCCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTg                                                                                                      <  3:222624/103‑1 (MQ=255)                                                                                               aCCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTg                                                                                                      >  4:222624/1‑103 (MQ=255)                                                                                               aCCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGt                                                      >  4:226023/1‑149 (MQ=255)                                                                                                 ccGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGtt                                                     >  3:180376/1‑148 (MQ=255)                                                                                                   gTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTAt                                                   <  8:5907/148‑1 (MQ=255)                                                                                                   gTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTAt                                                   >  7:5907/1‑148 (MQ=255)                                                                                                   gTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATc                                                  >  2:25109/1‑149 (MQ=255)                                                                                                   gTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATc                                                  >  8:174361/1‑149 (MQ=255)                                                                                                        tcatcCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCgg                                             <  5:186646/149‑1 (MQ=255)                                                                                                             ctccACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTAt                                         <  6:161490/148‑1 (MQ=255)                                                                                                             ctccACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATg                                        >  6:104942/1‑149 (MQ=255)                                                                                                              tccACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATgg                                       >  5:167634/1‑149 (MQ=255)                                                                                                                cacaACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATCAAGAGGCTACACATGAAGGCGCCCCGGGGGGGTGTGGTAATGCATGGGCGGTGTACTCGGGGTTATGAGCGGGttt                                         >  1:112151/1‑147 (MQ=255)                                                                                                                 acaACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAATA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGc                                                                    <  2:9594/117‑1 (MQ=255)                                                                                                                 acaACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAATA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGc                                                                    >  1:9594/1‑117 (MQ=255)                                                                                                                  caACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTg                                   <  2:200296/149‑1 (MQ=255)                                                                                                                       ggCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGggcg                               >  6:51206/1‑148 (MQ=255)                                                                                                                         cTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGggcgtag                            <  4:181365/149‑1 (MQ=255)                                                                                                                          ttCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGggcgtagg                           >  7:251266/1‑149 (MQ=255)                                                                                                                             ggcggcTCCGGCCGCAACCCGGGCAATGGCAGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGggcgtaggcg                         >  1:28359/1‑148 (MQ=255)                                                                                                                                   tCCGGCCGCAACCCGGCCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTgg                                  >  3:166546/1‑133 (MQ=255)                                                                                                                                   tCCGGCCGCAACCCGGCCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTgg                                  >  3:161076/1‑133 (MQ=255)                                                                                                                                   tCCGGCCGCAACCCGGCCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTgg                                  <  4:161076/133‑1 (MQ=255)                                                                                                                                   tCCGGCCGCAACCCGGCCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTgg                                  <  4:166546/133‑1 (MQ=255)                                                                                                                                     cGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGt                 <  3:108278/148‑1 (MQ=255)                                                                                                                                         cGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCgg                                                         >  7:5025/1‑104 (MQ=255)                                                                                                                                         cGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCgg                                                         <  8:5025/104‑1 (MQ=255)                                                                                                                                         cGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCAc             >  1:297331/1‑148 (MQ=255)                                                                                                                                            aaCCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCACCAgg         <  2:297331/149‑1 (MQ=255)                                                                                                                                                cGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCACAAGGCACt     >  5:264273/1‑149 (MQ=255)                                                                                                                                                   gCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCACCAGGCACTGGc  <  7:174361/149‑1 (MQ=255)                                                                                                                                                   gCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCACCAGGCACTGGc  <  2:28359/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                         GATGTGGCTCTGGAGTATGCGTTCGTGTTTGGCGAGCTGGGTATTCTTTGAGATAGCGGGGGGCGCGTTGCACCCCCCAGCCCCCTCAGCTCGCACCCGTCAATCATCCTCCACAACCGGCTTCGGCGGCTCCGGCCGCAACCCGGGCAATGGCTGATCCTTCGGCGGGCCGGGCACTGGCATC‑GGCGGCAACA‑GTG‑GCGC‑CCCCGGCTGGCTGTCGTAAGCCTTGGGCGGTGTACTCGGGGTTATCTGCGGGTATGGCGTGGGCGTAGGCGTACCCGGTGCACCAGGCACTGGC  >  NC_002947/5203686‑5203980

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