Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,883,595	C→G	11.3%	A60G (GCC→GGC)	PP_5157 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,883,595	0	C	G	11.3%	174.7 / 12.5	62	A60G (GCC→GGC)	PP_5157	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base G (4/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

ATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAA  >  NC_002947/5883453‑5883738                                                                                                                                                |                                                                                                                                                 aTCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTgg                                                                                                                                             <  7:294841/148‑1 (MQ=255)        cATTGGCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGaa                                                                                                                                     >  8:140360/1‑149 (MQ=255)         aTTGGCTCAGGCGGAGCTGATCGACGAAATCGCGCGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGaa                                                                                                                                     >  1:229049/1‑149 (MQ=255)         aTTGGCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaaattcg                                                                                                                                    <  7:336154/149‑6 (MQ=255)            ggCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGGCCTGGCCaaattcgat                                                                                                                                  <  5:102657/148‑8 (MQ=255)             gCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGa                                                                                                                                <  4:136379/149‑1 (MQ=255)                  ggCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGCCATCGAATTTGGCCCGGCCCTGGCCACGGAACTGGATTTgc                                                                                                                           >  4:111218/1‑149 (MQ=255)                        gCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCg                                                                                                                     <  6:116254/149‑1 (MQ=255)                           gATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGt                                                                                                                  >  5:8891/1‑149 (MQ=255)                               gacgaAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgt                                                                                                              >  6:6789/1‑149 (MQ=255)                               gacgaAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcg                                                                                                               >  6:186448/1‑148 (MQ=255)                               gacgaAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcg                                                                                                               <  5:186448/148‑1 (MQ=255)                                acgaAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGTAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtc                                                                                                             >  2:84963/1‑149 (MQ=255)                                acgaAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtc                                                                                                             >  6:117772/1‑149 (MQ=255)                                  gaAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTATAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTtcgtcg                                                                                                            >  7:185528/1‑148 (MQ=255)                                    aaTCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAg                                                                                                         <  7:140360/149‑1 (MQ=255)                                               cGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAc                                                                                               <  5:117772/149‑1 (MQ=255)                                               cGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAc                                                                                               >  3:91956/1‑149 (MQ=255)                                                   gCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGaa                                                                                           <  6:8891/149‑1 (MQ=255)                                                             cATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                 <  1:117922/149‑1 (MQ=255)                                                             cATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCgg                                                                                 >  2:117922/1‑149 (MQ=255)                                                              aTTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCATGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCggg                                                                                >  8:28002/1‑149 (MQ=255)                                                                  ccGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTc                                                                            >  6:358791/1‑149 (MQ=255)                                                                   cGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGt                                                                                                                  >  3:139366/1‑110 (MQ=255)                                                                   cGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGt                                                                                                                  <  4:139366/110‑1 (MQ=255)                                                                    gTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCga                                                                          <  2:187625/149‑1 (MQ=255)                                                                       cAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga                                                                        <  3:111218/148‑1 (MQ=255)                                                                         ggCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                     <  8:185528/149‑1 (MQ=255)                                                                         ggCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTg                                                                     <  5:6789/149‑1 (MQ=255)                                                                          gCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTgg                                                                    >  5:348184/1‑149 (MQ=255)                                                                               gaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga                                                                        >  6:284430/1‑140 (MQ=255)                                                                               gaACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCgaga                                                                        <  5:284430/140‑1 (MQ=255)                                                                                                   agcaagACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCg                                          >  7:339204/1‑149 (MQ=255)                                                                                                       agACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTc                                        >  6:127425/1‑147 (MQ=255)                                                                                                         acgacCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCATGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGt                                    <  7:28002/149‑1 (MQ=255)                                                                                                         acgacCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGt                                    <  8:339204/149‑1 (MQ=255)                                                                                                           gacCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg                                                                                                   <  4:28027/85‑1 (MQ=255)                                                                                                           gacCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCgg                                                                                                   >  3:28027/1‑85 (MQ=255)                                                                                                              caccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATc                               >  1:47180/1‑149 (MQ=255)                                                                                                              caccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATc                               <  4:283675/149‑1 (MQ=255)                                                                                                                ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAc                                                                                               >  6:108188/1‑84 (MQ=255)                                                                                                                ccTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTAc                                                                                               <  5:108188/84‑1 (MQ=255)                                                                                                                        gTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGc                     >  4:334081/1‑149 (MQ=255)                                                                                                                         tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGc                                                                                   >  3:193354/1‑87 (MQ=255)                                                                                                                         tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGc                                                                                   <  4:193354/87‑1 (MQ=255)                                                                                                                         tttCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGc                     >  1:157351/1‑148 (MQ=255)                                                                                                                             gACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGATGACGCGGCTACCGAAGTGATGGCGGGTGTCGAGAGTGGCA‑AGTACGAAATCGCATTCACCACGTCGAGTGAATCGATCACTGGCGaagt                >  6:341404/1‑146 (MQ=255)                                                                                                                              aCATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCACAGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACgg                <  1:84963/149‑1 (MQ=255)                                                                                                                               cATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGcc              >  2:83881/1‑149 (MQ=255)                                                                                                                                cagttccttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCCCCCCGTCGAGTGAATCGCTCACAGGCGACGGccc             >  1:316905/8‑149 (MQ=255)                                                                                                                                aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGt                                    >  2:157990/1‑126 (MQ=255)                                                                                                                                aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGt                                    <  1:157990/126‑1 (MQ=255)                                                                                                                                aTCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGccc             >  6:71318/1‑149 (MQ=255)                                                                                                                                   ttccttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCt          >  6:235428/5‑149 (MQ=255)                                                                                                                                   gAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGc           <  4:91956/148‑1 (MQ=255)                                                                                                                                   gAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGAGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCt          >  7:66094/1‑149 (MQ=255)                                                                                                                                    tccttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTc         >  3:281111/4‑149 (MQ=255)                                                                                                                                      cttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGa       >  7:117338/2‑149 (MQ=255)                                                                                                                                      cttGGCCAGGGCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGa       <  5:235428/148‑1 (MQ=255)                                                                                                                                          gCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGt    <  6:348184/148‑1 (MQ=255)                                                                                                                                           ccAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTaa  <  3:334081/149‑1 (MQ=255)                                                                                                                                           ccAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTaa  <  2:47180/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 ATCAGTCCATTGGCTCAGGCGGAGCTGATCGACGAAATCGC‑CGACCGGGGCGAACTGCGCATTGCCGTACAGGCCGAGAACTCGCCATACGCCTTCAAGCAAGACGACCACCTGACAGGTTTCGACATCGAATTTGGCCAGGCCCTGGCCAAGGAACTGGATTTGCGCGCCGAGTTCGTCGAAGCCCCGGCTACCGAAGTGCTGGCGGGTGTCGAGAGTGGCA‑AGTACGACATCGCATTCACCCCGTCGAGTGAATCGCTCACAGGCGACGGCCCGCTCGATGTAA  >  NC_002947/5883453‑5883738

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