Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,512,830	C→T	10.0%	noncoding (139/376 nt)	PP_mr18 →	RNaseP_bact_a ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,512,830	0	C	T	10.0%	104.0 / 2.9	40	noncoding (139/376 nt)	PP_mr18	RNaseP_bact_a ncRNA
Reads supporting (aligned to +/- strand):  ref base C (18/18);  new base T (2/2);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.01e-01

TTGGCGGCGGAGAGTCGATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGG  >  NC_002947/1512686‑1512970                                                                                                                                                 |                                                                                                                                             ttGGCGGCGGAGAGTCGATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCca                                                                                                                                            >  5:162837/1‑147 (MQ=255)         ggAGAGTCGATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCca                                                                                                                                            <  1:146840/139‑1 (MQ=255)         ggAGAGTCGATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCca                                                                                                                                            >  2:146840/1‑139 (MQ=255)            gagTCGATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCg                                                                                                                               >  8:143254/1‑149 (MQ=255)                 gATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCt                                                                                                                            >  4:263748/1‑147 (MQ=255)                 gATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTaa                                                                                                                          <  7:99633/149‑1 (MQ=255)                   tCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGc                                                                                                                        <  6:162837/149‑1 (MQ=255)                           tCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGt                                                                                                                <  8:161136/149‑1 (MQ=255)                           tCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCAGGGCTCCACAGGGCGGATTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGt                                                                                                                >  3:158974/1‑149 (MQ=255)                                 ccTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCGGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCgg                                                                                                           >  7:141775/1‑148 (MQ=255)                                   tGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGt                                                                                                          >  6:126124/1‑147 (MQ=255)                                    gTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAg                                                                                                       >  2:229212/1‑149 (MQ=255)                                     tttttGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTaa                                                                                                        <  4:60625/147‑1 (MQ=255)                                     tttttGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAgg                                                                                                      >  2:189849/1‑149 (MQ=255)                                             ccGCAACTTTTTGTTGTGTTAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGaaaa                                                                                               <  7:216109/148‑1 (MQ=255)                                              cGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAgg                                                                                             <  5:263050/149‑1 (MQ=255)                                              cGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAg                                                                                              <  2:47742/148‑1 (MQ=255)                                               gCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGt                                                                                            >  2:204721/1‑149 (MQ=255)                                                           tgtgtgAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAgc                                                                                 <  5:126124/148‑1 (MQ=255)                                                               tgAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGcgcacc                                                                            >  6:108123/1‑149 (MQ=255)                                                                   cagcagGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGcgcaccgcac                                                                        <  7:143254/149‑1 (MQ=255)                                                                     gcagGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACgg                                                                      <  1:204721/149‑1 (MQ=255)                                                                     gcagGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACgg                                                                      >  7:187847/1‑149 (MQ=255)                                                                         ggggAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAAGTGCGGTAAGAGCGCACCGCACGGCTg                                                                   <  1:229212/148‑1 (MQ=255)                                                                          gggAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTg                                                                                           >  4:65790/1‑123 (MQ=255)                                                                          gggAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTg                                                                                           <  3:65790/123‑1 (MQ=255)                                                                            gaggaAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGTaa                                                               <  8:141775/149‑1 (MQ=255)                                                                                              acaGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACa                                                             >  7:88976/1‑133 (MQ=255)                                                                                              acaGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACa                                                             <  8:88976/133‑1 (MQ=255)                                                                                              acaGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTa                                              >  5:245894/1‑148 (MQ=255)                                                                                                         tGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGTCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGGGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCg                                   >  4:49505/1‑148 (MQ=255)                                                                                                           ccAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGTCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGGGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGa                                 <  3:49505/148‑1 (MQ=255)                                                                                                            cAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAg                                >  3:95655/1‑148 (MQ=255)                                                                                                                    gCCTGGGAGGCGTGAGCCTACGGAAAGTGTCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTc                                    <  6:62730/136‑1 (MQ=255)                                                                                                                    gCCTGGGAGGCGTGAGCCTACGGAAAGTGTCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTc                                    >  5:62730/1‑136 (MQ=255)                                                                                                                    gCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCa                        >  4:193206/1‑148 (MQ=255)                                                                                                                     ccTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCaaa                      <  5:108123/149‑1 (MQ=255)                                                                                                                         ggAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAgg                  >  3:27698/1‑149 (MQ=255)                                                                                                                                     cTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGc      <  8:187847/149‑1 (MQ=255)                                                                                                                                         ggAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGAAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTgg  <  1:189849/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             TTGGCGGCGGAGAGTCGATCGGACAGTCGCTGCCTGTTTTTGTTCCGCAACTTTTTGTTGTGTGAGCAGCAGGGGGAGGAAAGTCCGGGCTCCACAGGGCGGAGTGCCAGGTAACGCCTGGGAGGCGTGAGCCTACGGAAAGTGCCACAGAAAATAACCGCCTAAGCACTTCGGTGCCGGTAAGGGTGAAAAGGTGCGGTAAGAGCGCACCGCACGGCTGGCAACAGTTCGTGGCTAGGTAAACCCCACTCGGAGCAAGACCAAATAGGGTTCCATCAGGCGTGG  >  NC_002947/1512686‑1512970

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