Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,226,707	T→A	8.9%	intergenic (‑217/‑114)	gltI ← / → PP_1072	glutamate/aspartate ABC transporter substrate‑binding protein/leucine‑rich repeat‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,226,707	0	T	A	8.9%	127.2 / 3.2	45	intergenic (‑217/‑114)	gltI/PP_1072	glutamate/aspartate ABC transporter substrate‑binding protein/leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base A (2/2);  total (20/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 = 8.97e-01

AGCGGGTGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTC  >  NC_002947/1226561‑1226849                                                                                                                                                   |                                                                                                                                               agCGGGTGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTa                                                                                                                                              <  7:149156/149‑1 (MQ=255)    gggTGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATcc                                                                                                                                           <  2:305346/149‑1 (MQ=255)          aTTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAAc                                                                                                                                       >  3:129553/1‑147 (MQ=255)          aTTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAACga                                                                                                                                     <  4:129553/149‑2 (MQ=255)             cGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCAt                                                                                                                                   <  5:177146/148‑1 (MQ=255)                    cAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGc                                                                                                                            >  7:16433/1‑148 (MQ=255)                     aaGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGcc                                                                                                                           >  2:81962/1‑148 (MQ=255)                       gtgtAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAc                                                                                                                        <  1:81962/149‑1 (MQ=255)                       gtgtAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAc                                                                                                                        <  2:316967/149‑1 (MQ=255)                       gtgtAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAc                                                                                                                        <  6:163916/149‑1 (MQ=255)                            gagCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATc                                                                                                                   >  3:273416/1‑149 (MQ=255)                              gCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATcc                                                                                                                  <  8:16433/148‑1 (MQ=255)                                  gACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAACga                                                                                                                                     >  1:16295/1‑124 (MQ=255)                                  gACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAACga                                                                                                                                     <  2:16295/125‑2 (MQ=255)                                           cAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCg                                                                                                    >  8:210747/1‑149 (MQ=255)                                                         gACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGt                                                                                       <  7:210747/148‑1 (MQ=255)                                                            cTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAc                                                                                                                        >  3:127051/1‑112 (MQ=255)                                                            cTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAc                                                                                                                        <  4:127051/112‑1 (MQ=255)                                                            cTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTgcgc                                                                                   >  1:255169/1‑149 (MQ=255)                                                            cTTAAACCACTCGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGACAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGAATGGGCTGCCAACCTATCCAGGTAAAGCTTGCGGTAACTATGCAGTacg                                                                                    >  8:72333/1‑148 (MQ=255)                                                             ttAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGcgcg                                                                                  >  5:312497/1‑149 (MQ=255)                                                                 aCCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGc                                                                                >  4:98524/1‑147 (MQ=255)                                                                      tGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCAt                                                                                                                                   >  4:248931/1‑91 (MQ=255)                                                                      tGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCAt                                                                                                                                   <  3:248931/91‑1 (MQ=255)                                                                          aTTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAgg                                                                            >  6:72656/1‑142 (MQ=255)                                                                          aTTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAgg                                                                            <  5:72656/142‑1 (MQ=255)                                                                            tGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGAc                                                                    <  4:273416/148‑1 (MQ=255)                                                                              cTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACg                                                                   >  7:237161/1‑147 (MQ=255)                                                                              cTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGt                                                                  <  6:29792/148‑1 (MQ=255)                                                                               tCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAAGTATGCAGTGCGCGGCAAGGGATCAGACGAAg                                                                >  6:252679/1‑149 (MQ=255)                                                                                      tGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGa                                                                     <  4:75622/137‑1 (MQ=255)                                                                                      tGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGa                                                                     >  3:75622/1‑137 (MQ=255)                                                                                        aCAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTg                                                       <  6:312497/149‑1 (MQ=255)                                                                                          aGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGct                                                     <  8:237161/149‑1 (MQ=255)                                                                                                    tctctcCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCaa                                                                                                                         <  5:107276/71‑1 (MQ=255)                                                                                                    tctctcCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCaa                                                                                                                         >  6:107276/1‑71 (MQ=255)                                                                                                     cgcatcagggTTCACCAGGCCGAATAGACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGga                                          <  3:168125/142‑1 (MQ=255)                                                                                                             ggTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCca                                  <  7:72333/149‑1 (MQ=255)                                                                                                                           aaTAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcac                                 <  2:242615/136‑1 (MQ=255)                                                                                                                           aaTAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcac                                 >  1:242615/1‑136 (MQ=255)                                                                                                                           aaTAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCAAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcac                                 <  6:291613/136‑1 (MQ=255)                                                                                                                           aaTAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCAAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcac                                 >  5:291613/1‑136 (MQ=255)                                                                                                                                   gCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTcc            <  3:98524/149‑1 (MQ=255)                                                                                                                                          cTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCttt                                                 <  8:243520/105‑1 (MQ=255)                                                                                                                                              ttCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTc  >  7:81682/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                               AGCGGGTGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTC  >  NC_002947/1226561‑1226849

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑