Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,376,482	G→C	10.8%	intergenic (+14/‑134)	PP_2975 → / → PP_2976	insertion element ISR1 protein A3/transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,376,482	0	G	C	10.8%	125.7 / 5.9	44	intergenic (+14/‑134)	PP_2975/PP_2976	insertion element ISR1 protein A3/transposase
Reads supporting (aligned to +/- strand):  ref base G (18/21);  new base C (2/3);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGGCGGCATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGAGC  >  NC_002947/3376336‑3376621                                                                                                                                                   |                                                                                                                                             cggcggCATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGgcgc                                                                                                                                           <  1:468155/150‑1 (MQ=255)   gcggcATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑ga                                                                                                                                        <  1:1079104/150‑1 (MQ=255)       cATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATg                                                                                                                                    <  2:1043827/150‑1 (MQ=255)        aTGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATgg                                                                                                                                   <  1:1017849/150‑1 (MQ=255)                       atgaCCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGgttc                                                                                                                    <  1:34294/150‑1 (MQ=255)                       atgaCCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGgttc                                                                                                                    <  1:360757/150‑1 (MQ=255)                           ccAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATg                                                                                                                                    <  1:8511/130‑1 (MQ=255)                           ccAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATg                                                                                                                                    >  2:8511/1‑130 (MQ=255)                            cAGGCTTCGTGAACTGGAAGAGAAGACCCGTCGGCTCGCGAGCGTGTTGGCTTAGTAGCGCCTCAATTCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCa                                                                                                               <  2:1030602/150‑1 (MQ=255)                            cAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATTTCTCAGCAGGCGGTTCGTTCa                                                                                                               <  2:98835/150‑1 (MQ=255)                                            gAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCGAGTATCa                                                                                               >  1:1006944/1‑150 (MQ=255)                                               gaggagAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAAt                                                                                            >  1:928495/1‑150 (MQ=255)                                                aggagAACCGTCGGCTCAAGAAGATGTACGCCGAGGGGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAAtt                                                                                           >  2:781927/1‑150 (MQ=255)                                                  gagaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGc                                                                                                       <  1:1011527/136‑1 (MQ=255)                                                  gagaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGc                                                                                                       >  2:1011527/1‑136 (MQ=255)                                                  gagaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAg                                                                                         >  1:814532/1‑150 (MQ=255)                                                   agaACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGc                                                                                        <  2:123265/150‑1 (MQ=255)                                                     aaCCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGt                                                                                      >  1:355489/1‑150 (MQ=255)                                                       ccGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGc                                                                                    <  1:61784/150‑1 (MQ=255)                                                           cGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGt                                                                                <  2:814532/150‑1 (MQ=255)                                                            ggCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATGgc                                                                                                                                  <  2:1008297/99‑1 (MQ=255)                                                            ggCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTC‑GCGCCGCTGAGATGgc                                                                                                                                  >  1:1008297/1‑99 (MQ=255)                                                                   gaagaTGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGt                                                                        >  1:146715/1‑150 (MQ=255)                                                                            cGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACg                                                               >  2:99365/1‑150 (MQ=255)                                                                            cGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACg                                                               >  1:1071552/1‑150 (MQ=255)                                                                             gCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGgttcgtt                                                                                                                 >  1:811271/1‑99 (MQ=255)                                                                             gCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGgttcgtt                                                                                                                 <  2:811271/99‑1 (MQ=255)                                                                               cGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACt                                                            <  2:147248/150‑1 (MQ=255)                                                                                          tCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAgg                                                                                                                          >  1:841057/1‑77 (MQ=255)                                                                                          tCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAgg                                                                                                                          <  2:841057/77‑1 (MQ=255)                                                                                           cAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCgg                                                <  2:239850/150‑1 (MQ=255)                                                                                             aGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGt                                              <  2:928495/150‑1 (MQ=255)                                                                                             aGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCGCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGt                                              >  2:214719/1‑150 (MQ=255)                                                                                              ggCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGtt                                             >  2:656677/1‑150 (MQ=255)                                                                                                        atcCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGc                                   <  2:684152/150‑1 (MQ=255)                                                                                                        atcCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGc                                   >  1:614623/1‑150 (MQ=255)                                                                                                           cAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTcc                                >  2:392503/1‑150 (MQ=255)                                                                                                                 gCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAAAATGACGGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCtt                          >  1:57516/1‑150 (MQ=255)                                                                                                                  ccATGGCAAAAAAGTGGTCCAGCCTTCTCGGCGCTGC‑GATATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCttt                         <  2:10545/150‑1 (MQ=255)                                                                                                                  ccATCGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGCGTTTCGCGGAGGCTCCACCCttt                         <  2:552216/150‑1 (MQ=255)                                                                                                                              aGTGGTAAAGACATCTCAGCGGAGC‑GAGAAGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCACCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAg             <  2:57516/150‑1 (MQ=255)                                                                                                                                tGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGcc           >  1:1038776/1‑150 (MQ=255)                                                                                                                                      cgCCATCTC‑GCGCCGCTGAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCa                                                                    <  2:255366/86‑1 (MQ=38)                                                                                                                                      cgCCATCTC‑GCGCCGCTGAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCa                                                                    >  1:255366/2‑87 (MQ=38)                                                                                                                                         cATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGagc  <  1:408516/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             CGGCGGCATGGATGCATCATTGATGACCAGGCTGCGTGAACTGGAAGAGGAGAACCGTCGGCTCAAGAAGATGTACGCCGAGGAGCGCCTCAAGGCCGAGATCATCCAGGAGGCCATGGCAAAAAAGTGGTAAAGCCATCTCAGCGGCGC‑GAGATGGCGCAGCAGGCGGTTCGTTCAGGCAAAGCCAGTATCAAATTAGCGTGCCTGTCGTTTGGTATCAGCACGACTGAACCGCCCCGGGTTTCGCGGAGGCTCCACCCTTTGAGAGAATGGAGCCCATATGAGC  >  NC_002947/3376336‑3376621

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