Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,051,670	T→G	11.0%	intergenic (+19/‑38)	ridA → / → PP_5304	aminoacrylate/iminopropionate hydrolase/deaminase/lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,051,670	0	T	G	11.0%	115.9 / 5.9	45	intergenic (+19/‑38)	ridA/PP_5304	aminoacrylate/iminopropionate hydrolase/deaminase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (22/17);  new base G (2/3);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCAAGGTAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCTGGATCAACCAGGAAGCCATCA  >  NC_002947/6051525‑6051813                                                                                                                                                    |                                                                                                                                                cGCCAAGGTAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGc                                                                                                                                               <  1:303635/148‑1 (MQ=255)   ccAAGGTAAACGAGGTGATGGGCCGCTACTTCGAACTGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCg                                                                                                                                            >  1:175325/1‑149 (MQ=255)     aaGGTAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCGGGTGACGGGCAGGGCCCGTc                                                                                                                                          <  4:295014/149‑1 (MQ=255)        gTAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCAcc                                                                                                                                       >  7:199245/1‑149 (MQ=255)         tAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCGACATGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCAc                                                                                                                                        >  1:282087/1‑149 (MQ=255)         tAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCAccc                                                                                                                                      <  8:140762/149‑1 (MQ=255)         gaaaCGAGGTGATGGGCCGCTACGTCGACCAGCCCTACCCGGCCCGCGCCGCCCTCGGCG‑‑TGGCCGCACTGCCGACCGGCGCCCAGGGCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCAGGGCCCGTCAcca                                                                                                                                      <  1:327010/148‑2 (MQ=255)                ggTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCGGCCCGTAAcccaacatc                                                                                                                                >  2:350025/1‑142 (MQ=255)                   gATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGCCGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCAGGGCCCGTCAccag                                                                                                                                     <  5:348424/140‑3 (MQ=255)                   gATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCAGGGCCCGTCAccag                                                                                                                                     >  6:348424/1‑138 (MQ=255)                         ccGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCg                                                                                                                                            <  6:223156/127‑1 (MQ=255)                         ccGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCg                                                                                                                                            >  5:223156/1‑127 (MQ=255)                               cTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTAccc                                                                                                                >  3:228015/1‑149 (MQ=255)                                         cccTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCATGGTTACCCCGTCATGCGt                                                                                                      <  2:282087/149‑1 (MQ=255)                                         cccTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCAGCTGCAAGGTTACCCCGTCATGct                                                                                                       >  6:25853/1‑147 (MQ=255)                                            tACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTc                                                                                                     >  5:24295/1‑147 (MQ=255)                                               ccGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAg                                                                                                  >  4:90937/1‑147 (MQ=255)                                               ccGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAg                                                                                                  <  3:90937/147‑1 (MQ=255)                                                cGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCAc                                                                                               >  2:95050/1‑149 (MQ=255)                                                  gCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTg                                                                                             <  6:24295/149‑1 (MQ=255)                                                          cgccATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCCCCCTCATTCCTCTTCAAGGTTACCACGGCATGCGTCCAGCAATGCCCTTcgc                                                                                     >  6:259589/1‑149 (MQ=255)                                                              aTCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCAc                                                                                               >  2:123469/1‑135 (MQ=255)                                                              aTCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTg                                                                                 <  1:350025/149‑1 (MQ=255)                                                               tCGGCG‑‑TTGCCGCCCTGCCGAAAGGCGCCCAGGTCGAAATGGCCGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTgg                                                                                >  7:145174/1‑149 (MQ=255)                                                                          ccGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGc                                                                       <  8:199245/149‑1 (MQ=255)                                                                           tgCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGc                                                                       >  2:331642/2‑148 (MQ=255)                                                                                  ccGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGggc                                                               <  8:109846/149‑1 (MQ=255)                                                                                  ccGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTCCTACTggg                                                                <  2:175325/148‑1 (MQ=255)                                                                                   cGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGggcg                                                              >  1:242646/1‑149 (MQ=255)                                                                                    gAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTg                                                                                             >  4:319152/1‑117 (MQ=255)                                                                                    gAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTg                                                                                             <  3:319152/117‑1 (MQ=255)                                                                                         gcgcCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCa                                                        <  1:331642/149‑1 (MQ=255)                                                                                              cAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGc                                                   >  8:86635/1‑149 (MQ=255)                                                                                                      aaTGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAAcc                                           <  5:354015/149‑1 (MQ=255)                                                                                                          gACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGaa                                       <  2:242646/149‑1 (MQ=255)                                                                                                           aCGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCg                                         >  1:83851/1‑146 (MQ=255)                                                                                                            cGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCg                                         >  7:273318/1‑145 (MQ=255)                                                                                                                          tCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCt                       >  2:272456/1‑149 (MQ=255)                                                                                                                                 atgCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGGCAAGCACTGCCCTTAGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCTGGATCaa                >  1:250844/1‑149 (MQ=255)                                                                                                                                    cccTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTTCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCTGGATCAAcc              <  5:25853/148‑1 (MQ=255)                                                                                                                                      cTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTg                                                                  >  7:17214/1‑94 (MQ=255)                                                                                                                                      cTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTg                                                                  <  8:17214/94‑1 (MQ=255)                                                                                                                                           gACGGGCCCTGGCCGTCAACCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCTGGATCAACCAGGAAGcc      <  8:145174/149‑1 (MQ=255)                                                                                                                                             cGGGCCCTGCCCGTCACCCTCCTTCCCCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCTGGATCAACCAGGAAGCCa     >  8:222475/1‑148 (MQ=255)                                                                                                                                               ggCCCCGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCGCAAACCCGAAGATTTCAACGGCATCTGGATCAACCAGGAAGCCATCa  >  1:175335/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                CGCCAAGGTAAACGAGGTGATGGGCCGCTACTTCGAACAGCCCTACCCGGCCCGCGCCGCCATCGGCG‑‑TTGCCGCACTGCCGAAAGGCGCCCAGGTCGAAATGGACGCCATCCTGGTCATCGAATGATGCCCTGGTGACGGGCCCTGCCCGTCACCCTCCTTCCTCTGCAAGGTTACCCCGTCATGCGTCAAGCACTGCCCTTCGCGCTGGCAGCCCTGCTACTGGGCGGCTGCACCAGCCACAAACCCGAAGATTTCAACGGCATCTGGATCAACCAGGAAGCCATCA  >  NC_002947/6051525‑6051813

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