Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,869,613	G→C	5.1%	G41R (GGG→CGG)	PP_5144 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,869,613	0	G	C	5.1%	257.3 / 3.1	78	G41R (GGG→CGG)	PP_5144	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/36);  new base C (2/2);  total (40/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 = 9.15e-01

TGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTG  >  NC_002947/5869466‑5869759                                                                                                                                                     |                                                                                                                                                   tGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTgg                                                                                                                                                   >  5:265552/1‑149 (MQ=255)   aCCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTgggg                                                                                                                                                 >  5:309448/1‑149 (MQ=255)    ccAAAAAACACAGAGGCCCATCCTATGTTCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCTGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGAGGCCTGGGAGGATGCACCTGGGGt                                                                                                                                                <  6:26096/149‑1 (MQ=255)    ccAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGt                                                                                                                                                <  2:156577/149‑1 (MQ=255)     cAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTg                                                                                                                                               >  4:133246/1‑149 (MQ=255)      aaaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGc                                                                                                                                              <  1:299176/149‑1 (MQ=255)       aaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCa                                                                                                                                             >  5:5161/1‑149 (MQ=255)       aaaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCa                                                                                                                                             >  2:198404/1‑149 (MQ=255)        aaaaCACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCAt                                                                                                                                            >  8:117125/1‑149 (MQ=255)           acacAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGcc                                                                                                                                         <  5:245970/149‑1 (MQ=255)            cacaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCg                                                                                                                                        <  7:91979/149‑1 (MQ=255)             acaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACTCCAGCCCCTGGCGGCCTGGGAGGATGCACCCCAGGTGCAtcctc                                                                                                                                        <  1:69532/149‑5 (MQ=255)             acaGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCgg                                                                                                                                       <  7:186428/149‑1 (MQ=255)               agagGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGgc                                                                                                                                     <  3:133246/149‑1 (MQ=255)                    cccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAc                                                                                                                                <  1:198404/149‑1 (MQ=255)                     ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGTCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                               >  3:229837/1‑149 (MQ=255)                     ccAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGAcc                                                                                                                               <  7:117125/149‑1 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGaggcaggc                                                                                                                    >  6:141425/1‑148 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGaggcaggc                                                                                                                    <  5:141425/148‑1 (MQ=255)                                 ccTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCg                                                                                                                   >  8:174360/1‑149 (MQ=255)                                  cTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCgg                                                                                                                  >  6:335703/1‑149 (MQ=255)                                        gTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACct                                                                                                             <  2:99831/148‑1 (MQ=255)                                              gCCTGGCGGTCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggctgg                                                                                                       <  4:229837/148‑1 (MQ=255)                                              gCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggctgg                                                                                                       <  5:335703/148‑1 (MQ=255)                                               ccTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCggg                                                                                                                                      <  7:291791/116‑1 (MQ=255)                                               ccTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCggg                                                                                                                                      >  8:291791/1‑116 (MQ=255)                                                cTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                    >  5:100417/1‑149 (MQ=255)                                                cTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCg                                                                                                    >  6:189371/1‑149 (MQ=255)                                                        cGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGgcgc                                                                                                                                   <  4:191091/110‑1 (MQ=255)                                                        cGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGc                                                                                                     <  1:322319/140‑1 (MQ=255)                                                        cGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCca                                                                                            <  6:265552/149‑1 (MQ=255)                                                        cGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCAACTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGc                                                                                                     >  2:322319/1‑140 (MQ=255)                                                         ggggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCcac                                                                                           <  4:81291/149‑1 (MQ=255)                                                         ggggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCca                                                                                            >  1:23758/1‑148 (MQ=255)                                                           ggCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAg                                                                                         <  3:329660/149‑1 (MQ=255)                                                                    tGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggctgg                                                                                                       <  1:151680/126‑1 (MQ=255)                                                                    tGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACctggctgg                                                                                                       >  2:151680/1‑126 (MQ=255)                                                                    tGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCttt                                                                                <  8:118514/149‑1 (MQ=255)                                                                    tGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCttt                                                                                >  5:18129/1‑149 (MQ=255)                                                                      ccGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACct                                                                                                             >  1:307763/1‑118 (MQ=255)                                                                      ccGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACct                                                                                                             <  2:307763/118‑1 (MQ=255)                                                                            aTCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGGACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCt                                                                        <  7:174360/149‑1 (MQ=255)                                                                            aTCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGTGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCt                                                                        >  5:225750/1‑149 (MQ=255)                                                                               gTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGa                                                                      <  2:23758/148‑1 (MQ=255)                                                                                    ggcCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACa                                                                >  6:255392/1‑149 (MQ=255)                                                                                     gcCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACAt                                                               >  6:341658/1‑149 (MQ=255)                                                                                        aaCCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCcg                                                            >  2:353786/1‑149 (MQ=255)                                                                                            gcgACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCa                                                                   >  4:114706/1‑138 (MQ=255)                                                                                            gcgACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCa                                                                   <  3:114706/138‑1 (MQ=255)                                                                                               acgaATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCg                                                     <  6:100417/149‑1 (MQ=255)                                                                                                cgaATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCgg                                                    <  6:18129/149‑1 (MQ=255)                                                                                                 gaATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGTGGTACCTGGCTGGGCTTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCggg                                                   <  6:225750/149‑1 (MQ=255)                                                                                                    ttCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAg                                                                                         >  7:341025/1‑108 (MQ=255)                                                                                                    ttCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAg                                                                                         <  8:341025/108‑1 (MQ=255)                                                                                                       tATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGgc                                             >  8:284227/1‑149 (MQ=255)                                                                                                         tGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTggg                                                                                                                                                  <  3:301895/46‑1 (MQ=255)                                                                                                         tGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTggg                                                                                                                                                  >  4:301895/1‑46 (MQ=255)                                                                                                           cccgcccCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAg                                         >  4:1464/1‑149 (MQ=255)                                                                                                            ccgcccCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                        <  6:309448/149‑1 (MQ=255)                                                                                                            ccgcccCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAgg                                        >  2:131531/1‑149 (MQ=255)                                                                                                               ccccACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGc                                     <  5:255392/149‑1 (MQ=255)                                                                                                                 ccACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGc                                     >  4:252412/1‑147 (MQ=255)                                                                                                                 ccACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGc                                     <  3:252412/147‑1 (MQ=255)                                                                                                                 ccACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCct                                   <  6:5161/149‑1 (MQ=255)                                                                                                                  cACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctc                                  <  7:284227/149‑1 (MQ=255)                                                                                                                        aGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctcgc                                >  4:337100/1‑146 (MQ=255)                                                                                                                        aGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCctcgc                                <  3:337100/146‑1 (MQ=255)                                                                                                                          ccccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGCCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCACACAGGGCCGCTTTGCCGCCCTGACCAACATCCCCGACCCGGGCCACTCGCAGGGGCCTCGCTcgcg                           >  4:203743/1‑149 (MQ=255)                                                                                                                           cccTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCgg                          >  2:133345/1‑149 (MQ=255)                                                                                                                                gcggcCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGt                       >  6:346106/1‑147 (MQ=255)                                                                                                                                     cTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtgg                >  5:249388/1‑149 (MQ=255)                                                                                                                                      tGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCgg                          >  2:168998/1‑138 (MQ=255)                                                                                                                                      tGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCgg                          <  1:168998/138‑1 (MQ=255)                                                                                                                                      gcccggcaTGCACCCCAGGTGCATCCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGAtggt               >  4:243426/8‑149 (MQ=255)                                                                                                                                        ccggcaTGCACCCCAGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCggg                         <  4:159428/132‑1 (MQ=255)                                                                                                                                        ccggcaTGCACCCCAGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCggg                         >  3:159428/6‑137 (MQ=255)                                                                                                                                         gaggaTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCtggtgg             >  6:69954/1‑148 (MQ=255)                                                                                                                                             aTGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCg        >  5:156798/1‑149 (MQ=255)                                                                                                                                                 aCCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGCCCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTc    >  6:248100/1‑149 (MQ=255)                                                                                                                                                  ccTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCt   >  8:298799/1‑149 (MQ=255)                                                                                                                                                   cTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg  >  3:163459/1‑149 (MQ=255)                                                                                                                                                   cTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTg  >  1:281962/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TGACCAAAAAACACAGAGGCCCAACCTATGTGCCTGATCGTATTTGCCTGGCGGCCGGGGCACGCCCTGCCGCTGATCGTGGCGGCCAACCGCGACGAATTCTATGCCCGCCCCACT‑CAGCCCCTGGCGGCCTGGGAGGATGCACCTGGGGTGCATGCCGGGCGCGACCTGGAGGCAGGCGGTACCTGGCTGGGCGTAGGCCCACAGGGGCGCTTTGCCGCGCTGACCAACATCCGCGACCCGGGCCAGGCGCAGGGGCCTCGCTCGCGGGGTGAGCTGGTGGCAGCGTTTCTG  >  NC_002947/5869466‑5869759

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