Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	654,406	T→G	6.8%	K114Q (AAG→CAG)	ctpL ←	methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	654,406	0	T	G	6.8%	166.9 / 4.2	59	K114Q (AAG→CAG)	ctpL	methyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base G (2/2);  total (27/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 = 9.78e-01

GATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGT  >  NC_002947/654259‑654552                                                                                                                                                     |                                                                                                                                                   gATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGc                                                                                                                                                  <  2:1207489/150‑1 (MQ=255)    gCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCcgg                                                                                                                                               <  2:328953/150‑1 (MQ=255)       gggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCAAGCTGGCcggcga                                                                                                                                            <  2:317739/150‑2 (MQ=255)        ggTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCggcggc                                                                                                                                           <  1:603249/150‑1 (MQ=255)           ggCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCcag                                                                                                                                        <  2:1208734/150‑1 (MQ=255)            gCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAATTCGGCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCcagc                                                                                                                                       >  1:221785/1‑150 (MQ=255)                 tAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTc                                                                                                                                  <  2:120753/150‑1 (MQ=255)                        ggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCg                                                                                                                           <  2:775980/150‑1 (MQ=255)                        ggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCg                                                                                                                           <  1:933328/150‑1 (MQ=255)                            gcaAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGa                                                                                                                       >  2:992211/1‑150 (MQ=255)                                gTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGc                                                                                                                   <  1:43491/150‑1 (MQ=255)                                     gATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGACCCCAGCTCCCGCTCGGCTTCCAGCTTGACGTC‑TGTGAGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCa                                                                                                              >  1:611213/1‑150 (MQ=255)                                      aTTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAg                                                                                                             <  2:285316/150‑1 (MQ=255)                                            cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGCCGCCGGCAAGCTGGCcggcgacc                                                                                                                                          >  1:696722/1‑112 (MQ=255)                                            cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGCCGCCGGCAAGCTGGCCGGCGAcc                                                                                                                                          <  2:696722/115‑1 (MQ=255)                                                   ttGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTg                                                                                                <  1:574943/150‑1 (MQ=255)                                                    tGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                               <  1:1119045/150‑1 (MQ=255)                                                     gCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGc                                                                                                               <  1:460330/133‑1 (MQ=255)                                                     gCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGc                                                                                                               >  2:460330/1‑133 (MQ=255)                                                          tGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGcc                                                                                                                                                 <  1:1135812/94‑1 (MQ=255)                                                          tGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGcc                                                                                                                                                 >  2:1135812/1‑94 (MQ=255)                                                            tCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAAc                                                                                       >  2:825213/1‑150 (MQ=255)                                                                       gCTGGCCATCTGCTCGAAGTTCGCCCCCACCTACCGCTAGGCTAGCCACAGCAGGGCGTTTG‑GGTCGCCGGCCCGTATTCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAg                                                                            <  1:606880/150‑1 (MQ=255)                                                                        cTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGc                                                                           >  2:1250112/1‑150 (MQ=255)                                                                           gCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGACCAGctct                                                                        >  1:416693/1‑150 (MQ=255)                                                                                cTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCgg                                                                   >  1:398428/1‑150 (MQ=255)                                                                                 tGCTCGAAGTTCGACCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                                  <  2:611213/150‑1 (MQ=255)                                                                                    tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGc                                                                                                                                                  >  1:494466/1‑67 (MQ=255)                                                                                    tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGc                                                                                                                                                  <  2:494466/67‑1 (MQ=255)                                                                                      gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTc                                                             <  1:663780/150‑1 (MQ=255)                                                                                      gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTc                                                             <  1:133425/150‑1 (MQ=255)                                                                                           tCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAACTCTACCGGCCAGTCCGTCa                                                        >  2:70999/1‑150 (MQ=255)                                                                                                cccAGCTCCCGCTCTGCTGGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                       <  1:1115481/98‑1 (MQ=255)                                                                                                cccAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                       >  2:1115481/1‑98 (MQ=255)                                                                                                  cAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                                 <  2:1134323/150‑1 (MQ=255)                                                                                                   aGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGcc                                                <  2:350542/150‑1 (MQ=255)                                                                                                     cTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTg                                              >  1:833385/1‑150 (MQ=255)                                                                                                          gCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                                 >  1:276246/1‑142 (MQ=255)                                                                                                          gCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGc                                                 <  2:276246/142‑1 (MQ=255)                                                                                                          gCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACt                                         >  1:568500/1‑150 (MQ=255)                                                                                                                 ttGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGtt                                  >  2:624567/1‑150 (MQ=255)                                                                                                                    cagcagcagGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGc                               >  2:881250/1‑150 (MQ=255)                                                                                                                       cagcagGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtc                                                                                    <  2:862714/94‑1 (MQ=255)                                                                                                                       cagcagGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtc                                                                                    >  1:862714/1‑94 (MQ=255)                                                                                                                          cagGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                             >  2:1258435/1‑82 (MQ=255)                                                                                                                          cagGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgggg                                                                                             <  1:1258435/82‑1 (MQ=255)                                                                                                                           ggggCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATg                        >  2:882665/2‑150 (MQ=255)                                                                                                                             ggCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCg                                                                                                                           >  1:303011/1‑49 (MQ=255)                                                                                                                             ggCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCg                                                                                                                           <  2:303011/49‑1 (MQ=255)                                                                                                                                ttGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTTAAGGCCTGAAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggc                   >  1:670376/1‑150 (MQ=255)                                                                                                                                    g‑ggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTgg                                    <  1:267922/129‑1 (MQ=255)                                                                                                                                    g‑ggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTgg                                    >  2:267922/1‑129 (MQ=255)                                                                                                                                          gCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGgcg          >  2:895784/1‑150 (MQ=255)                                                                                                                                            cGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCa        <  2:221785/150‑1 (MQ=255)                                                                                                                                              gCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg      >  1:1144587/1‑150 (MQ=255)                                                                                                                                              gCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg      <  1:992211/150‑1 (MQ=255)                                                                                                                                              gCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg      <  1:895784/150‑1 (MQ=255)                                                                                                                                                cAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGc    <  2:1144587/150‑1 (MQ=255)                                                                                                                                                  gCTGGCCGCCGGCAAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGt  >  1:516410/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGCTTGTG‑GGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGGGCGT  >  NC_002947/654259‑654552

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑