Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,667,908	A→C	12.7%	I304L (ATC→CTC)	ffh →	signal recognition particle GTPase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,667,908	0	A	C	12.7%	141.2 / 14.5	55	I304L (ATC→CTC)	ffh	signal recognition particle GTPase subunit
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base C (3/4);  total (27/28)
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

CGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAC  >  NC_002947/1667761‑1668054                                                                                                                                                    |                                                                                                                                                    cGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGAt                                                                                                                                                    >  4:160106/1‑149 (MQ=255)             gcgcTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGcc                                                                                                                                          >  6:327437/1‑147 (MQ=255)                      gTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCAt                                                                                                                               <  5:327437/149‑1 (MQ=255)                       tACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATc                                                                                                                              <  3:243469/149‑1 (MQ=255)                          gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                                   <  3:118960/125‑1 (MQ=255)                          gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATc                                                                                                                                                   >  4:118960/1‑125 (MQ=255)                          gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAg                                                                                                                                       >  6:85102/1‑137 (MQ=255)                          gTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAg                                                                                                                                       <  5:85102/137‑1 (MQ=255)                           tGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAc                                                                                                                           >  2:144984/1‑148 (MQ=255)                               aTCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggc                                                                                                                      <  7:116026/149‑1 (MQ=255)                                tCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaaggc                                                                                                                      >  1:181153/1‑148 (MQ=255)                                      gCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaag                                                                                                                        <  1:284855/140‑1 (MQ=255)                                      gCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGAcaag                                                                                                                        >  2:284855/1‑140 (MQ=255)                                         taCCGATGAACTTGATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGAACGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGa                                                                                                             >  2:79029/2‑148 (MQ=255)                                               tCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCa                                                                                                                                              >  4:177958/1‑109 (MQ=255)                                               tCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCa                                                                                                                                              <  3:177958/109‑1 (MQ=255)                                                cAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCa                                                                                                                                              >  6:337990/1‑108 (MQ=255)                                                 aaGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTgg                                                                                                     <  2:181153/148‑1 (MQ=255)                                                   gTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGc                                                                                                                                      <  4:82474/113‑1 (MQ=255)                                                       aTCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGaa                                                                                              >  4:230119/1‑149 (MQ=255)                                                              tGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaag                                                                                       <  3:160106/149‑1 (MQ=255)                                                              tGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTg                                                                                          >  5:286566/1‑146 (MQ=255)                                                                ggTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTgaagaa                                                                                     <  1:144984/149‑1 (MQ=255)                                                                    agaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaag‑ggc                                                                                <  2:201302/149‑1 (MQ=255)                                                                     gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGCacg                                                                                                                                   >  3:250016/1‑96 (MQ=38)                                                                     gaagaCCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCGGCCTGCTCGATCAGGCTGAGCacg                                                                                                                                   <  4:250016/98‑3 (MQ=38)                                                                           cGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaag‑ggcaag                                                                             <  6:355184/145‑1 (MQ=255)                                                                           cGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGaag‑ggcaag                                                                             >  5:355184/1‑145 (MQ=255)                                                                                  cTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCt                                                                  >  2:277829/1‑149 (MQ=255)                                                                                   tCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGaaa                                                                                             <  8:156243/122‑1 (MQ=255)                                                                                   tCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGaaa                                                                                             >  7:156243/1‑122 (MQ=255)                                                                                             ccACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTc                                                        <  1:32247/148‑1 (MQ=255)                                                                                             ccACCCCGACCGCCTCGCTTCGCGCATCCGCGGCATGGGCGACGTGCTCCGCCTGATCGAGCAGGACGAGCAGACCATCGACAAGTCCAAGGCCGACAACCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTc                                                        <  2:43855/148‑1 (MQ=255)                                                                                               aCCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCcgc                                                     <  6:286566/149‑1 (MQ=255)                                                                                                ccccGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGc                                                                                                    <  2:332225/102‑1 (MQ=255)                                                                                                   cGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                  >  7:278977/1‑148 (MQ=255)                                                                                                   cGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGAc                                                  <  8:278977/148‑1 (MQ=255)                                                                                                                 cgcgCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAGTGGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCaaa                                     <  1:79029/148‑1 (MQ=255)                                                                                                                  gcgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGTCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATga                                  <  1:277829/149‑1 (MQ=255)                                                                                                                  gcgcATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATg                                   >  2:179769/1‑148 (MQ=255)                                                                                                                          tCGGCGTGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGaa                                                              >  8:267520/1‑113 (MQ=255)                                                                                                                          tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCTAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGAc                                                           >  1:109592/1‑116 (MQ=255)                                                                                                                          tCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGAc                                                           <  2:109592/116‑1 (MQ=255)                                                                                                                                   gCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTc                 >  8:78018/1‑149 (MQ=255)                                                                                                                                       tctGCTCGGCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTcggcg             >  7:106493/3‑149 (MQ=255)                                                                                                                                           cTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACTAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCt         >  6:301176/1‑149 (MQ=255)                                                                                                                                            tcggCCTGCTCGATCAGGCTGAGCAGTCCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCt                  <  2:143332/136‑1 (MQ=255)                                                                                                                                            tcggCCTGCTCGATCAGGCTGAGCAGGCCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCt                  >  1:143332/4‑139 (MQ=255)                                                                                                                                            tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc           >  6:64244/1‑146 (MQ=255)                                                                                                                                            tCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCggcggc           <  5:64244/146‑1 (MQ=255)                                                                                                                                               gCCTGCTCGATCAGGCTGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg     <  8:106493/149‑1 (MQ=255)                                                                                                                                               gCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATg     <  7:78018/149‑1 (MQ=255)                                                                                                                                                ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAg                                               <  2:4756/106‑1 (MQ=255)                                                                                                                                                ccTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAg                                               >  1:4756/1‑106 (MQ=255)                                                                                                                                                   gATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAc  >  6:41810/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CGAGGCGGTGCTGCGCTGTCGGTACGTGCAATCACCGGCAAACCGATCAAGTTCATCGGTATGGGTGAGAAGACCGAGGCCCTCGAGCCGTTCCACCCCGACCGCGTCGCTTCGCGCATCCTCGGCATGGGCGACGTGCTCAGCCTGATCGAGCAGGCCGAGCAGACCATCGACAAGGCCAAGGCCGACAAGCTGGCCAAGAAACTGAAGAAG‑GGCAAGGGCTTCGACCTCGAAGACTTCCGCGACCAGCTGCAGCAAATGAAGAACATGGGCGGCCTCGGCGGCCTGATGGAC  >  NC_002947/1667761‑1668054

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