Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,339,132	C→G	7.4%	H120Q (CAC→CAG)	pcnB →	poly(A) polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,339,132	0	C	G	7.4%	170.9 / 2.3	54	H120Q (CAC→CAG)	pcnB	poly(A) polymerase
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base G (2/2);  total (27/27)
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.54e-01

GCCGAAGGACTTCGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGTACTACGATCCGGTCAGC  >  NC_002947/5338985‑5339279                                                                                                                                                    |                                                                                                                                                     gCCGAAGGACTTCGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCGTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCccac                                                                                                                                                      >  4:268988/1‑148 (MQ=255)    gAAGGACTTCGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACt                                                                                                                                                  >  3:45312/1‑149 (MQ=255)         aCTTCGACGTCGCCACCAGCGCCACCCCCGAAAAGCTCCGTCCCGAGTTCCGCAACGCCCGCATCAAAGTTCGTCGCTTCAAATTCGTCCATGTTCATATCGGCAGAGAGATCATTGAAGTCGCCACCTTCCGAGCCCACCAATcc                                                                                                                                                >  4:43833/1‑146 (MQ=255)           ttCGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAgg                                                                                                                                            >  3:100895/1‑148 (MQ=255)           ttCGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGga                                                                                                                                           >  3:163799/1‑149 (MQ=255)             cGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAgg                                                                                                                                            <  4:100895/146‑1 (MQ=255)             cGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGgac                                                                                                                                          <  7:57286/148‑1 (MQ=255)               aCGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGgacgac                                                                                                                                       >  7:277255/1‑149 (MQ=255)                  tCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCccac                                                                                                                                                      >  4:339727/1‑131 (MQ=255)                  tCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCccac                                                                                                                                                      <  3:339727/131‑1 (MQ=255)                         caGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGAATCATTGGTCGTCGCTTCAAATTGGTCCATGTTCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCATGGCGATg                                                                                                                             <  1:83135/149‑1 (MQ=255)                                 cccccGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCAGCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGt                                                                                                                    >  1:12217/1‑149 (MQ=255)                                 cccccGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCAGCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGt                                                                                                                    >  1:6356/1‑149 (MQ=255)                                  ccccGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGtt                                                                                                                   >  2:324352/1‑149 (MQ=255)                                   cccGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtc                                                                                                                  <  6:149650/149‑1 (MQ=255)                                     cGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTTAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtcgt                                                                                                                >  4:181591/1‑149 (MQ=255)                                     cGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtcgt                                                                                                                <  4:45312/149‑1 (MQ=255)                                      gAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtcgtc                                                                                                               <  2:315687/149‑1 (MQ=255)                                      gAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtcgtc                                                                                                               >  1:315687/1‑149 (MQ=255)                                                   ccGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCAGCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGt                                                                                                  <  2:6356/149‑1 (MQ=255)                                                   ccGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCAGCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGt                                                                                                  <  2:12217/149‑1 (MQ=255)                                                     gAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtcgtc                                                                                                               <  8:172591/134‑1 (MQ=255)                                                     gAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTtcgtc                                                                                                               >  7:172591/1‑134 (MQ=255)                                                      aGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGc                                                                                               >  7:84873/1‑149 (MQ=255)                                                        ttCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCg                                                                                             >  5:325588/1‑149 (MQ=255)                                                          ccGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGc                                                                                            <  3:162329/148‑1 (MQ=255)                                                                   ccGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGc                                                                                     >  8:219506/1‑146 (MQ=255)                                                                             ggTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTa                                                                        >  8:57302/1‑149 (MQ=255)                                                                              gtcgtcGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTa                                                                        <  6:325588/148‑1 (MQ=255)                                                                              gtcgtcGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTAc                                                                       <  7:219506/149‑1 (MQ=255)                                                                              gtcgtcGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTAc                                                                       <  7:57302/149‑1 (MQ=255)                                                                                    gCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCAcc                                                                 <  8:84873/149‑1 (MQ=255)                                                                                     cTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCAccc                                                                <  1:324352/149‑1 (MQ=255)                                                                                     cTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTACGCGCCCACCACTACGAGGACGACCAGGCCGATGCCAA‑CCGTTCGTCGCATAACGACAGTGGCCGCATTCTGCGAGACAACGTGTACGGCAccc                                                                >  3:100526/1‑149 (MQ=255)                                                                                      ttCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGc                                                                                     <  2:11114/127‑1 (MQ=255)                                                                                      ttCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGc                                                                                     >  1:11114/1‑127 (MQ=255)                                                                                      ttCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCgaga                                                                                 <  5:221/131‑1 (MQ=255)                                                                                      ttCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCgaga                                                                                 >  6:221/1‑131 (MQ=255)                                                                                       tCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGa                                                                                                                                              <  5:49738/70‑1 (MQ=255)                                                                                       tCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGa                                                                                                                                              >  6:49738/1‑70 (MQ=255)                                                                                         aaaTTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCTCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTTGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTgg                                                             <  8:275917/148‑1 (MQ=255)                                                                                         aaaTTGGTCCATGTGCATTTCGGCAGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCAGCATTCTGCGAGACAACGTGTACTGCACCCTgg                                                             >  1:201652/1‑148 (MQ=255)                                                                                            ttGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAaga                                                         <  8:277255/149‑1 (MQ=255)                                                                                            ttGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAaga                                                         <  3:268988/149‑1 (MQ=255)                                                                                              ggTCCATGTGCATTTCTGCCTTTATATCATTGAAGTCGCCAGCTTCCGCGCCCACCTCTCCGAGGACGACCAGGGCGATGCTCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCGTGGAagag                                                        <  3:43833/148‑1 (MQ=255)                                                                                              ggTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAaga                                                         >  7:92264/1‑147 (MQ=255)                                                                                              ggTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAaga                                                         <  8:92264/147‑1 (MQ=255)                                                                                                    tgtgCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCACCCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCga                                                                                   <  1:244340/116‑1 (MQ=255)                                                                                                    tgtgCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCga                                                                                   >  2:244340/1‑115 (MQ=255)                                                                                                             cGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCa                                                                                                    <  6:57365/89‑1 (MQ=255)                                                                                                             cGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCa                                                                                                    >  5:57365/1‑89 (MQ=255)                                                                                                              ggCCGTGAGATCATTTAAGTAGCCACCTTCCGCGCCCACCAATCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGTCACCCTTGAAGAGGACGCGCAACGCCGCGa                                       <  4:100526/149‑1 (MQ=255)                                                                                                                            tGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACg                         >  6:189868/1‑149 (MQ=255)                                                                                                                                   gCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCACATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGt                   >  8:134721/1‑148 (MQ=255)                                                                                                                                    ccACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGtac                 >  5:126453/1‑149 (MQ=255)                                                                                                                                     cACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGtac                 >  3:146198/1‑148 (MQ=255)                                                                                                                                             gcgcCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGTACTACGATCCg        <  5:189868/149‑1 (MQ=255)                                                                                                                                                 ccaccaCTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTTGCCGCATTCAGCGAGACAACGTGGACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGTACTACGATCCGGTCa    <  8:243845/149‑1 (MQ=255)                                                                                                                                                 ccaccaCTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCTACGCCCTGTACTACGATCCGGTCa    <  2:201652/149‑1 (MQ=255)                                                                                                                                                  caccacTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGTACTACGATCCGGTCag   >  2:208869/1‑149 (MQ=255)                                                                                                                                                   accacTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCAAAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGTACTACGATCCGGTCagc  >  2:59214/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCCGAAGGACTTCGACGTCGCCACCAGCGCCACCCCCGAACAGGTCCGTGCCGAGTTCCGCAACGCCCGCATCATTGGTCGTCGCTTCAAATTGGTCCATGTGCATTTCGGCCGTGAGATCATTGAAGTCGCCACCTTCCGCGCCCACCACTCCGAGGACGACCAGGGCGATGCCCA‑CCGTTCGTCGCATAACGCCAGTGGCCGCATTCTGCGAGACAACGTGTACGGCACCCTGGAAGAGGACGCGCAACGCCGCGACTTCACCATCAACGCCCTGTACTACGATCCGGTCAGC  >  NC_002947/5338985‑5339279

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