Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,798,278	G→T	5.7%	R529S (CGC→AGC)	PP_5077 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,798,278	0	G	T	5.7%	284.1 / 5.3	88	R529S (CGC→AGC)	PP_5077	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base T (2/3);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTG  >  NC_002947/5798138‑5798423                                                                                                                                              |                                                                                                                                                  aGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           <  4:222547/149‑1 (MQ=255)    cATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGc                                                                                                                                          >  5:21046/1‑147 (MQ=255)      tGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGc                                                                                                                                      >  3:81282/1‑149 (MQ=255)            cAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCTGGCCACGGCCCGctgata                                                                                                                                 <  3:61018/148‑5 (MQ=255)            cAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgt                                                                                                                                >  6:22678/1‑149 (MQ=255)             aaCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgtt                                                                                                                               <  4:81282/149‑1 (MQ=255)             aaCTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCttgtt                                                                                                                               <  3:288414/149‑1 (MQ=255)                    ggCTTGTGGACGGCAGTCCGGGACAGAAAAGCCTGCCGCTCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTTTTGGACGCt                                                                                                                        <  5:279071/149‑1 (MQ=255)                      cTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGc                                                                                                                      <  5:22678/149‑1 (MQ=255)                          tgGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGAc                                                                                                                  >  1:133590/1‑149 (MQ=255)                            gCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACg                                                                                                                 <  8:105872/148‑1 (MQ=255)                             ccGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGt                                                                                                               <  2:133590/149‑1 (MQ=255)                              cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                   <  8:63162/128‑3 (MQ=255)                              cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                   >  5:104106/1‑126 (MQ=255)                              cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                   <  6:104106/128‑3 (MQ=255)                              cGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCTGGCCACGGCCCGctga                                                                                                                                   >  7:63162/1‑126 (MQ=255)                                gCAGTCCGGGACAGAAAATCCTTCCGCGCTCGGGCGCAGCAGGCTATTTGAAATCTATAATTCTCTTATCTACTGGGCGGCGGGGTTATGCCGCGCTGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           <  7:50944/119‑1 (MQ=255)                                gCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           >  8:50944/1‑118 (MQ=255)                                gCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           >  4:298459/1‑118 (MQ=255)                                gCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           <  3:298459/118‑1 (MQ=255)                                gCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGc                                                                                                                                          <  5:195275/119‑1 (MQ=255)                                gCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGc                                                                                                                                          >  6:195275/1‑119 (MQ=255)                                  aGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGATGGACGCTGCCGACGGTACGTg                                                                                                          >  2:112242/1‑149 (MQ=255)                                            aGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGt                                                                                                  >  7:321836/1‑147 (MQ=255)                                            aGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGt                                                                                                  <  8:321836/147‑1 (MQ=255)                                              aaaaGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                              >  2:202414/1‑149 (MQ=255)                                              aaaaGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                              >  6:14150/1‑149 (MQ=255)                                                    cTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCt                                                                                        >  7:320042/1‑149 (MQ=255)                                                    cTGCCGCGCTCGCGCGCATCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCt                                                                                        <  6:224120/149‑1 (MQ=255)                                                      gCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTgg                                                                                      >  3:13273/1‑149 (MQ=255)                                                       ccGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCGAGAATTCGCTTAGCTACTGGGCGGCGGGGTGCTGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGa                                                                                     >  1:42801/1‑149 (MQ=255)                                                        cgcgCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGAc                                                                                    <  2:42603/149‑1 (MQ=255)                                                         gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           <  3:310635/93‑1 (MQ=255)                                                         gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTgg                                                                                                                                           >  4:310635/1‑93 (MQ=255)                                                         gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGAcc                                                                                   >  7:245320/1‑149 (MQ=255)                                                         gcgcTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGAcc                                                                                   >  8:268886/1‑149 (MQ=255)                                                             tCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCg                                                                               >  6:169264/1‑149 (MQ=255)                                                               gggCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCg                                                                               <  2:309403/147‑1 (MQ=255)                                                                ggCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTg                                                                            >  4:142237/1‑149 (MQ=255)                                                                 gcgcAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTa                                                                                                              >  2:135166/1‑114 (MQ=255)                                                                 gcgcAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTa                                                                                                              <  1:135166/114‑1 (MQ=255)                                                                     agcagGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAg                                                                       <  6:21046/149‑1 (MQ=255)                                                                     agcagGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAg                                                                       <  1:202414/149‑1 (MQ=255)                                                                              gTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGa                                                                <  8:245320/147‑1 (MQ=255)                                                                               tCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTg                                                             <  5:14150/149‑1 (MQ=255)                                                                                 gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                              <  5:128393/114‑1 (MQ=255)                                                                                 gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTAc                                                                                              >  6:128393/1‑114 (MQ=255)                                                                                 gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGcc                                                           >  8:117943/1‑149 (MQ=255)                                                                                 gAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGcc                                                           <  2:25482/149‑1 (MQ=255)                                                                                   aaGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGcc                                                           <  4:228524/147‑1 (MQ=255)                                                                                       tAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAAcc                                                     >  2:238861/1‑149 (MQ=255)                                                                                       tAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAAcc                                                     >  4:315651/1‑149 (MQ=255)                                                                                         gactccGCTTAGCTACTGGGCGGCGGGGGGATGCCGC‑CCGGAAGGTCTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGGGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGc                                                   <  2:306161/144‑1 (MQ=255)                                                                                         gAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCg                                                    <  4:13273/148‑1 (MQ=255)                                                                                            ttCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                <  5:169264/149‑1 (MQ=255)                                                                                            ttCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                <  2:42801/149‑1 (MQ=255)                                                                                               gCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTc                                             >  6:134488/1‑149 (MQ=255)                                                                                                 ttAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTcgc                                           >  7:153849/1‑149 (MQ=255)                                                                                                     cTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                       >  5:48824/1‑149 (MQ=255)                                                                                                     cTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                       >  2:315754/1‑149 (MQ=255)                                                                                                     cTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGtt                                       >  8:229390/1‑149 (MQ=255)                                                                                                      tACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCttta                                                                <  5:99948/123‑3 (MQ=255)                                                                                                      tACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGa                                                                >  6:99948/1‑123 (MQ=255)                                                                                                        cTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTg                                      <  6:48824/147‑1 (MQ=255)                                                                                                        cTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGc                                    >  4:192484/1‑149 (MQ=255)                                                                                                         tGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGc                                                   <  1:11965/133‑1 (MQ=255)                                                                                                         tGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGc                                                   >  2:11965/1‑133 (MQ=255)                                                                                                           ggcggcGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTg                                      >  8:132558/1‑144 (MQ=255)                                                                                                           ggcggcGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTg                                      <  7:132558/144‑1 (MQ=255)                                                                                                           ggcggcGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGa                                  >  2:334192/1‑148 (MQ=255)                                                                                                             cggcggGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAAc                               >  2:231873/1‑149 (MQ=255)                                                                                                                   ggTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                >  4:168708/1‑126 (MQ=255)                                                                                                                   ggTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGc                                                <  3:168708/126‑1 (MQ=255)                                                                                                                     tGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                       <  5:134488/149‑1 (MQ=255)                                                                                                                          ccgc‑cCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                  >  6:50179/1‑149 (MQ=255)                                                                                                                          ccgc‑cCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCac                  >  8:120520/1‑149 (MQ=255)                                                                                                                          ccgc‑cCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCa                   <  7:229390/148‑1 (MQ=255)                                                                                                                            gc‑cCGGAATGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTTCCGTCGGTACGTGGCCACGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCacga                <  1:112242/149‑1 (MQ=255)                                                                                                                            gc‑cCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                       <  6:170302/142‑1 (MQ=255)                                                                                                                            gc‑cCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTaa                       >  5:170302/1‑142 (MQ=255)                                                                                                                            gc‑cCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCacga                <  3:192484/149‑1 (MQ=255)                                                                                                                                  gAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTAc           <  7:268886/149‑1 (MQ=255)                                                                                                                                    aGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTa            >  7:295239/1‑146 (MQ=255)                                                                                                                                     ggTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACaga        <  7:117943/149‑1 (MQ=255)                                                                                                                                     ggTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACaga        <  1:231873/149‑1 (MQ=255)                                                                                                                                        aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCg                          >  5:113399/1‑128 (MQ=255)                                                                                                                                        aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCg                          <  6:113399/128‑1 (MQ=255)                                                                                                                                        aTCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGc     >  8:108803/1‑149 (MQ=255)                                                                                                                                              gggCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg  <  1:66215/146‑1 (MQ=255)                                                                                                                                              gggCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTg  >  2:66215/1‑146 (MQ=255)                                                                                                                                              |                                                                                                                                                  AGCCATGACTGCAACTAAGGGCTTGTGGCCGGCAGTCCGGGACAGAAAAGCCTGCCGCGCTCGGGCGCAGCAGGCTAGTCGAAAGCTAGAATTCGCTTAGCTACTGGGCGGCGGGGTGATGCCGC‑CCGGAAGGTATCAGCGGGCCGTGGCCAGCTCTTGTTGGACGCTGCCGACGGTACGTGGCCAAGGTTTACCAGCCTGGACCTTCGCTGGCAAGTTCTTGATTGCCGCAACCGCTGCGTCGCGGTTGGCGAAACTGCCGTAAGTCACGACGTACAGAGGCTTG  >  NC_002947/5798138‑5798423

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