Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,761,260	C→T	6.1%	G136G (GGC→GGT)	PP_1567 →	capsid protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,761,260	0	C	T	6.1%	197.3 / 2.2	65	G136G (GGC→GGT)	PP_1567	capsid protein
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (2/2);  total (32/33)
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.30e-01

CTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCA  >  NC_002947/1761116‑1761406                                                                                                                                                 |                                                                                                                                                   ctcCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACcaccggcac                                                                                                                                                  <  1:61970/147‑1 (MQ=255)    cGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTg                                                                                                                                             >  3:272684/1‑149 (MQ=255)    cGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTg                                                                                                                                             <  8:41396/149‑1 (MQ=255)    cGCCAATGATCCAGTTGACTTCAACTCGGCTGGCGAAGAACTGACCGAGGGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTg                                                                                                                                             <  8:267529/149‑1 (MQ=255)          tGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCg                                                                                                                                       <  7:19443/149‑1 (MQ=255)           gATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtaga                                                                                                                                        <  4:18728/147‑5 (MQ=255)           gATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtaga                                                                                                                                        >  3:18728/1‑143 (MQ=255)           gATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCgg                                                                                                                                      >  3:242220/1‑149 (MQ=255)            aTCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGc                                                                                                                                     <  5:1932/149‑1 (MQ=255)            aTCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGc                                                                                                                                     <  1:220614/149‑1 (MQ=255)             tCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGc                                                                                                                                     <  5:7016/148‑1 (MQ=255)               cAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTa                                                                                                                                  <  3:79586/149‑1 (MQ=255)                aGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAg                                                                                                                                 <  3:16886/149‑1 (MQ=255)                 gTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGc                                                                                                                                >  2:188118/1‑149 (MQ=255)                          aaGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTc                                                                                                                                        >  6:186538/1‑132 (MQ=255)                          aaGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTc                                                                                                                                        <  5:186538/132‑1 (MQ=255)                             tCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGAc                                                                                                                    <  6:11210/149‑1 (MQ=255)                                      gaagaaCTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTcc                                                                                                           >  6:187857/1‑149 (MQ=255)                                       aagaaCTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTccc                                                                                                          >  7:167823/1‑149 (MQ=255)                                        agaaCTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCg                                                                                                         <  4:272684/149‑1 (MQ=255)                                               aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtaga                                                                                                                                        <  2:124251/111‑5 (MQ=255)                                               aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGTGCCGGTggtaga                                                                                                                                        >  1:124251/1‑107 (MQ=255)                                               aCCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGt                                                                                                  <  1:188118/149‑1 (MQ=255)                                                ccGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAATAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcgc                                                                                                                >  6:303576/1‑134 (MQ=255)                                                ccGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAATAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcgc                                                                                                                <  5:303576/134‑1 (MQ=255)                                                 cGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTgg                                                                                                >  4:228190/1‑149 (MQ=255)                                                  gAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTggg                                                                                               >  8:139473/1‑148 (MQ=255)                                                   aGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTa                                                                                                                                  <  2:351207/113‑1 (MQ=255)                                                   aGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTa                                                                                                                                  >  1:351207/1‑113 (MQ=255)                                                           gacTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGa                                                                                      <  5:187857/149‑1 (MQ=255)                                                             cTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTTGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGtt                                                                                              <  6:322612/139‑1 (MQ=255)                                                             cTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTTGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGtt                                                                                              >  5:322612/1‑139 (MQ=255)                                                              ttCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAgcg                                                                                   >  7:226260/1‑149 (MQ=255)                                                                cAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAgcgc                                                                                  <  4:242220/148‑1 (MQ=255)                                                                 aaGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcgc                                                                                                                <  7:83238/117‑1 (MQ=255)                                                                 aaGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCcgc                                                                                                                >  8:83238/1‑117 (MQ=255)                                                                      gCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTc                                                                            <  8:167823/148‑1 (MQ=255)                                                                      gCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTc                                                                            >  3:291128/1‑148 (MQ=255)                                                                       cTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCAc                                                                          >  8:172800/1‑149 (MQ=255)                                                                       cTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCAc                                                                          <  8:320716/149‑1 (MQ=255)                                                                         gCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCgg                                                                                                                                      >  6:127936/1‑87 (MQ=255)                                                                          cAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCa                                                                        <  5:35091/148‑1 (MQ=255)                                                                               tcgcggtcgcggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTcgc                                                                   >  7:313557/1‑148 (MQ=255)                                                                                 gcggtcgcggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGAc                                                                >  4:221844/1‑149 (MQ=255)                                                                                       gcggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTcg                                                                    <  6:215402/139‑1 (MQ=255)                                                                                       gcggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCCTTCCCGGCATCGTGGTTTGGCCGGAGCGCGAGTTCACCATTcg                                                                    >  5:215402/1‑139 (MQ=255)                                                                                         ggCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCAt                                                        >  1:196300/1‑149 (MQ=255)                                                                                            atcgatcgAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACGGGTGGCGTAGGAGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGc                                                      >  5:251481/1‑148 (MQ=255)                                                                                            aTCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGCCCTGATCATGcc                                                     >  7:70447/1‑149 (MQ=255)                                                                                               gTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAg                                                   <  2:196300/148‑1 (MQ=255)                                                                                                   gAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCg                                               >  2:256695/1‑148 (MQ=255)                                                                                                     gTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCCGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTAc                                            <  1:256695/149‑1 (MQ=255)                                                                                                      tAGGCCGGAAGACCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTAc                                            <  5:38652/148‑1 (MQ=255)                                                                                                      tAGGCCGGAAGACCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTAc                                            >  6:38652/1‑148 (MQ=255)                                                                                                         gCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGc                                        >  7:16620/1‑149 (MQ=255)                                                                                                              aaGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGaa                                   <  3:228190/149‑1 (MQ=255)                                                                                                               aGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAAc                                  <  4:291128/149‑1 (MQ=255)                                                                                                                        ccAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAg                         <  8:313557/149‑1 (MQ=255)                                                                                                                            cATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCg                     >  5:77836/1‑149 (MQ=255)                                                                                                                             aTTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGt                    <  8:226260/149‑1 (MQ=255)                                                                                                                              ttACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAg                              <  3:273752/138‑1 (MQ=255)                                                                                                                              ttACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAg                              >  4:273752/1‑138 (MQ=255)                                                                                                                                 ccAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAg                >  1:221971/1‑149 (MQ=255)                                                                                                                                     ctctACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGgaga            >  6:26016/1‑149 (MQ=255)                                                                                                                                          ccaccGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATc                                                                                                    >  5:237222/1‑56 (MQ=255)                                                                                                                                          ccaccGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATc                                                                                                    <  6:237222/56‑1 (MQ=255)                                                                                                                                           caccggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGAc           <  6:148501/144‑1 (MQ=255)                                                                                                                                           caccggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGAc           >  5:148501/1‑144 (MQ=255)                                                                                                                                              cggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACg                                 >  8:121349/1‑119 (MQ=255)                                                                                                                                              cggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTc    >  4:283937/1‑148 (MQ=255)                                                                                                                                               ggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCa  <  2:221971/149‑1 (MQ=255)                                                                                                                                               ggcaccggTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCa  <  7:139473/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CTCCGCCAATGATCCAGTTGACTTCAAGTCGGCTGGCGAAGAACTGACCGAGTGTGACGACTTCAAGTCGCTGCAAACTCGCGGTCGCGGCATCGTTCGAGTAGGCCGGAAGGCCGTCACCAACATTACCAGCTCTACCACCGGCACCGGTGGCGTCGGCGTAGCCATCCAGCCGACCCGCGTTCCCGGCATCGTGGTTGGGCCGGAGCGCGAGTTCACCATTCGCGACCTGATCATGCCAGGCCGTACTGGCTCGAACGCAGTCGAGTTCGTGCAGGAGACCGGCTTCCA  >  NC_002947/1761116‑1761406

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