Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,659,576	Δ1 bp	6.9%	coding (760/1536 nt)	PP_2331 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,659,576	0	A	.	6.9%	179.5 / 6.8	58	coding (760/1536 nt)	PP_2331	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/24);  new base . (2/2);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGAGTATTCGAGGAAGTCGGCGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCTTTTCCATCGGCA  >  NC_002947/2659429‑2659722                                                                                                                                                    |                                                                                                                                                    gCGAGTATTCGAGGAAGTCGGCGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTccacca                                                                                                                                                     >  2:67286/1‑148 (MQ=255)           gAGGAAGTCGGCGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGt                                                                                                                                          >  3:298108/1‑149 (MQ=255)                  tCGGCGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGAAGGCGGTCGCTCgg                                                                                                                                   >  5:151128/1‑149 (MQ=255)                     gCGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAg                                                                                                                                >  7:300917/1‑149 (MQ=255)                      cGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAgg                                                                                                                               >  6:373171/1‑149 (MQ=255)                      cGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAgg                                                                                                                               <  1:67286/149‑1 (MQ=255)                       gTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGc                                                                                                                              <  6:150424/149‑1 (MQ=255)                       gTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGc                                                                                                                              <  7:215282/149‑1 (MQ=255)                            tgtTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGc                                                                                                                         <  8:82874/149‑1 (MQ=255)                               tctcGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACC‑CAGCAGGCGGTCGCTCGGCAGGCCTTGCTCAc                                                                                                                     >  5:368778/1‑149 (MQ=255)                               tctcGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACC‑CAGCAGGCGGTCGCTCGGCAGGCCTTGCTCAc                                                                                                                     >  5:363661/1‑149 (MQ=255)                                  cgtcgGTCAGCTTGGCCAGGCGGATGGCGTTCATCCCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCg                                                                                                                   <  5:338055/149‑1 (MQ=255)                                   gtcgGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCg                                                                                                                   <  4:338668/148‑1 (MQ=255)                                       gTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTAGGCAGGCCTTGCTCACCGCTGTc                                                                                                              >  8:203303/1‑149 (MQ=255)                                               ggCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGtt                                                                                                        <  5:373171/147‑1 (MQ=255)                                                 ccAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAg                                                                                                    >  2:39224/1‑149 (MQ=255)                                                  cAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAg                                                                                                    >  2:272152/1‑148 (MQ=255)                                                  cAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGt                                                                                                   <  8:300917/149‑1 (MQ=255)                                                          tGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAg                                                                                           >  7:287068/1‑149 (MQ=255)                                                          tGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGAAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCa                                                                                            <  6:151128/148‑1 (MQ=255)                                                                  aTCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGcc                                                                                   >  1:341302/1‑149 (MQ=255)                                                                      cGCTGTTGTTCATCCTGAAGGTGACGCTGGCCTTCATGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTg                                                                               <  4:298108/149‑1 (MQ=255)                                                                         tgttgtTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAg                                                                            >  7:126932/1‑149 (MQ=255)                                                                         tgttgtTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAg                                                                            >  5:225974/1‑149 (MQ=255)                                                                         tgttgtTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAg                                                                            >  5:191949/1‑149 (MQ=255)                                                                           ttgttCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAg                                                                            <  6:191949/147‑1 (MQ=255)                                                                              ttCATGCTGATGGTGACGGTGGCCTTCTTGGCTCGATAGTCGGTGATCAGGTTGTCGTCTCCGGTCCACCAGAGGAGGCGGTCGCTCGGCAGGCCTTGCTGACCGGTGTCCGGGTTGAAGTACATGCAGTGGTCGCCGTTGAAGC‑Tgcg                                                                      <  8:287068/149‑1 (MQ=255)                                                                                 aTGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACC‑CAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGc                                                                  <  6:363661/149‑1 (MQ=255)                                                                                 aTGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACC‑CAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGc                                                                  <  6:368778/149‑1 (MQ=255)                                                                                    cTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCca                                                                >  7:243066/1‑149 (MQ=255)                                                                                    cTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTACGGGTTGAAGTACAGCCAGTCGTCGCAGTTCAAGC‑TGCGAAGCca                                                                >  7:134623/1‑149 (MQ=255)                                                                                     tGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCca                                                                >  3:143556/1‑148 (MQ=255)                                                                                     tGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGGCGTCACCGGTCCCCCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCca                                                                >  1:266467/1‑148 (MQ=255)                                                                                      gAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCcaca                                                              >  6:244552/1‑149 (MQ=255)                                                                                       aaGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAg                                                             >  3:335515/1‑149 (MQ=255)                                                                                                     cttcttGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCg                                               >  7:109776/1‑149 (MQ=255)                                                                                                           gccgccATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGCGTGCGCAGCCACAGCTCCGGGGTTTGCGGGgt                                           <  1:39224/148‑1 (MQ=255)                                                                                                            ccgccATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTc                                        <  8:126932/149‑1 (MQ=255)                                                                                                                       ggTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCg                                                                                  >  5:40388/1‑97 (MQ=255)                                                                                                                       ggTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCg                                                                                  <  6:40388/97‑1 (MQ=255)                                                                                                                        gTGATCAGCTTGTCGTCACAGGTCCAACACAGCAGGCGGTCGATCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAATC‑TGAGCAGCCACAGCTCCGGGGTTTGCGGGGTGTAGGCCAGCAGGCg                            >  8:273548/1‑149 (MQ=255)                                                                                                                             cAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCggg                                             <  6:88603/127‑1 (MQ=255)                                                                                                                             cAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCggg                                             >  5:88603/1‑127 (MQ=255)                                                                                                                             cAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATgg                       >  3:247101/1‑149 (MQ=255)                                                                                                                              aGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCTCAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGgt                      >  8:169441/1‑149 (MQ=255)                                                                                                                              aGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATgg                       <  5:244552/148‑1 (MQ=255)                                                                                                                                gTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCgg                                              >  1:310562/1‑123 (MQ=255)                                                                                                                                gTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCgg                                              <  2:310562/123‑1 (MQ=255)                                                                                                                                   gtcgtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGc                  <  8:243066/148‑1 (MQ=255)                                                                                                                                   gtcgtcACCGGGCCACCACAGCAGGCGGTCGCTCGGCAGGACGTGCTGTCCGGTGTCAGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCATCCACAGCTCCGGGGGTTGCGGGGTGTCGGCCAGCAGGCGGTTGGTGTGc                  <  8:134623/148‑1 (MQ=255)                                                                                                                                    tcgtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGCTTTGCGGGGTGTCGGCCAGCAGGGGGATGGAGTGCAc                >  8:342866/1‑149 (MQ=255)                                                                                                                                     cgtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGgt                      >  2:186113/1‑142 (MQ=255)                                                                                                                                     cgtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGgt                      <  1:186113/142‑1 (MQ=255)                                                                                                                                      gtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCcac                                                               <  3:61894/100‑1 (MQ=255)                                                                                                                                      gtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCcac                                                               >  4:61894/1‑100 (MQ=255)                                                                                                                                      gtcACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCCGCAGGCGGATGGtgtggccct              >  6:259509/1‑144 (MQ=255)                                                                                                                                         aCCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCtttt           <  4:335515/149‑1 (MQ=255)                                                                                                                                         aCCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCttt            <  8:109776/148‑1 (MQ=255)                                                                                                                                         aCCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGCGCACCtttt           >  2:2401/1‑149 (MQ=255)                                                                                                                                               ccaccaCAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCTTTTCCATCg     <  4:247101/149‑1 (MQ=255)                                                                                                                                                 accacAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCTTTTCCATCgg    >  2:50029/1‑148 (MQ=255)                                                                                                                                                 accacAGCAGGCGGTCGCTCCGCAGGCATTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGACAGCAGGCGGATGGTGTGCACCTTTTCCATCgg    >  8:195619/1‑148 (MQ=255)                                                                                                                                                  ccacAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCTTTTCCATCGGCa  >  8:152482/1‑149 (MQ=255)                                                                                                                                                  ccacAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCTTTTCCATCGGCa  >  6:311903/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCGAGTATTCGAGGAAGTCGGCGTCGGTGTTCTCGTCGGTCAGCTTGGCCAGGCGGATGGCGTTCATCTCGCTGTTGTTCATGCTGAAGGTGACGGTGGCCTTCTTGCCGCCATCGACGGTGATCAGGTTGTCGTCACCGGTCCACCACAGCAGGCGGTCGCTCGGCAGGCCTTGCTCACCGGTGTCCGGGTTGAAGTACAGCCAGTCGTCGCCGTTGAAGC‑TGCGCAGCCACAGCTCCGGGGTTTGCGGGGTGTCGGCCAGCAGGCGGATGGTGTGCACCTTTTCCATCGGCA  >  NC_002947/2659429‑2659722

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