Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,595,670	T→A	10.1%	intergenic (‑248/+6)	kgtP ← / ← dctD‑III	alpha‑ketoglutarate permease/C4‑dicarboxylate transport transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,595,670	0	T	A	10.1%	137.9 / 4.5	49	intergenic (‑248/+6)	kgtP/dctD‑III	alpha‑ketoglutarate permease/C4‑dicarboxylate transport transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base A (3/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCATAGCGGTGGGTGTTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGCCCCCATGCTGACGCAAGGCCGCACGCAGGCA  >  NC_002947/1595543‑1595813                                                                                                                                |                                                                                                                                                 cGGCATAGCGGTGGGTGTTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATa                                                                                                                            <  2:256931/149‑1 (MQ=255)    cATAGCGGTGGGTGTTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATaaaa                                                                                                                         <  7:313493/149‑1 (MQ=255)    cATAGCGGGGTTTGTTGCGACGCACTTTTCATTATCTGTCCTGACCGCTAAGGGCTTCAAAGCAGTCCAAAACCTTGAGCTTATTAGATGTTTTTAGCCTATCAGTTTCAGTTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATaaaa                                                                                                                         <  8:293897/149‑1 (MQ=255)           gtgggtgTTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTcgc                                                                                                                  <  8:201243/149‑1 (MQ=255)           gtgggtgTTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTcgc                                                                                                                  >  7:182459/1‑149 (MQ=255)                 gtTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCAtg‑cctgc                                                                                                                                  >  6:231250/1‑127 (MQ=255)                 gtTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCAtg‑cctgc                                                                                                                                  <  5:231250/127‑1 (MQ=255)                        cgccCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATg                                                                                                     <  8:182459/149‑1 (MQ=255)                        cgccCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATg                                                                                                     <  6:79372/149‑1 (MQ=255)                            cTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGc                                                                                                          <  1:156535/140‑1 (MQ=255)                            cTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGc                                                                                                          >  2:156535/1‑140 (MQ=255)                             tttGCAGGAGCGGTCCTGACCTCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGtt                                                                                                <  5:229806/149‑1 (MQ=255)                                  aGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCAtc                                                                                           >  5:340176/1‑149 (MQ=255)                                  aGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCAt                                                                                            <  6:25413/148‑1 (MQ=255)                                  aGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCAt                                                                                            >  5:25413/1‑148 (MQ=255)                                        ggTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGccc                                                                                                        >  2:161392/1‑130 (MQ=255)                                        ggTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGccc                                                                                                        <  1:161392/130‑1 (MQ=255)                                          tCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGt                                                                                    >  7:246247/1‑148 (MQ=255)                                          tCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGG‑CATGAGCGGAAATCCGCTTATgacctgc                                                                                                                                  <  1:263365/102‑1 (MQ=37)                                          tCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGG‑CATGAGCGGAAATCCGCTTATgacctgc                                                                                                                                  >  2:263365/1‑102 (MQ=37)                                            cTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTaa                                                                                 <  8:246247/149‑1 (MQ=255)                                                          cTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGtt                                                                                                <  6:134980/120‑1 (MQ=255)                                                          cTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGtt                                                                                                >  5:134980/1‑120 (MQ=255)                                                           tGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGccc                                                                                                                                <  8:144203/87‑1 (MQ=255)                                                            gCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTc                                                                                                                    >  7:189040/1‑98 (MQ=255)                                                            gCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTc                                                                                                                    <  8:189040/98‑1 (MQ=255)                                                            gCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAgg                                                                    >  2:226647/1‑146 (MQ=255)                                                            gCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCa                                                                  <  1:226647/148‑1 (MQ=255)                                                             cAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAg                                                                 <  4:144289/148‑1 (MQ=255)                                                              aaaGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGa                                                                      >  1:349344/1‑142 (MQ=255)                                                              aaaGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGa                                                                      <  2:349344/142‑1 (MQ=255)                                                                   aGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGc                                                             <  4:95027/146‑1 (MQ=255)                                                                   aGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGc                                                             >  3:95027/1‑146 (MQ=255)                                                                       cAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCt                                                      >  6:11404/1‑149 (MQ=255)                                                                        aaaaCCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGc                                                                               >  8:20122/1‑123 (MQ=255)                                                                        aaaaCCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGc                                                                               <  7:20122/123‑1 (MQ=255)                                                                         aaaCCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTcc                                                    <  6:340176/149‑1 (MQ=255)                                                                                  gCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGAc                                           <  5:11404/149‑1 (MQ=255)                                                                                          gCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGAtctc                                   <  2:126388/149‑1 (MQ=255)                                                                                          gCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGAtctc                                   >  1:126388/1‑149 (MQ=255)                                                                                           cGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCGGTTGCAAGGCCTCCATCACCGACTTGATCTCg                                  >  2:63240/1‑149 (MQ=255)                                                                                           cGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCg                                  >  2:196724/1‑149 (MQ=255)                                                                                                        aTCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGCCCCCATGCTGAc                     >  2:159932/1‑149 (MQ=255)                                                                                                           aGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGCCCCCATGCTGACGc                   >  5:192160/1‑148 (MQ=255)                                                                                                              tGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGAGTTGATCTCGCCCCCATGCTGACGCAAgg               <  1:196724/149‑1 (MQ=255)                                                                                                                 caggcATGAGCGGAAATCCGCTTATGACCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTGGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGc                                 <  1:266941/125‑1 (MQ=38)                                                                                                                 caggcATGAGCGGAAATCCGCTTATGACCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGc                                 >  2:266941/5‑129 (MQ=38)                                                                                                                 caggcATGAGCGGAAATCCGCTTATGACCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGCCCCCATGCTGACGCAAGGcc             >  8:255186/5‑149 (MQ=255)                                                                                                                           cGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGCCCCCATGCTGACGCAAGGCCGCACGCAGGCa  >  4:90781/1‑149 (MQ=255)                                                                                                                                |                                                                                                                                                 CGGCATAGCGGTGGGTGTTGCGCCGCCCTTTGCAGGAGCGGTCCTGACCGCGAAGGGCTGCAAAGCAGTCCAAAACCTTGAGCTTATGAGCGGTTTTTTGCCTATCAGGTGCAGGTCATAAGCGGATTTCCGCTCATG‑CCTGCCCAATAAAATCCTCGCGCACCAGCCCATGCCGTTGCATCTTCTCGTTAAGCGTGCGCCGAGGCAGTTGCAAGGCCTCCATCACCGACTTGATCTCGCCCCCATGCTGACGCAAGGCCGCACGCAGGCA  >  NC_002947/1595543‑1595813

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