Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,274,268	T→A	5.3%	H330L (CAC→CTC)	cioB ←	cyanide insensitive ubiquinol oxidase subunit II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,274,268	0	T	A	5.3%	225.4 / 2.7	75	H330L (CAC→CTC)	cioB	cyanide insensitive ubiquinol oxidase subunit II
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGACGGCGGGAT  >  NC_002947/5274122‑5274410                                                                                                                                                    |                                                                                                                                                gCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTgg                                                                                                                                               >  1:177333/1‑149 (MQ=255) gCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTgg                                                                                                                                               <  8:336312/149‑1 (MQ=255) gCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTgg                                                                                                                                               >  8:15276/1‑149 (MQ=255)  cATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTggg                                                                                                                                              <  6:166228/149‑1 (MQ=255)      acacaGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCa                                                                                                                                           >  1:145416/1‑148 (MQ=255)           ggCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGc                                                                                                                                     <  2:177333/149‑1 (MQ=255)                ccACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCcgcg                                                                                                                                 >  5:385470/1‑148 (MQ=255)                   cacCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAac                                                                                                                             <  2:4242/149‑1 (MQ=255)                   cacCCAGCGCCGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAac                                                                                                                             <  7:28662/149‑1 (MQ=255)                      ccAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACg                                                                                                                          >  5:377274/1‑149 (MQ=255)                        aGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACg                                                                                                                          >  7:329353/1‑147 (MQ=255)                               cGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGc                                                                                                                   <  1:9583/147‑1 (MQ=255)                                       aTCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAgg                                                                                                          >  5:225157/1‑148 (MQ=255)                                       aTCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAg                                                                                                           <  4:268552/147‑1 (MQ=255)                                       aTCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAg                                                                                                           <  6:304701/147‑1 (MQ=255)                                       aTCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGACTTCGTGGGTCACCTTGCCGCGGAACACCTAGTAGCTCCAGAAgg                                                                                                          >  5:200669/1‑148 (MQ=255)                                       aTCGCCCAGATCAGCAACACCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGt                                                                                                         <  6:385470/149‑1 (MQ=255)                                          gCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGtat                                                                                                       <  3:205860/148‑1 (MQ=255)                                               gATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACAt                                                                                                   >  2:24817/1‑147 (MQ=255)                                               gATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACAt                                                                                                   >  3:263400/1‑147 (MQ=255)                                               gATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGa                                                                                                 <  4:194617/149‑1 (MQ=255)                                                cacAGCAACAGGCAACCCAGGCTCTGCCACATCGGTTTCTCGTCTTCCACACGCTGGTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCTTCTTCGTGGGGCACCTTGCCGCGGAACACGTAGTAGCTCCAGAATGTATACCTGAg                                                                                                <  4:237066/147‑1 (MQ=255)                                                aTCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAg                                                                                                <  2:145416/149‑1 (MQ=255)                                                   aGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGa                                                                                              <  4:109687/148‑1 (MQ=255)                                                          gCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACgg                                                                                       >  7:138521/1‑148 (MQ=255)                                                           ccAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGaa                                                                                                            <  5:283424/126‑1 (MQ=255)                                                           ccAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGaa                                                                                                            >  6:283424/1‑126 (MQ=255)                                                           ccAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACgg                                                                                       >  4:732/1‑147 (MQ=255)                                                           ccAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCa                                                                                     <  6:377274/149‑1 (MQ=255)                                                             aaCCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAg                                                                                    >  3:247129/1‑148 (MQ=255)                                                              aCCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTgg                                                                                                                                               >  4:340179/1‑88 (MQ=255)                                                              aCCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGa                                                                                  >  3:137815/1‑149 (MQ=255)                                                               cccAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGAt                                                                                 >  8:78696/1‑149 (MQ=255)                                                                ccAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGAGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGaa                                                                                                            <  7:64226/121‑1 (MQ=255)                                                                ccAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGAGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGaa                                                                                                            <  7:69736/121‑1 (MQ=255)                                                                ccAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGAGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGaa                                                                                                            >  8:69736/1‑121 (MQ=255)                                                                ccAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGAGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGaa                                                                                                            >  8:64226/1‑121 (MQ=255)                                                                     cGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATaa                                                                                           >  6:245850/1‑133 (MQ=255)                                                                     cGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATaa                                                                                           <  5:245850/133‑1 (MQ=255)                                                                         gCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGt                                                                                                         <  4:224770/115‑1 (MQ=255)                                                                         gCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGt                                                                                                         >  3:224770/1‑115 (MQ=255)                                                                          ccACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTAc                                                                      >  4:265398/1‑149 (MQ=255)                                                                              aGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCa                                                                                     <  1:8389/130‑1 (MQ=255)                                                                              aGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCa                                                                                     >  2:8389/1‑130 (MQ=255)                                                                              aGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACccac                                                                  <  7:15276/149‑1 (MQ=255)                                                                                cGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACccacca                                                                <  4:247129/149‑1 (MQ=255)                                                                                  gTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACa                                                                            >  1:360858/1‑135 (MQ=255)                                                                                  gTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACa                                                                            <  2:360858/135‑1 (MQ=255)                                                                                         gTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGaa                                                         >  5:338798/1‑147 (MQ=255)                                                                                         gTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGaa                                                         <  6:338798/147‑1 (MQ=255)                                                                                             tCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTg                                                   >  1:168075/1‑149 (MQ=255)                                                                                              ccAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCTGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCt                                                    <  6:225157/147‑1 (MQ=255)                                                                                               cAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAGCCCGTGTTCGTGGGTCACATTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGCATGAACAGGGTACCCACCAGCATCAAGCCATg                                                   >  1:288474/1‑148 (MQ=255)                                                                                                aaaCCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCt                                                >  6:410084/1‑149 (MQ=255)                                                                                                 aaCCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAg                                                        >  8:364755/1‑140 (MQ=255)                                                                                                 aaCCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAg                                                        <  7:364755/140‑1 (MQ=255)                                                                                                        cTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTg                                        <  3:265398/149‑1 (MQ=255)                                                                                                              agccaTCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑g                                 <  2:288474/149‑1 (MQ=255)                                                                                                                ccaTCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCg                               <  4:137815/149‑1 (MQ=255)                                                                                                                 catcatGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCgg                              >  2:208976/1‑149 (MQ=255)                                                                                                                  atcatGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGaa                                                                              >  6:366576/1‑101 (MQ=255)                                                                                                                  atcatGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGaa                                                                              <  5:366576/101‑1 (MQ=255)                                                                                                                   tcatGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCa                            >  7:36206/1‑149 (MQ=255)                                                                                                                     atGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGGCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACCGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATc                          <  8:329353/149‑1 (MQ=255)                                                                                                                             ctcctAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATc                   >  5:145084/1‑148 (MQ=255)                                                                                                                                ctAGTGATAG‑CTGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTATCTCCAGAAGTTTTACATGAGGATAAACGGCAGGATGAACAGTTTACCGACCAGCATGAAGCCCTGTCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGaaa               <  6:200669/149‑1 (MQ=255)                                                                                                                                 tAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCAGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGTATAAAAGGCAGGATGAACAGGGTACACACCAGCATGAAGCCCTGTCTTTGAGGTGGCGCT‑GCg                               >  7:229343/1‑132 (MQ=255)                                                                                                                                   gTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCg            <  2:168075/149‑1 (MQ=255)                                                                                                                                      aTAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGAcg         >  1:298838/1‑149 (MQ=255)                                                                                                                                      aTAG‑CCGTCTTCCTGGCTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATCAGGATAAACCCCAGGATGAACAGCGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCACATGCAAATCTAcg         >  6:242175/1‑149 (MQ=255)                                                                                                                                        aG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGAcggc       >  1:276412/1‑149 (MQ=255)                                                                                                                                         g‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGAcggcg      <  8:36206/149‑1 (MQ=255)                                                                                                                                           ccGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCATAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTAGCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCTAGCGGCATCCCAGATCGAAATCGAcggcg      <  1:24817/149‑1 (MQ=255)                                                                                                                                              tCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGACGGCGGgat  >  6:177366/1‑149 (MQ=255)                                                                                                                                              tCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGACGGCGGgat  >  4:69541/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                GCATCACACAGGCCGCCACACCCAGCGCTGCGACGCTCATCGCCCAGATCAGCAACAGCCAACCCAGTCGCTGCCACAGCGGTTTCTTGTCTTCCAAACCATGCTTGCCAGCCATCATGCAGCCCTCCTAGTGATAG‑CCGTCTTCGTGGGTCACCTTGCCGCGGAACACGTAGTAGCTCCAGAAGGTATACATGAGGATAAACGGCAGGATGAACAGGGTACCCACCAGCATGAAGCCCTGGCTTTGCGGTGGCGCT‑GCGGCATCCCAGATCGAAATCGACGGCGGGAT  >  NC_002947/5274122‑5274410

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