Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	833,979	C→G	6.1%	H219D (CAC→GAC)	PP_0717 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	833,979	0	C	G	6.1%	205.3 / 3.7	66	H219D (CAC→GAC)	PP_0717	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base G (2/2);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

CGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTG  >  NC_002947/833834‑834122                                                                                                                                                  |                                                                                                                                                cGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACg                                                                                                                                              >  1:280789/1‑149 (MQ=255)     cATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAg                                                                                                                                          <  4:50038/149‑1 (MQ=255)               gggTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGa                                                                                                                                    <  1:312901/145‑1 (MQ=255)               gggTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGa                                                                                                                                    >  2:312901/1‑145 (MQ=255)                 gTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGAt                                                                                                                              >  6:299084/1‑149 (MQ=255)                    gcgcTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCgg                                                                                                                           <  1:365135/149‑1 (MQ=255)                     cgcTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGt                                                                                                                          <  1:38993/149‑1 (MQ=255)                     cgcTGCGTGCCACGTTGATCGCGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGt                                                                                                                          <  3:147250/149‑1 (MQ=255)                      gcTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGtt                                                                                                                         >  8:217971/1‑149 (MQ=255)                        tGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTc                                                                                                                        >  4:149707/1‑148 (MQ=255)                              ccACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACcttg                                                                                                                 >  1:150869/1‑149 (MQ=255)                              ccACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACcttg                                                                                                                 <  2:280789/149‑1 (MQ=255)                               cACGGTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACcttgc                                                                                                                <  8:141224/149‑1 (MQ=255)                                aCGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACcttgct                                                                                                               >  4:43154/1‑149 (MQ=255)                                           gATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGgtgt                                                                                                    >  6:95309/1‑149 (MQ=255)                                             tGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGGTCAACCTTGCTTGCCGGGGTGTTg                                                                                                  <  5:95309/149‑1 (MQ=255)                                              gCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCTCGCGCGGCACGCCAGCCGGGAGCt                                                                                                                                 <  7:234720/117‑1 (MQ=255)                                              gCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCt                                                                                                                                 >  8:234720/1‑117 (MQ=255)                                                cGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGGACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGc                                                                                                >  2:143616/1‑148 (MQ=255)                                                cGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGGACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGc                                                                                                >  2:138103/1‑148 (MQ=255)                                                 ggCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTgg                                                                                                 >  3:164351/1‑146 (MQ=255)                                                 ggCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTgg                                                                                                 <  4:164351/146‑1 (MQ=255)                                                  gCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGAGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCATGCTTGCAGCGGGCTTGGCGa                                                                                              >  4:44968/1‑148 (MQ=255)                                                    cTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGAGCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAACCCCGCGCGCGGCACGCCAGCCGGGAGCTGAACGGTACAACATTGCTTCCCGGGTTGTTGGCGATAg                                                                                           >  5:185546/1‑149 (MQ=255)                                                     tttCTGCAGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAgc                                                                                          >  3:210105/1‑149 (MQ=255)                                                           ctgGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGTGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCCGGACCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCtgtg                                                                                    >  2:179604/1‑149 (MQ=255)                                                           ctgGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCggg                                                                                                        >  1:111644/1‑129 (MQ=255)                                                           ctgGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCggg                                                                                                        <  2:111644/129‑1 (MQ=255)                                                           ctgGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCtg                                                                                      <  5:333859/147‑1 (MQ=255)                                                              gCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCCACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGt                                                                                 >  8:260183/1‑149 (MQ=255)                                                              gCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTgg                                                                                   <  2:289883/147‑1 (MQ=255)                                                                   gATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAgcg                                                                             <  5:299084/148‑1 (MQ=255)                                                                         cccTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCtgt                                                                                     <  3:80508/134‑1 (MQ=255)                                                                         cccTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCtgt                                                                                     >  4:80508/1‑134 (MQ=255)                                                                              ctCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTgt                                                                   <  7:197848/146‑1 (MQ=255)                                                                              aTCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTgt                                                                   <  3:149707/147‑1 (MQ=255)                                                                                cAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTgttcgt                                                               >  5:49979/1‑149 (MQ=255)                                                                                         ccGTTGATCATGAAGGCCGTGAGTGTGGGGCAACTAGCATGCGAAACCCGCGCGCGGCACGCCAGCCGGTTGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGTCGATAGGGGTGTGTGTAGCGCTGAGCCTGTTCGTTCACCtgtg                                                      <  1:179604/149‑1 (MQ=255)                                                                                              gATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCggg                                                                                                        <  5:198068/94‑1 (MQ=255)                                                                                              gATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCggg                                                                                                        >  6:198068/1‑94 (MQ=255)                                                                                              gATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGt                                                 >  4:26714/1‑149 (MQ=255)                                                                                                tCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCt                                                                           <  5:83657/121‑1 (MQ=255)                                                                                                tCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCt                                                                           >  6:83657/1‑121 (MQ=255)                                                                                                 cATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGa                                                                         >  8:353968/1‑122 (MQ=255)                                                                                                 cATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGa                                                                         <  7:353968/122‑1 (MQ=255)                                                                                                   tGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGGACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGttcttc                                            <  1:138103/149‑1 (MQ=255)                                                                                                   tGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGGACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGttcttc                                            <  1:143616/149‑1 (MQ=255)                                                                                                         ccGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTgg                                      <  4:210105/149‑1 (MQ=255)                                                                                                          cGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTggg                                     >  6:200533/1‑149 (MQ=255)                                                                                                          cGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTggg                                     <  3:43154/149‑1 (MQ=255)                                                                                                               gtgtGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTcttctt                                          <  4:199970/139‑1 (MQ=255)                                                                                                               gtgtGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTcttctt                                          >  3:199970/1‑139 (MQ=255)                                                                                                                tgtgGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTActgc                               >  5:337448/1‑149 (MQ=255)                                                                                                                tgtgGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTActgc                               >  1:256463/1‑149 (MQ=255)                                                                                                                       caacaaGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTActg                                <  3:138584/141‑1 (MQ=255)                                                                                                                       caacaaGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTActg                                >  4:138584/1‑141 (MQ=255)                                                                                                                       caacaaGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGc                        >  8:277139/1‑149 (MQ=255)                                                                                                                        aacaaGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGcc                       <  2:256463/149‑1 (MQ=255)                                                                                                                        aacaaGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGGGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGcc                       <  2:36355/149‑1 (MQ=255)                                                                                                                          caaGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCt                     >  1:72026/1‑149 (MQ=255)                                                                                                                            aGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCttt                   <  6:185546/149‑1 (MQ=255)                                                                                                                                 gCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTtgt                                        >  7:334544/1‑123 (MQ=255)                                                                                                                                 gCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTtgt                                        <  8:334544/123‑1 (MQ=255)                                                                                                                                  cGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAgg             <  7:217971/149‑1 (MQ=255)                                                                                                                                      aCCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCg         >  8:347083/1‑149 (MQ=255)                                                                                                                                              gcgGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTg  >  2:265479/1‑148 (MQ=255)                                                                                                                                                  cACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTtgt                                        <  6:383829/106‑1 (MQ=255)                                                                                                                                                  cACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTtgt                                        >  5:383829/1‑106 (MQ=255)                                                                                                                                                  |                                                                                                                                                CGAGCATGTCAGCTGGGTGGCGCTGCGTGCCACGTTGATCGTGATGCCGGCCTTTCTGCTGGCATTGATCGCCCCTGATCAGTACATGCCGTTGATCATGAAGGCCGTGAGTGTGGGGCAACAAGCAGGCGAAACCCGCGCGCGGCACGCCAGCCGGGAGCTGATCGGTTCAACCTTGCTTGCCGGGGTGTTGGCGATAGCGCTGTGGGTAGCGCTGAGCCTGTTCGTTCACCTGTGGATGTTCTTCTTGTGGGTACTGCTGTTCGCCCTTTGGCAGGCTCGGCGCCTG  >  NC_002947/833834‑834122

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