Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,063,002	T→A	5.1%	E446V (GAG→GTG) ‡	PP_3581 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,063,002	0	T	A	5.1%	252.7 / 2.6	78	E446V (GAG→GTG) ‡	PP_3581	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/41);  new base A (2/2);  total (35/43)
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.34e-01

CCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTC  >  NC_002947/4062855‑4063150                                                                                                                                                    |                                                                                                                                                     ccGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTc                                                                                                                                                     <  6:332738/149‑1 (MQ=255)  cGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTc                                                                                                                                                     >  5:307714/1‑148 (MQ=255)     tCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgc                                                                                                                                                 <  7:270459/149‑1 (MQ=255)       ttttCCGCCAGGCCTATGAATCGGTCAGGCAGCAGGAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTc                                                                                                                                                     <  3:246744/143‑1 (MQ=255)       ttttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAgcgc                                                                                                                                               >  6:256351/1‑149 (MQ=255)       ttttCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAgcgc                                                                                                                                               >  2:370557/1‑149 (MQ=255)       ttttCCGCCAGGCCGATGAAGCGCTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTCGGCTGCTAATGAAAGCGGAACTGGCCCTGATc                                                                                                                                                     >  4:246744/1‑143 (MQ=255)          tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCa                                                                                                                                             <  5:256351/148‑1 (MQ=255)          tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCa                                                                                                                                             <  1:342056/148‑1 (MQ=255)          tCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCa                                                                                                                                             <  1:370557/148‑1 (MQ=255)            cGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAgg                                                                                                                                           <  8:346307/148‑1 (MQ=255)             gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgttc                                                                                                                                         >  8:324989/1‑146 (MQ=255)             gCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGgcg                                                                                                                                         >  1:114601/1‑149 (MQ=255)                 ggCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGGAGCAGCGCCAgttccgct                                                                                                                                     <  7:324989/149‑8 (MQ=255)                  gCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGt                                                                                                                                    >  4:179531/1‑149 (MQ=255)                     gATGAAGCGGTCAGGCAGCAGCTGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCt                                                                                                                                 <  4:242971/149‑1 (MQ=255)                          aGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgcag                                                                                                                             >  6:125783/1‑148 (MQ=255)                           gCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTgcagca                                                                                                                           <  5:125783/149‑1 (MQ=255)                             ggTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCagcagc                                                                                                                         <  2:114601/149‑1 (MQ=255)                                        agcagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAActgctg                                                                                                              >  8:214855/1‑149 (MQ=255)                                          cagGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGgc                                                                                                            <  4:75074/149‑1 (MQ=255)                                               ccTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCa                                                                                                       >  6:70704/1‑149 (MQ=255)                                                ctgctgCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTCGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAt                                                                                                      <  3:179531/149‑1 (MQ=255)                                                   ctgcGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCAtt                                                                                                     <  5:151274/147‑1 (MQ=255)                                                      cGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                  <  6:349916/147‑1 (MQ=255)                                                      cGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTg                                                                                                 <  2:331400/148‑1 (MQ=255)                                                      cGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGt                                                                                                <  7:214855/149‑1 (MQ=255)                                                       gggTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGtt                                                                                               >  1:330941/1‑149 (MQ=255)                                                         gTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCa                                                                                             >  3:167968/1‑149 (MQ=255)                                                          tGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGTGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCAt                                                                                            >  7:99043/1‑149 (MQ=255)                                                            cTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                           <  6:307714/148‑1 (MQ=255)                                                             tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTgg                                                                                                             <  1:265532/129‑1 (MQ=255)                                                             tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTgg                                                                                                             >  2:265532/1‑129 (MQ=255)                                                             tCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGa                                                                                         <  7:254226/149‑1 (MQ=255)                                                               ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                  <  7:110360/138‑1 (MQ=255)                                                               ccAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCt                                                                                                  >  8:110360/1‑138 (MQ=255)                                                                cAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCAt                                                                                            >  7:20143/1‑143 (MQ=255)                                                                cAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCAt                                                                                            <  8:20143/143‑1 (MQ=255)                                                                cAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGacac                                                                                      <  3:25603/149‑1 (MQ=255)                                                                cAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGaca                                                                                       >  2:203038/1‑148 (MQ=255)                                                                 aGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATc                                                                                           <  7:293712/143‑1 (MQ=255)                                                                 aGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGCCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGCCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCCCGCATTGCTCTTCATc                                                                                           >  8:293712/1‑143 (MQ=255)                                                                    gcagcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCg                                                                                  <  3:349570/149‑1 (MQ=255)                                                                      agcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                   >  2:100353/1‑98 (MQ=255)                                                                      agcaACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTc                                                                                                                                   <  1:100353/98‑1 (MQ=255)                                                                           cGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGc                                                                           <  1:139164/149‑1 (MQ=255)                                                                            gCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaa                                                                                >  2:241228/1‑143 (MQ=255)                                                                            gCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaa                                                                                <  1:241228/143‑1 (MQ=255)                                                                              ttCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTg                                                                         >  3:30015/1‑148 (MQ=255)                                                                              ttCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTg                                                                         >  8:141887/1‑148 (MQ=255)                                                                                   cACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCg                                                                   <  8:46552/149‑1 (MQ=255)                                                                                      cGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTc                                                                <  4:30015/149‑1 (MQ=255)                                                                                            gTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGTGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                          >  4:82238/1‑149 (MQ=255)                                                                                               ggCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGTGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGgc                                                        <  8:99043/148‑1 (MQ=255)                                                                                               ggCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                            >  8:288818/1‑128 (MQ=255)                                                                                               ggCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                            <  7:288818/128‑1 (MQ=255)                                                                                                tcAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTgg                                                                                                             <  8:9163/93‑1 (MQ=255)                                                                                                gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCtt                                                              <  3:71281/141‑1 (MQ=255)                                                                                                gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCtt                                                              >  4:71281/1‑141 (MQ=255)                                                                                                gCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTACTGCAGCAGCTGCAACTGCTgg                                                                                                             >  7:9163/1‑94 (MQ=255)                                                                                                    gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaagaa                                                                             >  6:10932/1‑122 (MQ=255)                                                                                                    gCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCgaagaa                                                                             <  5:10932/122‑1 (MQ=255)                                                                                                     ccTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                 <  1:203038/149‑1 (MQ=255)                                                                                                     ccTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTg                                                 <  4:167968/149‑1 (MQ=255)                                                                                                               gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                          <  4:342525/130‑1 (MQ=255)                                                                                                               gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                          >  3:342525/1‑130 (MQ=255)                                                                                                               gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                          >  3:359191/1‑130 (MQ=255)                                                                                                               gAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCg                                                          <  4:359191/130‑1 (MQ=255)                                                                                                                     gggCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCa                                  <  7:141887/148‑1 (MQ=255)                                                                                                                        cTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                            >  8:187349/1‑103 (MQ=255)                                                                                                                        cTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGaagaag                                                                            <  7:187349/103‑1 (MQ=255)                                                                                                                         tGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTc                             >  2:133090/1‑149 (MQ=255)                                                                                                                           gTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCgg                           >  3:142712/1‑149 (MQ=255)                                                                                                                                  aaGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                    <  4:142712/149‑1 (MQ=255)                                                                                                                                  aaGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                    >  6:373499/1‑149 (MQ=255)                                                                                                                                    gCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACg                    >  5:88859/1‑147 (MQ=255)                                                                                                                                       gAACTGGTGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCATCGGCGTTACgcagca               <  3:82238/149‑1 (MQ=255)                                                                                                                                          cTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCt            >  5:151237/1‑149 (MQ=255)                                                                                                                                               gctgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGc       <  5:373499/149‑1 (MQ=255)                                                                                                                                               gctgctCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGc       >  1:90177/1‑149 (MQ=255)                                                                                                                                                    tCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTc  >  8:289011/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCGGTCTTTTCCGCCAGGCCGATGAAGCGGTCAGGCAGCAGCAGGCCCTGCTGCGGGTGCTCCCAGCGCAGCAACGCTTCGGCACCGATCGGTTGGCAGGCCTGGGCGTCGAACTTGGGCTGGTAATGCAAGCGGAACTGGCGCTGCTCCAGCGCCAGGCGCAAGTCCTGCAGCAGCTGCAACTGCTGGCGCGCATTGCTGTTCATCGACACATCGAAGAAGCTGTAACCGTTCTTGCCGGCGCTTTTGGCATGGTACATGGCAGCGTCGGCGTTACGCAGCAGCTCGTGCTGGTC  >  NC_002947/4062855‑4063150

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