Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,447,895	T→C	5.5%	Q1063R (CAG→CGG)	PP_2146 ←	SNF2/RAD54 family helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,447,895	0	T	C	5.5%	213.5 / 2.1	72	Q1063R (CAG→CGG)	PP_2146	SNF2/RAD54 family helicase
Reads supporting (aligned to +/- strand):  ref base T (28/40);  new base C (2/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

GCTCTATCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACGGTGTCCGCCGCCGT  >  NC_002947/2447747‑2448043                                                                                                                                                     |                                                                                                                                                        gCTCTATCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACt                                                                                                                                                         <  3:345069/149‑1 (MQ=255)     tATCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGa                                                                                                                                                      >  5:311282/1‑148 (MQ=255)       tCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGAtct                                                                                                                                                   >  1:171319/1‑149 (MQ=255)       tCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGAtct                                                                                                                                                   >  2:210565/1‑149 (MQ=255)       tCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGAtc                                                                                                                                                    <  1:187993/148‑1 (MQ=255)        cccccTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGAtctt                                                                                                                                                  >  7:1263/1‑149 (MQ=255)         ccccATTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACCGGGCCTCGATGTCTACATCGCCAAGCCGCCATTGCCCACCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATCCTGCGTATTCCTGCTGCAACAACTGGAtcttc                                                                                                                                                 >  5:322868/1‑149 (MQ=255)            cTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGa                                                                                                                                                      <  2:37597/141‑1 (MQ=255)            cTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGAtagatcg                                                                                                                                               >  1:37597/1‑142 (MQ=255)               tGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGAtcttctctt                                                                                                                                             <  8:83338/147‑1 (MQ=255)                 gCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACCGGATCTTCTCTT‑CCAc                                                                                                                                        >  8:249863/1‑149 (MQ=255)                      ccgcccgcGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCAAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACCACTGGATCTTCTCTT‑CCACCGTAc                                                                                                                                   >  6:383313/1‑149 (MQ=255)                         cccgcGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTAcccc                                                                                                                                >  6:64028/1‑149 (MQ=255)                          ccgcGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCg                                                                                                                               <  1:107790/149‑1 (MQ=255)                            gcgcTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCggg                                                                                                                             <  4:318588/149‑1 (MQ=255)                                    cGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCggg                                                                                                                             <  4:38331/141‑1 (MQ=255)                                    cGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCggg                                                                                                                             >  3:38331/1‑141 (MQ=255)                                    cGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGc                                                                                                                   <  8:1263/149‑1 (MQ=255)                                     gggCAATGGCGCGTACAGGTCCTCGATTTCTTCATCGCCAAGCCGGCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTCTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCTCCGTACCCCGGGTA‑A‑TCAGCt                                                                                                                  <  6:322868/149‑1 (MQ=255)                                     gggCAATGGCGCGAACAGGGCGTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCt                                                                                                                  >  5:239618/1‑149 (MQ=255)                                     gggCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACATCTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCt                                                                                                                  <  3:280633/149‑1 (MQ=255)                                      ggCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCtt                                                                                                                 <  8:147177/149‑1 (MQ=255)                                      ggCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCtt                                                                                                                 <  3:336072/149‑1 (MQ=255)                                      ggCAAAGGCGCGAACCTGGCCGCCCTTTCATCCTCGCCCAGCCGCCATTGCCCACGCTGCCCGCCGTCCATCATGCGTGCGGCCAAATCTGCCTTTTCCTTGTGCAACCACTTGCTCTTCTCCT‑CCACCGTACCCCGGGTA‑A‑TCAGCtt                                                                                                                 <  8:176903/149‑1 (MQ=255)                                         aaTGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTgaa                                                                                                              <  7:306728/149‑1 (MQ=255)                                         aaTGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTgaa                                                                                                              <  1:92239/149‑1 (MQ=255)                                               gcgAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTgaacacgaa                                                                                                        <  2:171319/149‑1 (MQ=255)                                                  aaCAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTgaacacgaacac                                                                                                     >  1:350358/1‑149 (MQ=255)                                                  aaCAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTgaacacgaacac                                                                                                     <  1:210565/149‑1 (MQ=255)                                                      gggCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGt                                                                                                 <  1:366775/149‑1 (MQ=255)                                                      gggCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGt                                                                                                 <  5:139661/149‑1 (MQ=255)                                                         ccTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTg                                                                                              >  2:9997/1‑149 (MQ=255)                                                         ccTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACAGGTTTg                                                                                              >  1:205957/1‑149 (MQ=255)                                                              aTTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTg                                                                                         <  6:311282/149‑1 (MQ=255)                                                              aTTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTg                                                                                         <  6:239618/149‑1 (MQ=255)                                                               tttCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTTCCCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTg                                                                                         <  1:274901/149‑1 (MQ=255)                                                                 tcttcATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGc                                                                                       <  4:308176/148‑1 (MQ=255)                                                                  cttcATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGcc                                                                                      <  5:64028/148‑1 (MQ=255)                                                                           cAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACCGGATCTTCTCTT‑CCACCGTAcccc                                                                                                                                >  1:61980/1‑99 (MQ=255)                                                                           cAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACCGGATCTTCTCTT‑CCACCGTAcccc                                                                                                                                <  2:61980/99‑1 (MQ=255)                                                                            aaGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAgg                                                                           >  1:214161/1‑149 (MQ=255)                                                                             aGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCACGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAgg                                                                           >  8:276073/1‑148 (MQ=255)                                                                                cgccATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCCCCGTACCCCGCGAA‑AGGGGGCGAGAACACGCACACCGGGTTTTGCAGGGCGCTGGGGTAgggcc                                                                        >  1:315809/1‑146 (MQ=255)                                                                                cgccATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTACA‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAg                                                                            <  1:269628/145‑1 (MQ=255)                                                                                cgccATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAg                                                                            >  2:269628/1‑144 (MQ=255)                                                                                 gccATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTa                                                                             >  8:342110/1‑142 (MQ=255)                                                                                 gccATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTa                                                                             <  7:342110/142‑1 (MQ=255)                                                                                   cATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCAcggtc                                                                    >  8:48585/1‑149 (MQ=255)                                                                                    aTTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTTGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCAcggtcg                                                                   <  2:205957/149‑1 (MQ=255)                                                                                    aTTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACCGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCAcggtcg                                                                   <  7:249863/149‑1 (MQ=255)                                                                                     ttGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCAcggtcgg                                                                  <  2:350358/149‑1 (MQ=255)                                                                                        cccAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGCTCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTa                                                                <  2:243999/148‑1 (MQ=255)                                                                                         ccAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGc                                                              >  8:18934/1‑149 (MQ=255)                                                                                                   cgccgTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTtcg                                                    <  7:18934/149‑1 (MQ=255)                                                                                                      cgTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTtcgctc                                                 <  2:214161/149‑1 (MQ=255)                                                                                                                      ggcCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCggg                                           <  4:21234/139‑1 (MQ=255)                                                                                                                      ggcCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCggg                                           >  3:21234/1‑139 (MQ=255)                                                                                                                       gcCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCAcggtcgg                                                                  >  6:200138/1‑115 (MQ=255)                                                                                                                       gcCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCAcggtcgg                                                                  <  5:200138/115‑1 (MQ=255)                                                                                                                       gcCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAggg                                <  1:9997/149‑1 (MQ=255)                                                                                                                        ccAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAgggg                               >  5:171607/1‑149 (MQ=255)                                                                                                                         cAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGt                              <  7:251184/149‑1 (MQ=255)                                                                                                                          aaTGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTc                             <  5:383313/149‑1 (MQ=255)                                                                                                                               tgcCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTc                                        <  6:221580/133‑1 (MQ=255)                                                                                                                               tgcCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTc                                        >  5:221580/1‑133 (MQ=255)                                                                                                                                 ccTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGt                        <  6:171607/147‑1 (MQ=255)                                                                                                                                 ccTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCACCGGGTTCCACCAGGGGTCGAAg                         <  5:174364/146‑1 (MQ=255)                                                                                                                                 ccTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCACCGGGTTCCACCAGGGGTCGAAg                         >  6:174364/1‑146 (MQ=255)                                                                                                                                        ctgctgCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACg                <  7:48585/148‑1 (MQ=255)                                                                                                                                             tcataaacTGGATCTTCTCTT‑CCAGCGTACCCCGTGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCTATGCGGTAGGCACTGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACGGTGTcc          <  4:260175/144‑1 (MQ=255)                                                                                                                                              caacaaCTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGAGGTAGTCACGGTCGGTAGCCTGTTTTTCGCTAGCCTGGTTCCACCAGGGGTCGAAGTGGATCACCGTTTccg         >  2:142803/1‑149 (MQ=255)                                                                                                                                               aacaacTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAAAAACGGTTTGTCCTGGCCGATGCGGTAGGAACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACGGTGTccga        >  3:260175/1‑148 (MQ=255)                                                                                                                                                acaagTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTTGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACGGTCTccgcc       >  8:293775/1‑149 (MQ=255)                                                                                                                                                     tGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACGGTGTCCGCCGCCGt  >  4:107299/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                        GCTCTATCCCCCTTTGGCCTACCGCCCGCGCTTCCCGGGCAATGGCGCGAACAGGGCCTCGATTTCTTCATCGCCAAGCCGCCATTGCCCAGCCTGCCCGCCGTCCAGCAGGCTGGCGGCCAATGCTGCCTTTTCCTGCTGCAACAACTGGATCTTCTCTT‑CCACCGTACCCCGGGTA‑A‑TCAGCTTGAACACGAACACCGGTTTGTCCTGGCCGATGCGGTAGGCACGGTCGGTAGCCTGGTTTTCGCTCGCCGGGTTCCACCAGGGGTCGAAGTGGATCACGGTGTCCGCCGCCGT  >  NC_002947/2447747‑2448043

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