Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	527,898	Δ1 bp	15.9%	noncoding (1133/2902 nt)	PP_23SC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	527,898	0	C	.	15.9%	168.3 / 31.0	63	noncoding (1133/2902 nt)	PP_23SC	23S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base . (3/7);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATG  >  NC_002947/527799‑528032                                                                                                    |                                                                                                                                          cTTAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑Cga                                                                                                    <  2:208338/136‑1 (MQ=32)   tAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑Cgatg                                                                                                  <  1:530552/136‑1 (MQ=32)    aGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑Cgatga                                                                                                 >  1:461877/1‑136 (MQ=32)             ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑Cga                                                                                                    >  1:128240/1‑124 (MQ=17)             ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTaga                                                                                        >  1:291841/1‑136 (MQ=17)             ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGCCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑Cga                                                                                                    <  2:128240/124‑1 (MQ=17)             ggatgTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGGAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAACCCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTagg                                                                                        >  1:278156/1‑135 (MQ=12)             ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTaga                                                                                        >  1:56245/1‑136 (MQ=32)             ggaggTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTaga                                                                                        >  1:382349/1‑136 (MQ=32)                       ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTcg                                                                              >  1:122733/1‑135 (MQ=17)                       ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCt                                                                              >  2:50036/1‑136 (MQ=32)                       ttAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCt                                                                              >  1:501373/1‑136 (MQ=32)                          gAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                           >  1:13756/1‑136 (MQ=12)                          gAAGCAGCCATCATTTAAAGAAAGCGTAATAGCTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                           >  1:262167/1‑136 (MQ=12)                          gAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTa                                                                           <  2:359802/136‑1 (MQ=32)                            agcagcCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAg                                                                         >  1:278815/1‑136 (MQ=32)                                  cACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCtgtg                                                                   <  2:204495/136‑1 (MQ=32)                                   aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCt                                                                      <  1:301441/132‑1 (MQ=32)                                   aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCt                                                                      >  2:301441/1‑132 (MQ=32)                                   aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                  >  1:369344/1‑136 (MQ=32)                                   aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGa                                                                  >  1:277712/1‑136 (MQ=32)                                   aCCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCtgtg                                                                   >  2:301882/1‑136 (MQ=0)                                        ttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                             >  2:194095/1‑136 (MQ=32)                                        ttAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                             >  2:403853/1‑136 (MQ=32)                                        atAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTTGAGGAGCGTTCTGTAAGCCTGTGAAGGTg                                                             >  2:114539/2‑136 (MQ=17)                                            agaaagCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGtt                                                         <  2:461877/136‑1 (MQ=32)                                               aagAGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAATCCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGTTGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTgag                                                      <  2:11725/136‑1 (MQ=12)                                                  cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                   >  2:441514/1‑136 (MQ=32)                                                  cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                   >  2:102094/1‑136 (MQ=32)                                                  cGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAg                                                   >  1:483121/1‑136 (MQ=32)                                                        aGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                             <  2:278815/136‑1 (MQ=32)                                                        aGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCt                                             <  2:12824/136‑1 (MQ=32)                                                          cTCACTGGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTgg                                           <  2:396659/136‑1 (MQ=17)                                                               tAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTa                                      >  1:412347/1‑136 (MQ=32)                                                                agtcgagtcgGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTAt                                     >  1:512455/1‑136 (MQ=32)                                                                      gtcgGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAg                               <  1:323495/136‑1 (MQ=32)                                                                          gCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCg                           >  1:139107/1‑136 (MQ=32)                                                                          gCCTGCGCGGAAGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGc                            >  1:347324/1‑136 (MQ=0)                                                                          gCCTGCGCGGAAGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGc                            >  1:285840/1‑136 (MQ=0)                                                                            cTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGaa                         >  1:296085/1‑136 (MQ=32)                                                                                 cgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                    <  1:71651/136‑1 (MQ=32)                                                                                 cgGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTg                    <  1:320157/136‑1 (MQ=32)                                                                                  ggAAGATGTAACGGGGCT‑AAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                  <  2:136313/136‑1 (MQ=17)                                                                                   gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                  <  2:501373/136‑1 (MQ=32)                                                                                   gAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                  <  2:510769/136‑1 (MQ=32)                                                                                    aaGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGAc                  <  2:199546/136‑1 (MQ=0)                                                                                    aaGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGTACTGTGAGGTATGCTGGAGGTATCAGAAGTGCGAATGCTGAc                  <  2:29172/136‑1 (MQ=0)                                                                                    aaGATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGAc                  <  2:32206/136‑1 (MQ=0)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:281105/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:91363/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:512455/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:501083/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:277712/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:475210/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:291841/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:388595/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:382349/136‑1 (MQ=32)                                                                                     aGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAATGCGAATGCTGACAt                <  2:369344/136‑1 (MQ=17)                                                                                      gATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:347324/136‑1 (MQ=0)                                                                                      gATGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACAt                <  2:473222/136‑1 (MQ=0)                                                                                       ttGTAACGGGGCT‑AAACCATGCACCGAAGCTGCGGCAGCGACGCTTATGCGTTGTTG‑GGTAGGGGAGCGTTCTGTAAGCCTGCGAAGGTGTGCTGTGAGGCATGCTGGAGGTATCAGAAGTGCGAATGCTGACAta               <  2:109993/135‑2 (MQ=0)                                                                                                   tCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg  <  2:400312/136‑1 (MQ=32)                                                                                                   tCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATg  >  1:78645/1‑136 (MQ=32)                                                                                                    |                                                                                                                                          CTTAGACAGCTAGGAGGTTGGCTTAGAAGCAGCCACCCTTTAAAGAAAGCGTAATAGCTCACTAGTCGAGTCGGCCTGCGCGGAAGATGTAACGGGGCTCAAACCATGCACCGAAGCTACGGGTATCA‑TCTTA‑‑CGATGATGCGGTAGAGGAGCGTTCTGTAAGCCTGTGAAGGTGAGTTGAGAAGCTTGCTGGAGGTATCAGAAGTGCGAATGCTGACATGAGTAACGACAATG  >  NC_002947/527799‑528032

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