Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,504,068	T→A	11.5%	intergenic (+35/‑228)	rpsI → / → petA	30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,504,068	0	T	A	11.5%	160.7 / 11.5	61	intergenic (+35/‑228)	rpsI/petA	30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (30/24);  new base A (3/4);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGA  >  NC_002947/1503920‑1504214                                                                                                                                                      |                                                                                                                                                    cTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                     <  7:276000/149‑1 (MQ=255)      tGGCGCTCTGCGTCGTGCTGGCTTCGTCACCCGCGACGATCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑ACCGTGAATACCGCTTCGGCAGTCGAAAAAAGCCCGCTTCGTTTGaag                                                                                                                                                >  8:239260/1‑148 (MQ=255)         cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                             >  3:934/1‑149 (MQ=255)          gcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCggg                                                                                                                                            >  5:333010/1‑149 (MQ=255)           ctctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGc                                                                                                                                           <  8:16627/149‑1 (MQ=255)           ctctGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGACAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGc                                                                                                                                           <  2:154512/149‑1 (MQ=255)              tGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑t                                                                                                                                         <  2:21973/147‑1 (MQ=255)                cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑ttttt                                                                                                                                     >  4:37716/1‑149 (MQ=255)                cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGC‑tt                                                                                                                                        <  6:333010/146‑1 (MQ=255)                  tcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGC‑ttttttc                                                                                                                                   <  6:219745/149‑2 (MQ=255)                   cgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑tttttttt                                                                                                                                  <  1:24648/149‑1 (MQ=255)                       cTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGTGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGc                                                                                                                              <  1:40692/149‑1 (MQ=255)                       cTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGc                                                                                                                              <  1:148940/149‑1 (MQ=255)                          gCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTGGcgt                                                                                                                            >  3:251391/1‑146 (MQ=255)                                       cgACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTc                                                                                                              <  3:127079/149‑1 (MQ=255)                                                tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                     <  6:33988/102‑1 (MQ=255)                                                tGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                     >  5:33988/1‑102 (MQ=255)                                                            tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                         >  7:1015/1‑149 (MQ=255)                                                             aaGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGc                                                                                                                                           <  6:20972/99‑1 (MQ=255)                                                               gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGaa                                                                                                                        >  6:205802/1‑115 (MQ=255)                                                               gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCgaa                                                                                      >  3:201764/1‑147 (MQ=255)                                                                 aaGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAGAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                      >  1:229143/1‑148 (MQ=255)                                                                       gTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTa                                                                              <  2:150989/149‑1 (MQ=255)                                                                              tAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                       <  2:229143/149‑1 (MQ=255)                                                                                     cgTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                       >  2:322020/1‑142 (MQ=255)                                                                                     cgTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAg                                                                       <  1:322020/142‑1 (MQ=255)                                                                                                tCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCtt                                                                                                                           >  3:163861/1‑79 (MQ=255)                                                                                                tCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCtt                                                                                                                           <  4:163861/79‑1 (MQ=255)                                                                                                tCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGc                                                     >  5:116833/1‑149 (MQ=255)                                                                                                      cTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGc                                               <  3:37716/149‑1 (MQ=255)                                                                                                       tCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCaa                                                                        >  6:305189/1‑123 (MQ=255)                                                                                                       tCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCaa                                                                        <  5:305189/123‑1 (MQ=255)                                                                                                       tCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGcc                                              >  4:224397/1‑149 (MQ=255)                                                                                                       tCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGcc                                              >  4:220135/1‑149 (MQ=255)                                                                                                       tCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGcc                                              >  7:208912/1‑149 (MQ=255)                                                                                                       tCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGATCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAATCAAAGCGTGGGTTGTCGCCCTTGcc                                              >  5:330567/1‑149 (MQ=255)                                                                                                         cA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTgg                                                             <  4:62653/132‑1 (MQ=255)                                                                                                         cA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTgg                                                             >  3:62653/1‑132 (MQ=255)                                                                                                         cA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgg                                            >  2:232112/1‑149 (MQ=255)                                                                                                          a‑aGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                           >  5:207623/1‑149 (MQ=255)                                                                                                            aGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAt                                                                           <  5:56865/116‑1 (MQ=255)                                                                                                            aGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAt                                                                           >  6:56865/1‑116 (MQ=255)                                                                                                              cGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGt                                         >  2:6391/1‑148 (MQ=255)                                                                                                               gTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTTGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTaa                                       <  4:934/149‑1 (MQ=255)                                                                                                                    tACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGCTTGTCGCCCTTGCCGGGGTAATCAcc                                  >  2:903/1‑149 (MQ=255)                                                                                                                     aCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTGAGCAGTCAGTGGCCGCGTGTGGCATAGACGTCGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTTATAACCt                                 >  2:255151/1‑149 (MQ=255)                                                                                                                     aCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGcagt                                                                                                        >  4:172925/1‑78 (MQ=255)                                                                                                                     aCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGcagt                                                                                                        <  3:172925/78‑1 (MQ=255)                                                                                                                      ccGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCaaa                                                                   <  8:19293/115‑1 (MQ=255)                                                                                                                      ccGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCaaa                                                                   >  7:19293/1‑115 (MQ=255)                                                                                                                            cGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGGCGCCCTTGCCGGGGTATTCACCTTGTGAg                           >  7:251037/1‑149 (MQ=255)                                                                                                                            cGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgt                          >  8:150709/1‑149 (MQ=255)                                                                                                                              gCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCCAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgt                        <  7:150709/149‑1 (MQ=255)                                                                                                                              gCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgt                        <  1:6391/149‑1 (MQ=255)                                                                                                                                    gAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAACAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATAAACTTGTCAGTGTGTGGAt                  >  6:21313/1‑149 (MQ=255)                                                                                                                                     aaaaaaGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGAtt                 >  2:196211/1‑149 (MQ=255)                                                                                                                                         aaGCCCGGTTCCTTGGAACCGGGA‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGtt             <  8:208912/149‑1 (MQ=255)                                                                                                                                           gCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAACAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCTCCTTGTCAGTGTGTGGATTTGTTTc           <  5:21313/149‑1 (MQ=255)                                                                                                                                           gCCCGGTTCCAAGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                            <  6:270719/84‑1 (MQ=255)                                                                                                                                           gCCCGGTTCCAAGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGa                                                                            >  5:270719/1‑84 (MQ=255)                                                                                                                                           gCCCGGTTCCAAGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTc           >  8:282438/1‑149 (MQ=255)                                                                                                                                               ggTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTAc       <  7:239260/149‑1 (MQ=255)                                                                                                                                                gTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTAcc      <  1:232112/149‑1 (MQ=255)                                                                                                                                                    cTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGa  >  5:146946/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCA‑AGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGC‑TTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGA  >  NC_002947/1503920‑1504214

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