Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,065,730	T→A	6.6%	I39F (ATC→TTC)	PP_0923 ←	acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,065,730	0	T	A	6.6%	228.3 / 5.1	76	I39F (ATC→TTC)	PP_0923	acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (35/36);  new base A (3/2);  total (38/38)
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.34e-01

GCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCGTG  >  NC_002947/1065586‑1065878                                                                                                                                                 |                                                                                                                                                     gCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGgc                                                                                                                                                  <  3:214233/149‑1 (MQ=255)   aGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGgcgc                                                                                                                                                >  6:111362/1‑149 (MQ=255)        gATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTg                                                                                                                                            >  4:148218/1‑148 (MQ=255)         aTGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGaa                                                                                                                                          <  1:233675/149‑1 (MQ=255)          tGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCATCGACCAGGCGCTTCATGAACAGGCAGGTCGAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAg                                                                                                                                         <  1:34402/149‑1 (MQ=255)          tGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAg                                                                                                                                         <  3:170977/149‑1 (MQ=255)           gTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGc                                                                                                                                        >  8:290840/1‑149 (MQ=255)               gTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCaa                                                                                                                                    >  3:122197/1‑149 (MQ=255)                   aGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGt                                                                                                                                <  5:42128/149‑1 (MQ=255)                     gAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTc                                                                                                                              >  8:50163/1‑149 (MQ=255)                      aaaCATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCg                                                                                                                             <  4:122197/149‑1 (MQ=255)                       aaCATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGc                                                                                                                            >  3:49649/1‑149 (MQ=255)                         cATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGcc                                                                                                                           >  7:281722/1‑148 (MQ=255)                         cATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGcc                                                                                                                           <  8:281722/148‑1 (MQ=255)                          aTGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGccc                                                                                                                          >  1:319316/1‑148 (MQ=255)                             gTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAg                                                                                                                       <  7:57124/148‑1 (MQ=255)                                gCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCg                                                                                                                    >  7:10458/1‑148 (MQ=255)                                 cTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGa                                                                                                                   <  5:280266/148‑1 (MQ=255)                                      ccACAACATCGGCCCTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGCCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCCGCGCTGGCGCCTTTCGATCGGCGCCATGAAGCCCAAGCGTTCGCCCAAGGCCATAATGc                                                                                                             >  7:236109/1‑149 (MQ=255)                                        acaactTCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGAACAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTg                                                                                                           <  3:52510/149‑1 (MQ=255)                                         caacaTCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTgg                                                                                                          <  5:111362/149‑1 (MQ=255)                                            caTCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCa                                                                                                       >  8:315533/1‑149 (MQ=255)                                                  ccGCTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCa                                                                                                 <  2:68502/149‑1 (MQ=255)                                                       gCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAACATGcc                                                                                             <  7:315533/148‑1 (MQ=255)                                                              ttCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGcca                                                                                     <  3:52573/149‑1 (MQ=255)                                                                 gCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCt                                                                                                            <  5:272967/123‑1 (MQ=255)                                                                 gCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCt                                                                                                            >  6:272967/1‑123 (MQ=255)                                                                     aCTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGccagcaccag                                                                              <  3:148218/149‑1 (MQ=255)                                                                          cTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACGGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACcagcag                                                                           >  5:71075/1‑147 (MQ=255)                                                                          cTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACcagcagca                                                                         >  6:346622/1‑149 (MQ=255)                                                                          cTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACcagcagca                                                                         >  4:97863/1‑149 (MQ=255)                                                                          cTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACcagcagca                                                                         >  7:307957/1‑149 (MQ=255)                                                                                 tCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGAt                                                                  <  2:303627/149‑1 (MQ=255)                                                                                        gcagggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCg                                                                     >  7:279895/1‑139 (MQ=255)                                                                                        gcagggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCg                                                                     <  8:279895/139‑1 (MQ=255)                                                                                          agggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGc                                                         >  6:229253/1‑149 (MQ=255)                                                                                            ggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGcc                                                        >  1:120781/1‑148 (MQ=255)                                                                                            ggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCa                                                       <  5:346622/149‑1 (MQ=255)                                                                                            ggcagCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCa                                                       <  2:120781/149‑1 (MQ=255)                                                                                                gCGACCAGGCGCTTGATGAACAGGCAGGGCCAGCGCTGGCGCCGTTCGATCGGTGCCTTGAAGCCCAAGCGTTCGCCCATGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCGAGTAc                                                   <  8:236109/149‑1 (MQ=255)                                                                                                     cAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGc                                               >  8:65528/1‑148 (MQ=255)                                                                                                        gcgcTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACt                                           >  7:237991/1‑149 (MQ=255)                                                                                                           cTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCg                                        >  1:151288/1‑149 (MQ=255)                                                                                                              cATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGTCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGc                                      <  7:255057/148‑1 (MQ=255)                                                                                                              cATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCgg                                       <  2:151288/147‑1 (MQ=255)                                                                                                               aTGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCAt                                    >  3:147426/1‑149 (MQ=255)                                                                                                                tGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCAt                                    >  5:103619/1‑148 (MQ=255)                                                                                                                tGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGACGCGGGCCAGTACCAGCAACATCGGCAt                                    >  7:140485/1‑148 (MQ=255)                                                                                                                  aaCAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCa                                  <  6:103619/148‑1 (MQ=255)                                                                                                                       gcaggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGAccgc                            >  4:119663/1‑149 (MQ=255)                                                                                                                       gcaggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGAccgc                            >  4:114769/1‑149 (MQ=255)                                                                                                                         aggTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGAccgccg                          >  2:340471/1‑149 (MQ=255)                                                                                                                            tCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCt                       >  3:41172/1‑149 (MQ=255)                                                                                                                            tCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCt                       >  3:36744/1‑149 (MQ=255)                                                                                                                             ccAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGTATCAGACCGCCGCCtt                      >  4:163051/1‑149 (MQ=255)                                                                                                                              cAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGTATCAGACCGCCGCCtt                      <  3:163051/148‑1 (MQ=255)                                                                                                                               aGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTg                     >  6:273438/1‑148 (MQ=255)                                                                                                                                ttcaaggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTg                     >  4:122504/6‑147 (MQ=255)                                                                                                                                gcgcTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGaa                   <  7:290840/149‑1 (MQ=255)                                                                                                                                gcgcTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGaa                   <  6:71075/149‑1 (MQ=255)                                                                                                                                   aaggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTg                >  4:65916/3‑149 (MQ=255)                                                                                                                                    tGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTgg               >  1:127275/1‑149 (MQ=255)                                                                                                                                    tGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTgg               <  3:114769/149‑1 (MQ=255)                                                                                                                                     ggCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggc              <  7:65528/149‑1 (MQ=255)                                                                                                                                     ggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGaa                   >  1:73325/1‑144 (MQ=255)                                                                                                                                     ggCGCCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGaa                   <  2:73325/144‑1 (MQ=255)                                                                                                                                      gcgcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggc              <  5:229253/148‑1 (MQ=255)                                                                                                                                      gcgcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTgg               <  3:119663/147‑1 (MQ=255)                                                                                                                                       cgcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgg            <  4:41172/149‑1 (MQ=255)                                                                                                                                       cgcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgg            <  4:36744/149‑1 (MQ=255)                                                                                                                                        gcCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcggg           <  4:49649/149‑1 (MQ=255)                                                                                                                                        gcCGATCGAACGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcggg           <  3:65916/149‑1 (MQ=255)                                                                                                                                          cGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgggcg         >  3:160706/1‑149 (MQ=255)                                                                                                                                          cGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTggcgggcg         <  2:127275/149‑1 (MQ=255)                                                                                                                                               gATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGc     <  8:237991/148‑1 (MQ=255)                                                                                                                                               gATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCg    >  3:249611/1‑149 (MQ=255)                                                                                                                                                 tCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGGCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCGt   <  8:10458/148‑1 (MQ=255)                                                                                                                                                 tCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCGTg  <  8:307957/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GCAGGGGGATGTCGGTCCAGGAAACATGGTTGCTGACCCACAACATCGGCCGTTGCGGAAGTTCGCCGACTACCTTGACGTCGAAGGGCAGGGCAGCGACCAGGCGCTTCATGAACAGGCAGGTCCAGCGCTGGCGCCGTTCGATCGGCGCCTTGAAGCCCAAGCGTTCGCCCAAGGCGATAATGCTGGCCATCAGCATGCCCTGCGCCAGCACCAGCAGCAGTCGGATGAGGCGGGCCAGTACCCGCAACTTCGGCATCAGACCGCCGCCTTGAAGTGGCGGGCGTAGCGTG  >  NC_002947/1065586‑1065878

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