Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,187,585	C→G	12.3%	R30R (CGC→CGG) ‡	PP_4567 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,187,585	0	C	G	12.3%	150.8 / 12.3	57	R30R (CGC→CGG) ‡	PP_4567	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base G (4/3);  total (31/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 = 6.29e-01

GGGGCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAG  >  NC_002947/5187438‑5187733                                                                                                                                                    |                                                                                                                                                     ggggCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGcc                                                                                                                                                     <  3:360201/149‑1 (MQ=255)  ggtcATGGCTGCTGAATGAGCGGTTGTCGATGGCGCATCTATACGGCGGACTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATCAACCGCGCGTACAAAAAGGTCCTGTTCCGCCt                                                                                                                                                    <  6:323606/146‑1 (MQ=255)        ggCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGaa                                                                                                                                                >  4:359294/1‑147 (MQ=255)          ctgctgAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAgg                                                                                                                                            <  1:71651/149‑1 (MQ=255)           tgctgAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGa                                                                                                                                           <  4:362960/149‑1 (MQ=255)              tgAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACg                                                                                                                                         >  3:239377/1‑148 (MQ=255)               gAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGAc                                                                                                                                          >  7:68748/1‑146 (MQ=255)                  tGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGACCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGAccttttt                                                                                                                                    <  7:225813/149‑7 (MQ=255)                    aGCCGTTGTCGATGGCTCATCTATACGGCGGCCTTCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATTAACACCGCGTACAAAAAGGTACTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                   <  7:43934/148‑1 (MQ=255)                     gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                   <  4:36018/147‑1 (MQ=255)                     gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACg                                                                                                                                   >  3:36018/1‑147 (MQ=255)                     gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGa                                                                                                                                           <  2:341867/139‑1 (MQ=255)                     gCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCAGGCGGAACAGGa                                                                                                                                           >  1:341867/1‑139 (MQ=255)                       cGTTGTCGATGGCCCATCTATACGGCGGCCTCCTGATTGGACTGGCCGTGTGGCTGGTGCTGCGCCCGGCGCGAAGTAGAAGCGGGGAACTGCATGAACAGCGCGTACAAAACGGTCCTGTTCCGCCTGGGACAGGCCTAGAACGGCTa                                                                                                                               >  4:333760/1‑149 (MQ=255)                         ttGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTAcc                                                                                                                             >  2:33905/1‑149 (MQ=255)                               aTGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCgg                                                                                                                       <  8:238327/149‑1 (MQ=255)                                  gcgcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCt                                                                                                                    >  2:335882/1‑149 (MQ=255)                                    gcATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCg                                                                                                                  >  8:20439/1‑149 (MQ=255)                                          ataCGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAg                                                                                                            >  5:9885/1‑149 (MQ=255)                                          ataCGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAg                                                                                                            >  3:79744/1‑149 (MQ=255)                                           taCGGCGGCCTGCTGATTGGCCTGGCCGTGTGGCGGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAg                                                                                                            <  4:119392/148‑1 (MQ=255)                                              ggcggcCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGcc                                                                                                        >  6:316302/1‑149 (MQ=255)                                                  gcCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCtgtg                                                                                                    <  6:363015/149‑1 (MQ=255)                                                        tgATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGcc                                                                                                >  5:259291/1‑147 (MQ=255)                                                          aTTGGCCTGGCCTTGGGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAg                                                                                             <  3:136391/148‑1 (MQ=255)                                                               cctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTg                                                                                       <  6:259291/149‑1 (MQ=255)                                                                ctggcctTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTg                                                                                       >  4:45770/1‑148 (MQ=255)                                                                     ctTGTGGCTGGTGCTGCGGCCGGGGCGGAGTTGATGCGGGGTACTGCATGAAGAGGGCGTAAAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAAGGGGTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAggc                                                                                 <  4:239377/149‑1 (MQ=255)                                                                        gtgGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAggcgg                                                                               <  7:175851/148‑1 (MQ=255)                                                                           gCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTAt                                                                           <  8:68748/149‑1 (MQ=255)                                                                             ttcaactgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGc                                                                         >  6:80257/6‑149 (MQ=255)                                                                             ttcaactgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGc                                                                         >  5:351803/6‑149 (MQ=255)                                                                                 actgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATc                                                                     >  4:19002/2‑149 (MQ=255)                                                                                 actgcGCGCCGGGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATc                                                                     >  1:365206/2‑149 (MQ=255)                                                                                     cgcCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGac                                                                  <  3:359294/148‑1 (MQ=255)                                                                                                                cTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATc                                       <  1:33905/148‑1 (MQ=255)                                                                                                                   cATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGccc                                   >  5:112478/1‑149 (MQ=255)                                                                                                                      gAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                >  7:335666/1‑149 (MQ=255)                                                                                                                       aaCAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCggg                               <  6:351803/149‑1 (MQ=255)                                                                                                                       aaCAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                <  3:19002/148‑1 (MQ=255)                                                                                                                         cAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGCCATCGCCCcgcgcg                             >  6:117093/1‑149 (MQ=255)                                                                                                                            cgcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGAc                           <  4:295844/148‑1 (MQ=255)                                                                                                                            cgcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACg                          >  2:222604/1‑149 (MQ=255)                                                                                                                             gcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGa                                                                   <  1:62831/107‑1 (MQ=255)                                                                                                                             gcgTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGa                                                                   >  2:62831/1‑107 (MQ=255)                                                                                                                               gcaCAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGa                       >  3:401765/3‑149 (MQ=255)                                                                                                                                tACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGTCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGAt                      <  1:335882/149‑1 (MQ=255)                                                                                                                                   aaaaaGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGa                   >  3:88747/1‑149 (MQ=255)                                                                                                                                     aaaGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTGGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGaca                 <  3:126436/149‑1 (MQ=255)                                                                                                                                         gTCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcg             >  2:77626/1‑149 (MQ=255)                                                                                                                                         gTCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcg             >  4:4269/1‑149 (MQ=255)                                                                                                                                          tCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCcgc            >  4:113915/1‑149 (MQ=255)                                                                                                                                          tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                <  2:103994/129‑1 (MQ=255)                                                                                                                                          tCCTGTTCCAGGCGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCgg                                >  1:103994/1‑129 (MQ=255)                                                                                                                                            cTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGa          <  2:365206/149‑1 (MQ=255)                                                                                                                                              gTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAg         >  8:110996/1‑148 (MQ=255)                                                                                                                                                 ccGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAg     >  8:112361/1‑149 (MQ=255)                                                                                                                                                  cGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAg     <  3:45770/148‑1 (MQ=255)                                                                                                                                                    ccTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAg  <  4:401765/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGGGCATGGCTGCTGAATGAGCCGTTGTCGATGGCGCATCTATACGGCGGCCTGCTGATTGGCCTGGCCTTGTGGCTGGTGCTGCGCCCGGCGCGCAGTTGAAGCGGGGTACTGCATGAACAGCGCGTACAAAAAGGTCCTGTTCCGCCTGGAACAGGACGAGAACGGCTACCCACCGGCCTCGGTCGAGGGCCTGTGGGCCAAGGCCGTGGAAGGCGGCTATGCCATCGACAACATACCCTTCCATGTCTATGGCATCGCCCCGGGCGACGTGATCGACACCCGCGAGGAAGCAG  >  NC_002947/5187438‑5187733

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