Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,761,927	A→T	5.7%	N359Y (AAC→TAC) ‡	PP_1567 →	capsid protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,761,927	0	A	T	5.7%	220.0 / 2.6	70	N359Y (AAC→TAC) ‡	PP_1567	capsid protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base T (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GACTGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCAGAGTCGTTCG  >  NC_002947/1761780‑1762068                                                                                                                                                    |                                                                                                                                               gACTGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGa                                                                                                                                                <  3:101264/148‑1 (MQ=255)         ggACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCgg                                                                                                                                       >  4:176858/1‑149 (MQ=255)           aCATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGcc                                                                                                                                     >  8:58995/1‑149 (MQ=255)             aTCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGcc                                                                                                                                     >  4:213233/1‑147 (MQ=255)             aTCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGcc                                                                                                                                     >  3:233623/1‑147 (MQ=255)                gAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGgt                                                                                                                                >  2:33320/1‑149 (MQ=255)                      aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCg                                                                                                                           <  5:128099/148‑1 (MQ=255)                      aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGa                                                                                                                          >  3:145807/1‑149 (MQ=255)                      aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGa                                                                                                                          >  7:321223/1‑149 (MQ=255)                      aCCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGa                                                                                                                          <  4:233623/149‑1 (MQ=255)                        cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGAcc                                                                                                                        >  4:154122/1‑149 (MQ=255)                        cAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGAcc                                                                                                                        <  4:145807/149‑1 (MQ=255)                         aaGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACcg                                                                                                                       >  8:223274/1‑149 (MQ=255)                         aaGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACcg                                                                                                                       >  8:223270/1‑149 (MQ=255)                           ggACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACcgcg                                                                                                                     <  8:332236/149‑1 (MQ=255)                            gACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCATcgcgc                                                                                                                                       <  1:22970/130‑2 (MQ=255)                            gACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCATcgcgc                                                                                                                                       >  2:22970/1‑129 (MQ=255)                            gACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGcc                                                                                                                                     >  8:33857/1‑132 (MQ=255)                            gACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGcc                                                                                                                                     <  7:33857/132‑1 (MQ=255)                                  aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc                                                                                                                            >  3:198853/1‑135 (MQ=255)                                  aaCGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc                                                                                                                            <  4:198853/135‑1 (MQ=255)                                    cGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCgg                                                                                                            <  8:321223/149‑1 (MQ=255)                                      gCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCggat                                                                                                          <  7:223270/149‑3 (MQ=255)                                      gCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCt                                                                                                          <  7:223274/149‑1 (MQ=255)                                          cTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTc                                                                                                       <  5:275569/148‑1 (MQ=255)                                           tACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGccc                                                                                                                                    <  6:100190/118‑1 (MQ=255)                                           tACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGccc                                                                                                                                    >  5:100190/1‑118 (MQ=255)                                            aCATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTc                                                                                                       <  4:8089/146‑1 (MQ=255)                                            aCATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTc                                                                                                       >  3:8089/1‑146 (MQ=255)                                             cATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGa                                                                                                     <  7:82238/147‑1 (MQ=255)                                             cATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                    <  4:336909/148‑1 (MQ=255)                                                  gggTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTcc                                                                                              >  8:13163/1‑149 (MQ=255)                                                        aTGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCTAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGaa                                                                                        <  7:58995/149‑1 (MQ=255)                                                          gTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAg                                                                                                                                  <  2:220372/105‑1 (MQ=255)                                                             cAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACgg                                                                                   >  8:283602/1‑149 (MQ=255)                                                             cAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACg                                                                                    <  3:176858/148‑1 (MQ=255)                                                              aGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGTTAGGCGCGATGAACATCGCGGCCCAg                                                                                                                                  <  8:287313/101‑1 (MQ=255)                                                              aGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGTTAGGCGCGATGAACATCGCGGCCCAg                                                                                                                                  >  7:287313/1‑101 (MQ=255)                                                                gAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                    <  6:87238/129‑1 (MQ=255)                                                                gAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg                                                                                                    >  5:87238/1‑129 (MQ=255)                                                                      ggcgTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc                                                                                                                            >  8:46087/1‑99 (MQ=255)                                                                      ggcgTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc                                                                                                                            <  7:46087/99‑1 (MQ=255)                                                                                  cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCAcct                                                                                          <  2:176419/121‑2 (MQ=255)                                                                                  cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACAGCGAGGAAGCGGCTGTCGAGGTCTCCACCg                                                                                          >  1:176419/1‑121 (MQ=255)                                                                                      ggAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGt                                                           >  4:341024/1‑148 (MQ=255)                                                                                       gAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa                                                         <  3:213233/149‑1 (MQ=255)                                                                                       gAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa                                                         >  4:103367/1‑149 (MQ=255)                                                                                            tGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATc                                                    >  3:274608/1‑149 (MQ=255)                                                                                             gCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACgg                                                                                   >  6:287698/1‑117 (MQ=255)                                                                                             gCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGACTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACgg                                                                                   <  5:287698/117‑1 (MQ=255)                                                                                                 gtggtgGACAGCAACGCGGTGCCAGAACGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCgcgc                                               >  1:205503/1‑149 (MQ=255)                                                                                                       gACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCgag                                           >  3:305586/1‑147 (MQ=255)                                                                                                              aCGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctg                                  <  7:13163/149‑1 (MQ=255)                                                                                                               cgcgGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctgg                                 <  3:341024/149‑1 (MQ=255)                                                                                                                    tGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctggcgct‑g                           <  7:283602/149‑1 (MQ=255)                                                                                                                        aGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctg                                  <  4:266371/139‑1 (MQ=255)                                                                                                                        aGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCctg                                  >  3:266371/1‑139 (MQ=255)                                                                                                                         gAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCga                                            >  4:92941/1‑128 (MQ=255)                                                                                                                               gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGcc                                    <  5:171074/130‑1 (MQ=255)                                                                                                                               gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGcc                                    >  6:171074/1‑130 (MQ=255)                                                                                                                               gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCg                >  2:160614/1‑149 (MQ=255)                                                                                                                                   tCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGcc              >  3:341414/1‑147 (MQ=255)                                                                                                                                     cTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCaga          <  3:103367/149‑1 (MQ=255)                                                                                                                                     cTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCaga          <  4:305586/149‑1 (MQ=255)                                                                                                                                     cTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCa            >  1:52209/1‑147 (MQ=255)                                                                                                                                        gTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCagag         <  4:274608/147‑1 (MQ=255)                                                                                                                                         tAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGGCGGTAGACCGccagcgcca       <  3:154122/149‑8 (MQ=255)                                                                                                                                            gcgcGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGATGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCCCCTGGCGCT‑GGCGGTCTACCGCCCAGAGtcgttc   <  2:205503/149‑1 (MQ=255)                                                                                                                                             cgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCct                                   <  8:235451/117‑1 (MQ=255)                                                                                                                                             cgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCct                                   >  7:235451/1‑117 (MQ=255)                                                                                                                                             cgcgATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCAGAGtcgttcg  <  1:33320/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GACTGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCT‑GGCGGTCTACCGCCCAGAGTCGTTCG  >  NC_002947/1761780‑1762068

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