Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	109,959	A→C	9.7%	intergenic (+22/‑22)	PP_0103 → / → ctaD	cytochrome c oxidase subunit 2/cytochrome c oxidase subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	109,959	0	A	C	9.7%	122.5 / 5.2	48	intergenic (+22/‑22)	PP_0103/ctaD	cytochrome c oxidase subunit 2/cytochrome c oxidase subunit 1
Reads supporting (aligned to +/- strand):  ref base A (17/26);  new base C (3/2);  total (20/28)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAAGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGTT  >  NC_002947/109811‑110107                                                                                                                                                     |                                                                                                                                                     cAAGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaa                                                                                                                                                     >  1:603601/1‑150 (MQ=255)  aaGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaat                                                                                                                                                    <  2:249493/150‑1 (MQ=255)   agcagcTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATg                                                                                                                                                   <  1:109322/150‑1 (MQ=255)          gTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGCAGCCGTCACAAGCTAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCa                                                                                                                                            <  2:710777/150‑1 (MQ=255)                  tCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAg                                                                                                                                    <  2:323507/150‑1 (MQ=255)                    gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGa                                                                                                                                  <  2:540705/150‑1 (MQ=255)                    gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGa                                                                                                                                  <  2:821713/150‑1 (MQ=255)                    gACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGACAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCCGAAGCCAAGTGAACCGGGTTACTCCCAGAAGACCAGAGATTGCAGGAGACAGGa                                                                                                                                  >  2:257178/1‑150 (MQ=255)                            tgctgTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTg                                                                                                                          >  1:559381/1‑150 (MQ=255)                              ctgTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCa                                                                                                                        <  1:878179/150‑1 (MQ=255)                                gTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGt                                                                                                                      <  2:1070975/150‑1 (MQ=255)                                gTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGt                                                                                                                      <  1:685207/150‑1 (MQ=255)                                 tAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTg                                                                                                                     >  1:541375/1‑150 (MQ=255)                                    tCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATc                                                                                                                  <  2:373287/150‑1 (MQ=255)                                      aCCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAg                                                                                                                                           <  2:694631/123‑1 (MQ=255)                                      aCCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAg                                                                                                                                           >  1:694631/1‑123 (MQ=255)                                       ccTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgac                                                                                                               <  2:559381/150‑1 (MQ=255)                                       ccTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgac                                                                                                               <  1:538028/150‑1 (MQ=255)                                       ccTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgac                                                                                                               <  1:305568/150‑1 (MQ=255)                                        cTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATcgacg                                                                                                              >  2:36534/1‑150 (MQ=255)                                             aGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCAc                                                                                                         <  2:603601/150‑1 (MQ=255)                                                     gCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGACCACgg                                                                                                 <  1:535161/150‑1 (MQ=255)                                                                 aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACa                                                                                                                                     <  1:234976/102‑1 (MQ=255)                                                                 aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACa                                                                                                                                     >  2:234976/1‑102 (MQ=255)                                                                 aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGt                                                                                                                           >  1:788142/1‑112 (MQ=255)                                                                 aacaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGt                                                                                                                           <  2:788142/112‑1 (MQ=255)                                                                  acaaAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCAAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATgcg                                                                                    >  2:942564/1‑150 (MQ=255)                                                                     aaGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCAc                                                                                 <  1:692542/150‑1 (MQ=255)                                                                             aTGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaa                                                                                                                                                     <  2:1109337/74‑1 (MQ=255)                                                                             aTGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGcaatcaa                                                                                                                                                     >  1:1109337/1‑74 (MQ=255)                                                                               ggTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCa                                                                       >  1:727968/1‑150 (MQ=255)                                                                                                tGCTGGCGCAGAAGCAGGCAGAAGACAAGTGAACCGGGTTACTCGCGACAATAAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCATGATGCGCTgg                                                      <  2:541375/150‑1 (MQ=255)                                                                                                 gCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTggg                                                     <  2:727968/150‑1 (MQ=255)                                                                                                         tGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGacc                                             >  1:537199/1‑150 (MQ=255)                                                                                                             gcaggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGAccacca                                         >  2:155781/1‑150 (MQ=255)                                                                                                             gcaggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGAccacca                                         >  2:153871/1‑150 (MQ=255)                                                                                                               aggcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaac                                       <  2:895590/150‑1 (MQ=255)                                                                                                                 gcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaa                                        >  1:1108096/1‑147 (MQ=255)                                                                                                                 gcagAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCaccaa                                        <  2:1108096/147‑1 (MQ=255)                                                                                                                     aaGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATaaa                                 <  2:537199/150‑1 (MQ=255)                                                                                                                               aaCCGGGTTACTCGCGGCAATCAATGGATTGAAGGAGACAGGACATAAGTACAGTGAGCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCa                       >  1:415148/1‑150 (MQ=255)                                                                                                                                   gggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      <  2:988428/99‑1 (MQ=255)                                                                                                                                   gggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCaa                                                                      >  1:988428/1‑99 (MQ=255)                                                                                                                                    ggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGAt                                                             <  2:1129226/107‑1 (MQ=255)                                                                                                                                    ggTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGAt                                                             >  1:1129226/1‑107 (MQ=255)                                                                                                                                          tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtg            >  2:684590/1‑150 (MQ=255)                                                                                                                                          tCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtg            <  1:36534/150‑1 (MQ=255)                                                                                                                                           tcctgCAATCCATTGATTGCCGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGa                                               <  1:850032/110‑1 (MQ=255)                                                                                                                                           tcctgCAATCCATTGATTGCCGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGa                                               >  2:850032/5‑114 (MQ=255)                                                                                                                                            gcgGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCtgtg          >  2:917891/1‑150 (MQ=255)                                                                                                                                             ctgCAATCCATTGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGac                                              >  2:213109/3‑113 (MQ=255)                                                                                                                                             ctgCAATCCATTGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGac                                              <  1:213109/111‑1 (MQ=255)                                                                                                                                                  atcaatGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAg    <  1:942564/150‑1 (MQ=255)                                                                                                                                                    caatGGATTGCAGGAGACAGGTCATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCCGCCACGCGCCTGATGCCCTAGCTGCTGACCACCAACCATAAAGACATCTGCACGATGTACCTGTGGTTCAGtt  <  2:49654/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CAAGCAGCTGTCGGAAGTCGACATCGCTGCTGTAGTCACCTACGAGCGCAATGCCTGGGGCAACAACAAAGGCGACATGGTGACGCCGAAGGACGTGCTGGCGCTGAAGCAGGCAGAAGCCAAGTGAACCGGGTTACTCGCGGCAATCAATGGATTGCAGGAGACAGGACATGAGTGCAGTGATCGACGACCACGCCCACGGCCACGACCATGCGCACGGCCCGGCCAAGGGCCTGATGCGCTGGGTGCTGACCACCAACCATAAAGACATCGGCACGATGTACCTGTGGTTCAGTT  >  NC_002947/109811‑110107

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