Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,674,821	A→G	11.1%	Q12R (CAG→CGG)	tnpT‑II →	cointegrate resolution protein T

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,674,821	0	A	G	11.1%	89.3 / 3.4	36	Q12R (CAG→CGG)	tnpT‑II	cointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base A (18/14);  new base G (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

CGATAAGCCGTGATTATTCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCCAGCAGCCGGTGGCGTCCGAAG  >  NC_002947/3674673‑3674964                                                                                                                                                     |                                                                                                                                                 cGATAAGCCGTGATTATTCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTgca                                                                                                                                                  <  3:105703/149‑1 (MQ=255)     aaGCCGTGATTATTCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGcagcag                                                                                                                                              >  1:47134/1‑149 (MQ=255)           tGATTATTCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGTGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGt                                                                                                                                        <  5:251719/149‑1 (MQ=255)                attCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAgg                                                                                                                                    <  1:138617/148‑1 (MQ=255)                 ttCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCa                                                                                                                                  >  3:272059/1‑149 (MQ=255)                            gAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCg                                                                                                                       <  6:280092/149‑1 (MQ=255)                            gAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCg                                                                                                                       >  4:154287/1‑149 (MQ=255)                              aaTGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCgg                                                                                                                      <  4:272059/148‑1 (MQ=255)                              aaTGTTAGCAGTAGGTGGTTGCTATCGTCATAATTTTTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCATGCGCGTCAGGCACTGATCGCCCggg                                                                                                                     <  4:264881/149‑1 (MQ=255)                                    tGCAGTAGGCGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTgag                                                                                                                <  8:313282/148‑1 (MQ=255)                                           ggggTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGTGTGAGAACCCCa                                                                                                         <  2:47134/148‑1 (MQ=255)                                                           aTAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACg                                                                                        >  7:162738/1‑148 (MQ=255)                                                                  gTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCGGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGt                                                                                >  3:322805/1‑149 (MQ=255)                                                                  gTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCGGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGt                                                                                >  3:328100/1‑149 (MQ=255)                                                                   tACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGtt                                                                               >  4:330760/1‑149 (MQ=255)                                                                   tACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCACGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCTCCCGGGGAGAGAACCCCAGCATCGATG‑CCATACGTATCGAGtt                                                                               >  5:251344/1‑149 (MQ=255)                                                                     cGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCGGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTgg                                                                             <  4:328100/149‑1 (MQ=255)                                                                     cGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCGGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTgg                                                                             <  4:322805/149‑1 (MQ=255)                                                                              aTGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCg                                                                    >  7:250615/1‑149 (MQ=255)                                                                                          acaGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCa                                                        <  8:250615/149‑1 (MQ=255)                                                                                              tACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGGGGTGAAAACGTATCGAGTTGGGCAACACCGCCTCCAAAACCACGa                                                     >  1:119963/1‑149 (MQ=255)                                                                                              tACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGtt                                                                               <  6:58917/122‑1 (MQ=255)                                                                                              tACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGtt                                                                               >  5:58917/1‑122 (MQ=255)                                                                                              tACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGAt                                                    <  5:115453/149‑1 (MQ=255)                                                                                                  aaaTCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCAc                                                <  1:23761/149‑1 (MQ=255)                                                                                                        tAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCCTACGTATCGAGTTGGGCAACTCCGGCTCCAAAACCACGATTCACCGCTAc                                          >  5:321741/1‑149 (MQ=255)                                                                                                        tAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTAc                                          >  2:220619/1‑149 (MQ=255)                                                                                                              caagcaTGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAg                                    <  3:330760/149‑1 (MQ=255)                                                                                                              caagcaTGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCACCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCTTCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAg                                    >  5:128233/1‑149 (MQ=255)                                                                                                                   aTGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACt                               >  3:84883/1‑149 (MQ=255)                                                                                                                       cccGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAAc                           <  3:154287/149‑1 (MQ=255)                                                                                                                        ccGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACg                          >  1:231254/1‑149 (MQ=255)                                                                                                                          gTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGa                                                     >  4:50129/1‑120 (MQ=255)                                                                                                                          gTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGa                                                     <  3:50129/120‑1 (MQ=255)                                                                                                                                      acaaGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCAc                                                       >  3:207590/1‑106 (MQ=255)                                                                                                                                           gtggtAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCCAGCAGCCGGTGGCGTc       >  1:212030/1‑149 (MQ=255)                                                                                                                                                aGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCCAGCAGCCGGTGGCGTCCGAAg  >  7:252291/1‑149 (MQ=255)                                                                                                                                                aGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCCAGCAGCCGGTGGCGTCCGAAg  >  7:37079/1‑149 (MQ=255)                                                                                                                                                    cagcagGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGtt                                                                               >  3:194237/1‑68 (MQ=255)                                                                                                                                                    cagcagGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGtt                                                                               <  4:194237/68‑1 (MQ=255)                                                                                                                                                    cagcagGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCcagca                   >  1:223048/1‑128 (MQ=255)                                                                                                                                                    cagcagGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCcagca                   <  2:223048/128‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 CGATAAGCCGTGATTATTCCACCTGAGGAAATGTTTGCAGTAGGGGTTTGCTATCGTCATAATTGGTACGTAATTACATGTTACATAACACAGTACAAAATCGTAGGAGCAAGCATGGCCCGTGGCGGTATCAACAAGGTGGTAGTGCAGCAGGCGCGTCAGGCACTGATCGCCCGGGGTGAGAACCCCAGCATCGATG‑CCATACGTATCGAGTTGGGCAACACCGGCTCCAAAACCACGATTCACCGCTACCTGAAGGAACTGAACGGCCAGCAGCCGGTGGCGTCCGAAG  >  NC_002947/3674673‑3674964

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