Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,226,710	T→A	8.7%	intergenic (‑220/‑111)	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,710	0	T	A	8.7%	135.2 / 3.1	46	intergenic (‑220/‑111)	gltI/PP_1072	glutamate/aspartate ABC transporter substrate‑binding protein/leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base A (2/2);  total (21/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

GGGTGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTC  >  NC_002947/1226564‑1226849                                                                                                                                                   |                                                                                                                                            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                                                                                                                        <  6:163916/149‑1 (MQ=255)                    gtgtAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAc                                                                                                                        <  2:316967/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                                                                            <  5:72656/142‑1 (MQ=255)                                                                       aTTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAgg                                                                            >  6:72656/1‑142 (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                                                                                                                         >  6:107276/1‑71 (MQ=255)                                                                                                 tctctcCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCaa                                                                                                                         <  5:107276/71‑1 (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)                                                                                                                                                  aTCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcaca                                <  7:249548/111‑1 (MQ=255)                                                                                                                                                  aTCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcaca                                >  8:249548/1‑111 (MQ=255)                                                                                                                                                   |                                                                                                                                            GGGTGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTCGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTC  >  NC_002947/1226564‑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: ‑