Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,068,527:1	+G	7.9%	coding (693/948 nt)	PP_0926 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,068,527	1	.	G	7.9%	241.3 / 16.1	63	coding (693/948 nt)	PP_0926	membrane protein
Reads supporting (aligned to +/- strand):  ref base . (30/28);  new base G (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CCATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCGGGCAGC  >  NC_002947/1068379‑1068673                                                                                                                                                      |                                                                                                                                                   ccATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAt                                                                                                                                                     <  1:198168/149‑1 (MQ=255)  cATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGG‑CATCGAT‑gg                                                                                                                                                  <  1:98993/149‑2 (MQ=255)       cTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGc                                                                                                                                              >  5:40851/1‑149 (MQ=255)       cTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGc                                                                                                                                              >  2:94809/1‑149 (MQ=255)        ttACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCt                                                                                                                                             >  8:171916/1‑149 (MQ=255)          aCGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGG‑CATCGATGGCCAGccg                                                                                                                                            <  8:44239/148‑3 (MQ=255)           cGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGG‑CATCGATGGCCAGccga                                                                                                                                           <  6:19656/148‑4 (MQ=255)            ggCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGt                                                                                                                                         <  1:82581/149‑1 (MQ=255)            ggCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGACCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGt                                                                                                                                         >  3:137492/1‑149 (MQ=255)            ggCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGACCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGt                                                                                                                                         <  4:137492/149‑1 (MQ=255)              ccAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGa                                                                                                                                       <  7:171916/149‑1 (MQ=255)              ccAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGG‑CATCGATGGCCAGc                                                                                                                                              <  4:50303/142‑1 (MQ=255)              ccAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGG‑CATCGATGGCCAGc                                                                                                                                              >  3:50303/1‑142 (MQ=255)                 ggTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGaaaa                                                                                                                                    <  1:291221/149‑1 (MQ=255)                     ccGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCCGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGct                                                                                                                                <  6:9431/149‑1 (MQ=255)                       gCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGctct                                                                                                                              <  6:40851/149‑1 (MQ=255)                           cccAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTgg                                                                                                                            <  2:198054/147‑1 (MQ=255)                             cAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGc                                                                                                                        <  6:31791/149‑1 (MQ=255)                                 tGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGAAACTGTGGTCGGCTGGCCATCGAT‑GCAAGCTTCGTGAAAATGCTCTGGCTGCCgct                                                                                                                    >  8:12161/1‑149 (MQ=255)                                          aCCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGc                                                                                                             <  4:272106/147‑1 (MQ=255)                                                        ccTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGt                                                                                             >  7:104625/1‑149 (MQ=255)                                                        ccTGGCACAGATGACCCTGCCACTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGt                                                                                             >  1:287438/1‑149 (MQ=255)                                                         cTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                            >  1:78380/1‑149 (MQ=255)                                                         cTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                            >  1:247530/1‑149 (MQ=255)                                                          tGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTtg                                                                                           >  5:296641/1‑149 (MQ=255)                                                              accgATGACCCTGCGGCTGGCGGTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCg                                                                                       <  2:247530/146‑1 (MQ=255)                                                              acaGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCg                                                                                       >  5:291128/1‑149 (MQ=255)                                                              acaGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCg                                                                                       >  5:291126/1‑149 (MQ=255)                                                                      cccTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgc                                                                               <  8:104625/149‑1 (MQ=255)                                                                      cccTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCGGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgc                                                                               >  2:91056/1‑149 (MQ=255)                                                                        cTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCgg                                                                             >  1:21358/1‑149 (MQ=255)                                                                        cTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCgg                                                                             <  1:94809/149‑1 (MQ=255)                                                                        cTGCCACTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCgg                                                                             <  2:287438/149‑1 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCgg                                                                                      <  4:31004/139‑1 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCgg                                                                                      >  3:31004/1‑139 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCggg                                                                            >  7:244082/1‑149 (MQ=255)                                                                            cGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCg                                                                         >  5:280624/1‑149 (MQ=255)                                                                               tGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAg                                                                      <  6:291128/149‑1 (MQ=255)                                                                                 gcgcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCTGAg                                                                      <  6:291126/147‑1 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                                  <  3:238415/149‑1 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATc                                                             <  2:57618/145‑1 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGAGGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATc                                                             <  7:12161/145‑1 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCCGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCtt                                                          <  1:91056/148‑1 (MQ=255)                                                                                               cggcggCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                       <  2:21358/148‑1 (MQ=255)                                                                                                           gTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                       >  7:83482/1‑136 (MQ=255)                                                                                                           gTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                       <  8:83482/136‑1 (MQ=255)                                                                                                            tCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCa                                          <  6:280624/148‑1 (MQ=255)                                                                                                               cTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCgg                                                                                      <  1:47972/101‑1 (MQ=255)                                                                                                               cTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCgg                                                                                      >  2:47972/1‑101 (MQ=255)                                                                                                                       gcTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCt                                                     <  6:268971/126‑1 (MQ=255)                                                                                                                       gcTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GACAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCt                                                     >  5:268971/1‑126 (MQ=255)                                                                                                                           cGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCg                            <  2:78380/147‑1 (MQ=255)                                                                                                                             tAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCa                        <  8:244082/149‑1 (MQ=255)                                                                                                                              aaCTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCGTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCGTGTCCAGATGACGCCCGACCGCTGCGGcctg                       >  8:223670/1‑147 (MQ=255)                                                                                                                               aCTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGc                      >  4:190208/1‑149 (MQ=255)                                                                                                                                 tgtgGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTa                    >  7:27496/1‑149 (MQ=255)                                                                                                                                    acgaagCTGG‑CATCGATGGCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGt                 >  5:213935/6‑149 (MQ=255)                                                                                                                                        ggCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTcc                                                      <  7:234080/108‑1 (MQ=255)                                                                                                                                        ggCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTcc                                                      >  8:234080/1‑108 (MQ=255)                                                                                                                                        ggCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATg             >  7:276851/1‑149 (MQ=255)                                                                                                                                        ggCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATg             >  5:44078/1‑149 (MQ=255)                                                                                                                                        ggCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCAGCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATg             >  7:48641/1‑149 (MQ=255)                                                                                                                                         gCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATgg            >  7:155907/1‑149 (MQ=255)                                                                                                                                                aTCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCggg      >  8:287980/1‑148 (MQ=255)                                                                                                                                                    aT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCGGGCAGc  >  6:35063/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CCATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGTCGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCGGGCAGC  >  NC_002947/1068379‑1068673

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