Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,068,527:1	+G	11.1%	coding (693/948 nt)	PP_0926 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,068,527	1	.	G	11.1%	185.2 / 20.4	54	coding (693/948 nt)	PP_0926	membrane protein
Reads supporting (aligned to +/- strand):  ref base . (25/23);  new base G (3/3);  total (28/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 = 7.44e-01

CCATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCGGGCAG  >  NC_002947/1068379‑1068672                                                                                                                                                         |                                                                                                                                                  ccATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAt                                                                                                                                                    >  3:348249/1‑149 (MQ=255)       cTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGG‑CATCGATGGCCAGc                                                                                                                                             <  1:104225/149‑1 (MQ=255)             gCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTg                                                                                                                                       <  5:190190/149‑1 (MQ=255)                           cccAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCt                                                                                                                         <  8:230171/149‑1 (MQ=255)                                cTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCgc                                                                                                                    <  7:302179/149‑1 (MQ=255)                                     tACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTc                                                                                                               >  4:228797/1‑149 (MQ=255)                                          aCCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGc                                                                                                            <  8:164687/147‑1 (MQ=255)                                            cTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGG‑CATCGATGGCCAGc                                                                                                                                             <  5:215065/112‑1 (MQ=255)                                            cTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGG‑CATCGATGGCCAGc                                                                                                                                             >  6:215065/1‑112 (MQ=255)                                                ggCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTgg                                                                                                    <  4:348249/149‑1 (MQ=255)                                                     cTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGTGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCt                                                                                               >  2:254285/1‑149 (MQ=255)                                                      tACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTg                                                                                              >  1:219937/1‑149 (MQ=255)                                                      tACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTg                                                                                              >  1:258298/1‑149 (MQ=255)                                                         cTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                           >  6:215896/1‑149 (MQ=255)                                                               caGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCgg                                                                                     <  7:105773/149‑1 (MQ=255)                                                                 gATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGt                                                                                   >  1:137673/1‑149 (MQ=255)                                                                      cccTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgc                                                                              >  4:38047/1‑149 (MQ=255)                                                                       ccTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgcg                                                                             <  2:137673/149‑1 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgcg                                                                             <  3:77485/147‑1 (MQ=255)                                                                         tGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTcgcg                                                                             >  4:77485/1‑147 (MQ=255)                                                                                ggCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGc                                                                    >  8:41396/1‑149 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGTGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                           >  2:366599/1‑123 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGTGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGtt                                                                                           <  1:366599/123‑1 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                                 <  2:258298/149‑1 (MQ=255)                                                                                   gcTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                                 >  1:293092/1‑149 (MQ=255)                                                                                     ggATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGTT‑GCCTGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGc                                                               <  3:274986/148‑1 (MQ=255)                                                                                       aTCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGGCTACCCGCTGGCCATCGAT‑GCCACCTGCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTgg                                                                 <  1:120895/148‑1 (MQ=255)                                                                                         cTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATcc                                                           >  4:202105/1‑149 (MQ=255)                                                                                         cTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATcc                                                           >  5:169070/1‑149 (MQ=255)                                                                                           gTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGCGCTGGGCATCCtt                                                         >  3:176238/1‑149 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCttt                                                        >  2:42944/1‑149 (MQ=255)                                                                                            tATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCGGGGc                                                               >  4:274986/1‑142 (MQ=255)                                                                                             aTCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCttt                                                        <  1:42944/148‑1 (MQ=255)                                                                                              tCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                      <  1:320325/149‑1 (MQ=255)                                                                                              tCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                      >  1:352560/1‑149 (MQ=255)                                                                                              tCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTc                                                      >  5:347630/1‑149 (MQ=255)                                                                                               cggcggCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTcc                                                     >  7:31800/1‑149 (MQ=255)                                                                                                        gCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAg                                            <  3:38047/149‑1 (MQ=255)                                                                                                        gCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAg                                            <  3:202105/149‑1 (MQ=255)                                                                                                              gCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGcc                                      <  3:350420/149‑1 (MQ=255)                                                                                                                      cgcTGTGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCt                              <  1:254285/149‑1 (MQ=255)                                                                                                                       gcTGCGTAACTGTGG‑T‑‑CGGCTGGCCACCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTg                             <  1:178459/149‑1 (MQ=255)                                                                                                                          gCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCgg                          >  3:270849/1‑149 (MQ=255)                                                                                                                            gTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGcc                        <  6:347630/149‑1 (MQ=255)                                                                                                                             tAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGc                         >  7:172428/1‑147 (MQ=255)                                                                                                                             tAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCa                       <  8:31800/149‑1 (MQ=255)                                                                                                                                 tgtgG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTa                   >  4:280150/1‑149 (MQ=255)                                                                                                                                   tgG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACg                 >  4:342981/1‑149 (MQ=255)                                                                                                                                     g‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGt                >  3:301241/1‑148 (MQ=255)                                                                                                                                        cgaagCTGG‑CATCGATGGCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGt                >  2:128684/5‑148 (MQ=255)                                                                                                                                            gCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATgg           >  2:206228/1‑149 (MQ=255)                                                                                                                                                gCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGgcgc        <  2:293092/148‑1 (MQ=255)                                                                                                                                                   aTCGATGGCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAg                                            >  4:256296/1‑110 (MQ=255)                                                                                                                                                   aTCGATGGCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAg                                            <  3:256296/110‑1 (MQ=255)                                                                                                                                                      gAT‑GCCAGCGTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCGGGCAg  <  5:215896/148‑1 (MQ=255)                                                                                                                                                        t‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCg                                                                        >  3:117054/1‑76 (MQ=255)                                                                                                                                                        t‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCg                                                                        <  4:117054/76‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                  CCATGGCTTACGGCCAGGTGCCGCTGCCCAACTGGCTACTGACCTCGGGCGACTACCTGGCACAGATGACCCTGCCGCTGGCGCTGATCTGTATCGGCGGCACGCTGTCGCTGGCAGCGCTGCGTAACTGTGG‑T‑‑CGGCTGGCCATCGAT‑GCCAGCTTCGTGAAAATGCTCTGGCTGCCGCTGCTCGGCACCCTTGGGGCCTGGTTGTGCGGGTTTCGCGGGGCGGAGCTGGGCATCCTTTTCCTGTACATAGGCAGCCCGACCGCTGCGGCCAGCTACGTGATGGCGCGGGCAG  >  NC_002947/1068379‑1068672

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