Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,053,362	T→A	5.6%	K80I (AAA→ATA)	PP_1828 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,053,362	0	T	A	5.6%	214.9 / 2.7	71	K80I (AAA→ATA)	PP_1828	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/34);  new base A (2/2);  total (35/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 = 9.58e-01

CACTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAGGCGTAGCAG  >  NC_002947/2053215‑2053507                                                                                                                                                       |                                                                                                                                                  cACTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCTTGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCGGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCttt                                                                                                                                                  <  3:301755/149‑1 (MQ=255) cACTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCttt                                                                                                                                                  <  7:60688/149‑1 (MQ=255)  aCTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGACGTCCAGGCGGTTCTGCTCTTtg                                                                                                                                                 <  8:33669/149‑1 (MQ=255)   cTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTtgt                                                                                                                                                >  1:271801/1‑149 (MQ=255)   cTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTtgt                                                                                                                                                >  2:87626/1‑149 (MQ=255)   cTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCACTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTtgt                                                                                                                                                <  6:339575/149‑1 (MQ=255)    ttGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTtgtg                                                                                                                                               >  6:4610/1‑149 (MQ=255)      gcgcTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCACGGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTtgt                                                                                                                                                <  1:278304/147‑1 (MQ=255)          tGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTATGTGCACttt                                                                                                                                         >  4:104464/1‑149 (MQ=255)          tGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTATGTGCACtt                                                                                                                                          >  4:99172/1‑148 (MQ=255)            cAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACttt                                                                                                                                         <  6:338991/147‑1 (MQ=255)               ccccTGTGCGGCCTGGGCCAGCATGAAAGGGCCGATGCCACGGTTCTCCGGTTGCAGTTGCAGGCCATGTTCATTTCCAATA‑CGCTCCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGcc                                                                                                                                    <  6:201246/146‑1 (MQ=255)                 ccTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCa                                                                                                                                  <  7:262445/148‑1 (MQ=255)                     tgCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCTTGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGc                                                                                                                               <  6:146949/148‑1 (MQ=255)                     tgCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCAAATA‑CGCACCA‑GGCACTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCGTTGTGCACTTTCTGCCCAGGc                                                                                                                               >  2:115716/1‑148 (MQ=255)                        ggCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACt                                                                                                                           <  4:202528/149‑1 (MQ=255)                            tGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCGTTGTGCACTTTCTGCCCAGGCAACTGGAt                                                                                                                       >  2:103438/1‑149 (MQ=255)                              ggCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGAt                                                                                                                       <  2:271801/147‑1 (MQ=255)                                ccAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTcc                                                                                                                   <  5:4610/149‑1 (MQ=255)                                  aggaTGAAACGGCCGATGCCACGGTTCAACGGTTCCAGTTGCAGGCGATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATTGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCGAGGCAACTGGATTTCCAt                                                                                                                 <  4:179194/146‑1 (MQ=255)                                  aGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCAt                                                                                                                 >  6:232430/1‑149 (MQ=255)                                         aaCGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCgcagc                                                                                                          <  5:119041/149‑1 (MQ=255)                                         aaCGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCgcag                                                                                                           >  4:327767/1‑148 (MQ=255)                                             gCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGc                                                                                                      >  3:71389/1‑149 (MQ=255)                                             gCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGc                                                                                                      >  8:275450/1‑149 (MQ=255)                                                  tGCCACGGTTCACCGGTTCCAGTTGCAGGCCTTGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCtgtgt                                                                                                 <  1:81612/149‑1 (MQ=255)                                                  tGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCTCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCtgtgt                                                                                                 <  3:279095/149‑1 (MQ=255)                                                     cACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGcc                                                                                              >  2:42806/1‑149 (MQ=255)                                                      aCGGTGCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCt                                                                                             >  5:122564/1‑149 (MQ=255)                                                       cGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTG‑CAATATCGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCtt                                                                                            <  1:315624/149‑1 (MQ=255)                                                        ggTGCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTg                                                                                           <  5:4837/149‑1 (MQ=255)                                                          ttCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTc                                                                                         <  4:163398/149‑1 (MQ=255)                                                            cACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCAt                                                                                       >  8:81224/1‑149 (MQ=255)                                                             aCCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATg                                                                                      <  7:144301/149‑1 (MQ=255)                                                             aCCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATg                                                                                      >  2:377181/1‑149 (MQ=255)                                                              ccGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATgg                                                                                     <  6:122564/149‑1 (MQ=255)                                                              ccGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATgg                                                                                     <  5:307671/149‑1 (MQ=255)                                                                 gTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGgcgc                                                                                  >  2:299687/1‑149 (MQ=255)                                                                 gTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGgcgc                                                                                  >  5:349658/1‑149 (MQ=255)                                                                  ttCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCt                                                                                 >  6:289548/1‑149 (MQ=255)                                                                        ttGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCAt                                                                           <  3:327767/149‑1 (MQ=255)                                                                          gCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAg                                                                         >  4:185551/1‑149 (MQ=255)                                                                              gCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGCCGCGGATGTCCAGGCGGTTCTGCTCTGTGTGCACTTTCTGCGCAGGCAACTGGATTTCGATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGcggc                                                                     >  4:248519/1‑146 (MQ=255)                                                                              gCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCg                                                                      >  1:341867/1‑148 (MQ=255)                                                                               ccATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGGACATGGTCGACGCGGATGTCCAGGCGGTTCTGCTATTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCg                                                                      >  1:242964/1‑148 (MQ=255)                                                                                 atcttCATTGCGAATA‑CCCACCA‑GGCTCTCGG‑CATGGTAGAAGCGGATGTACCGGCTGTGCTGCTCTTTCTGCACTCTCTGCCCCGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTCCCTTGTCCTGGCGCTCTGCATAGCCGAACt                                                                  <  4:284272/146‑1 (MQ=255)                                                                                 aTGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACt                                                                  <  1:87626/149‑1 (MQ=255)                                                                                   gTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTg                                                                 >  7:189526/1‑148 (MQ=255)                                                                                   gTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTg                                                                 >  7:130240/1‑148 (MQ=255)                                                                                       aTTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATCTCGAGCCGGTTCTGCTCTTTGTGCACTATCTGCCCAGGCAACTCCATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGcaca                                                             >  6:200389/1‑148 (MQ=255)                                                                                        ttGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACaaa                                                           <  1:377181/149‑1 (MQ=255)                                                                                        ttGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACaa                                                            >  7:240023/1‑148 (MQ=255)                                                                                          gCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACaaaa                                                          <  8:240023/148‑1 (MQ=255)                                                                                           ccAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAAc                                                         <  6:349658/148‑1 (MQ=255)                                                                                            cAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCTCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCAt                                                                           <  2:286135/129‑1 (MQ=255)                                                                                            cAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCAt                                                                           >  1:286135/1‑129 (MQ=255)                                                                                            cAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACg                                                        >  6:371477/1‑148 (MQ=255)                                                                                             aaTA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCa                                                      <  3:365084/149‑1 (MQ=255)                                                                                               tA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGc                                                    <  5:56395/149‑1 (MQ=255)                                                                                                     aCCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCcg                                               >  7:287854/1‑149 (MQ=255)                                                                                                      ccA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTACACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCcgc                                              >  5:114105/1‑149 (MQ=255)                                                                                                            cgcTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTcc                                         <  4:71389/149‑1 (MQ=255)                                                                                                                         tCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGcggtc                             >  8:309583/1‑149 (MQ=255)                                                                                                                               cgGATGTCCAGGCGGTTCTGCTCTATGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGtcggtcgg                       <  3:104464/149‑1 (MQ=255)                                                                                                                               cgGATGTCCAGGCGGTTCTGCTCTATGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGtcggtcgg                       <  3:99172/149‑1 (MQ=255)                                                                                                                                     tCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGtg                  <  8:130240/148‑1 (MQ=255)                                                                                                                                           cGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAg           <  1:299687/149‑1 (MQ=255)                                                                                                                                           cGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAg           <  2:141923/149‑1 (MQ=255)                                                                                                                                                   ctctTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAGGCGTag     >  6:150175/1‑147 (MQ=255)                                                                                                                                                   ctctTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAGGCGTag     >  8:271831/1‑147 (MQ=255)                                                                                                                                                    tctTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCggtcggtc                         >  6:146541/1‑126 (MQ=255)                                                                                                                                                    tctTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCggtcggtc                         <  5:146541/126‑1 (MQ=255)                                                                                                                                                    tctTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAGGCGTagcag  >  1:52687/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CACTTGCGCTGCAACCACTGTGCGGCCTGGGCCAGCATGAAACGGCCGATGCCACGGTTCACCGGTTCCAGTTGCAGGCCATGTTCATTGCCAATA‑CGCACCA‑GGCGCTCGG‑CATGGTCGACGCGGATGTCCAGGCGGTTCTGCTCTTTGTGCACTTTCTGCCCAGGCAACTGGATTTCCATGCGCAGCAGGCTGTGTGCCTTGTCATGGCGCTCGGCATAGCCGAACTGCACAAAACGCAGCGGCCGCGCTCCGGTGCTGCGGTCGGTCGGCAGTGCTGCCAGGCGTAGCAG  >  NC_002947/2053215‑2053507

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