Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,804,021	C→G	9.6%	intergenic (‑105/‑39)	PP_3361 ← / → PP_3362	aminotransferase/AsnC family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,804,021	0	C	G	9.6%	148.5 / 3.6	51	intergenic (‑105/‑39)	PP_3361/PP_3362	aminotransferase/AsnC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base G (3/2);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAG  >  NC_002947/3803881‑3804163                                                                                                                                             |                                                                                                                                                cGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGc                                                                                                                                         >  7:287497/1‑148 (MQ=255)      aCCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTCGGAA‑CCGGCATAATg                                                                                                                                   <  2:290785/149‑1 (MQ=255)      aCCGAGGTGCCGGGGGGGAAAGGTCTGTCATATCATCAGGGTCCGAGCGCGCTGATTGGTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTCGGAA‑CCGGCATAATg                                                                                                                                   <  6:87362/149‑1 (MQ=255)        cGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGAc                                                                                                                                 >  1:96953/1‑149 (MQ=255)              gCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAAt                                                                                                                                    <  4:10759/140‑1 (MQ=255)              gCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAAt                                                                                                                                    >  3:10759/1‑140 (MQ=255)               ccGTGTGGGTATGGTCTGTCTTATCATCAGGTTTCTAGCGCGTTGTTTGTTTTTCAATTGCTTGTAATGCTCCGCTGAATGGGACAGGGTGCATACGTTCGGGCGCAAGAGTCATTATGGCGGTTCCGA‑CCGGCATAATGAGGGctat                                                                                                                           <  8:114403/148‑4 (MQ=255)                 gggTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCt                                                                                                                         >  3:181899/1‑148 (MQ=255)                    tGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATc                                                                                                                     <  3:243203/149‑1 (MQ=255)                               tGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGG‑TCGGACCCGGCATAAt                                                                                                                                    >  2:141280/1‑123 (MQ=255)                                 tCATATCCTCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGATCCGGACACCCTGCATACGGTCGGGCGCTACAGGCATTATGCCGGTTCGGA‑CGGGCAGAAGGACGGCAAGCTGATCCCGGGGCTGCatt                                                                                                        >  2:8567/1‑147 (MQ=255)                                     atcatcAGGTTCCTAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCGTATCAAGGGGCTGCAATGCGt                                                                                                    <  8:287497/149‑1 (MQ=255)                                        atcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGAc                                                                                                 >  6:217300/1‑149 (MQ=255)                                         tcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGTGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGAcc                                                                                                >  5:34785/1‑149 (MQ=255)                                         tcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGAcc                                                                                                >  5:93647/1‑149 (MQ=255)                                           aGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCt                                                                                               <  4:152442/148‑1 (MQ=255)                                                      gcgcTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTa                                                                                       >  1:149806/1‑145 (MQ=255)                                                          tGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATcgc                                                                               >  2:87199/1‑149 (MQ=255)                                                            aTTGTTTTTGAATAGGTTGTCATGCGCCGCCGAAACGGACACCATGCATACGGTCGGTCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATcgc                                                                               <  7:347512/147‑1 (MQ=255)                                                            aTTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCga                                                                             <  5:217300/149‑1 (MQ=255)                                                            aTTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTACCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATcgc                                                                               <  2:96953/147‑1 (MQ=255)                                                             ttGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGTGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCga                                                                             <  6:34785/148‑1 (MQ=255)                                                             ttGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCga                                                                             <  2:177331/148‑1 (MQ=255)                                                             ttGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCga                                                                             <  2:183194/148‑1 (MQ=255)                                                              tGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATcgc                                                                               >  2:88459/1‑145 (MQ=255)                                                               gTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGAt                                                                          >  6:252747/1‑149 (MQ=255)                                                                  tttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCt                                                                       >  3:132192/1‑149 (MQ=255)                                                                  tttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCt                                                                       >  7:21653/1‑149 (MQ=255)                                                                    tCAATAGGTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCtt                                                                      <  5:175313/148‑1 (MQ=255)                                                                    tCAATAGCTTGTCATGCGCAGCCGAAACGGACACCCTTCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGTCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCtt                                                                      <  7:333013/148‑1 (MQ=255)                                                                              gTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACttt                                                            <  6:93647/148‑1 (MQ=255)                                                                                 aTGCGCCTCCTAAACGGTCACCCTGCATACGGTCGGGCTCAACATGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAATCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGaa                                                        <  5:212923/149‑1 (MQ=255)                                                                                 aTGCGCCGCCTAAACGGACACCCTGCATACGGTCGGGCGCAACAGTCATTATGCCGGTTCAGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGaa                                                        <  7:80627/149‑1 (MQ=255)                                                                                     gccgccGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAc                                                    >  8:308997/1‑149 (MQ=255)                                                                                       cgccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgc                                                  <  4:181899/149‑1 (MQ=255)                                                                                         ccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgcgc                                                <  4:132192/149‑1 (MQ=255)                                                                                                 gACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTgcg                                        >  5:38931/1‑149 (MQ=255)                                                                                                    acCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTc                                     <  7:308997/149‑1 (MQ=255)                                                                                                      ccTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCaa                                   <  1:88459/149‑1 (MQ=255)                                                                                                           aTACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGt                              <  5:252747/149‑1 (MQ=255)                                                                                                            tACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGt                              >  6:97241/1‑148 (MQ=255)                                                                                                            tACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGt                              >  8:181512/1‑148 (MQ=255)                                                                                                                 tCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccg                        >  5:173705/1‑149 (MQ=255)                                                                                                                        cAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCAt                 <  1:87199/149‑1 (MQ=255)                                                                                                                         aaCAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCa                  >  1:200100/1‑147 (MQ=255)                                                                                                                            aGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGc             >  4:142922/1‑149 (MQ=255)                                                                                                                              tcATTATGCCGGTCGGAA‑CCGGCATAATGACGGCAAGCTTATCAAGTGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCt           >  7:80652/2‑149 (MQ=255)                                                                                                                                    tGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTccc      >  8:113465/1‑148 (MQ=255)                                                                                                                                     gCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCg     >  3:253737/1‑148 (MQ=255)                                                                                                                                      ccGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGaa   <  7:181512/149‑1 (MQ=255)                                                                                                                                       cGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAg  >  8:49776/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                CGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCGA‑CCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAAG  >  NC_002947/3803881‑3804163

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