Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,108,664	T→C	6.4%	I134T (ATT→ACT)	PP_4494 →	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,108,664	0	T	C	6.4%	178.8 / 2.8	62	I134T (ATT→ACT)	PP_4494	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCGGCAGCCAAGGCTG  >  NC_002947/5108517‑5108811                                                                                                                                                    |                                                                                                                                                     ccACCTTGCAGGCACTGGCCTATCAATCGGGATTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGCCATGCTGGTGCACATTCATTCACGCAtt                                                                                                                                                     <  1:180896/149‑1 (MQ=255)    ccTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGaa                                                                                                                                                  <  3:253482/149‑1 (MQ=255)     cTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGaaa                                                                                                                                                 >  8:49773/1‑149 (MQ=255)     cTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGaaa                                                                                                                                                 <  3:91344/149‑1 (MQ=255)        gCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAAc                                                                                                                                                <  8:71046/147‑1 (MQ=255)        gCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAAc                                                                                                                                                <  8:18628/147‑1 (MQ=255)         cAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCAt                                                                                                                                             <  6:134763/149‑1 (MQ=255)          aGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCAt                                                                                                                                             >  7:117669/1‑148 (MQ=255)            gCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCaa                                                                                                                                          <  8:117669/149‑1 (MQ=255)              aCTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACt                                                                                                                                        >  8:15963/1‑149 (MQ=255)               cTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACtt                                                                                                                                       <  4:310382/149‑1 (MQ=255)                      aTCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGacac                                                                                                                                 >  7:191631/1‑148 (MQ=255)                      aTCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGacac                                                                                                                                 >  7:95915/1‑148 (MQ=255)                         aaTCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGc                                                                                                                             >  3:183601/1‑149 (MQ=255)                          aTCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGc                                                                                                                             >  4:183970/1‑148 (MQ=255)                             gggAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGaa                                                                                                                          <  8:21907/148‑1 (MQ=255)                               gAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGacac                                                                                                                                 <  7:189792/139‑1 (MQ=255)                               gAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGacac                                                                                                                                 >  8:189792/1‑139 (MQ=255)                                  tAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGGGCACATTCATTCACGCATTGAACCCATCACCTGCGACACCAGCGAAATCgct                                                                                                                    >  3:181732/1‑147 (MQ=255)                                     gCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATc                                                                                                                                            <  5:126166/122‑1 (MQ=255)                                     gCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATc                                                                                                                                            >  6:126166/1‑122 (MQ=255)                                     gCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATgc‑c                                                                                                                <  5:193188/149‑1 (MQ=255)                                      cACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATgc‑cg                                                                                                               <  3:183970/149‑1 (MQ=255)                                       aCGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATgc‑cgc                                                                                                              >  3:204861/1‑149 (MQ=255)                                                     cACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGgc                                                                                                  >  4:288574/1‑147 (MQ=255)                                                       cccTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGgccagc                                                                                              <  1:83749/149‑1 (MQ=255)                                                        ccTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGgccagc                                                                                              >  5:38942/1‑148 (MQ=255)                                                            ccAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACg                                                                                          <  8:191631/148‑1 (MQ=255)                                                              aaCCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCTGGCCACgccg                                                                                       <  7:49773/149‑4 (MQ=255)                                                                 cTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGc                                                                                    >  3:191014/1‑149 (MQ=255)                                                                      gTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGt                                                                               >  4:280786/1‑149 (MQ=255)                                                                          gAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCgg                                                                                 <  4:165507/143‑1 (MQ=255)                                                                          gAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCgg                                                                                 >  3:165507/1‑143 (MQ=255)                                                                          gAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTg                                                                           >  1:227126/1‑149 (MQ=255)                                                                              tGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTg                                                                       <  6:38942/149‑1 (MQ=255)                                                                                  tGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCAt                                                                   >  2:249316/1‑149 (MQ=255)                                                                                               gCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGcc                                                       >  2:182541/1‑148 (MQ=255)                                                                                                 aTGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACAAAAGCGAAATCGATGC‑CGCCATCGGCGTGGCCCGCCACGACTTGCCGGGTCTTGTCTGACAGCGCCTGCATGCCgag                                                    >  8:253765/1‑149 (MQ=255)                                                                                                   gCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCgag                                                    >  1:5637/1‑147 (MQ=255)                                                                                                     gTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCACTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGc                                                >  3:300058/1‑149 (MQ=255)                                                                                                     gTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCACTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGc                                                >  1:4812/1‑149 (MQ=255)                                                                                                           caccgcTCACCCCGTCATTCTGGTGCACATTCAGTCCCGCATTGAAACCATCAACTACTCCACCAGCGAACTCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCtggtg                                           <  4:181732/145‑1 (MQ=255)                                                                                                            aGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTggtggt                                         >  3:114590/1‑149 (MQ=255)                                                                                                                   aCCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTg                                  >  1:102690/1‑149 (MQ=255)                                                                                                                   aCCCGGGCATGCTGGTGCACATTCATTCACGCACTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCtt                                   <  4:300058/148‑1 (MQ=255)                                                                                                                   aCCCGGGCATGCTGGTGCACATTCATTCACGCACTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCtt                                   <  2:4812/148‑1 (MQ=255)                                                                                                                     ccGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACtt                                                                                                                                       >  4:253835/1‑47 (MQ=255)                                                                                                                     ccGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACtt                                                                                                                                       <  3:253835/47‑1 (MQ=255)                                                                                                                     ccGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTgcc                                >  4:124671/1‑149 (MQ=255)                                                                                                                         gCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGccgccg                            <  2:5637/149‑1 (MQ=255)                                                                                                                            tGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGGGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAg                         <  4:204861/149‑1 (MQ=255)                                                                                                                                 gTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTg                     <  2:227126/148‑1 (MQ=255)                                                                                                                                    cacaTTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCg                  <  3:124671/148‑1 (MQ=255)                                                                                                                                     acaTTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTc                <  2:102690/149‑1 (MQ=255)                                                                                                                                      cattcattCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCg               <  3:280786/149‑1 (MQ=255)                                                                                                                                       attcattcACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCtggtgg                                          <  7:136121/121‑1 (MQ=255)                                                                                                                                       attcattcACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCtggtgg                                          >  8:136121/1‑121 (MQ=255)                                                                                                                                       attcattcACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCgg              >  1:247309/1‑149 (MQ=255)                                                                                                                                       attcattcACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCgg              >  7:212335/1‑149 (MQ=255)                                                                                                                                               aCGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑TGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCGGCAGCCAAg      >  7:63204/1‑149 (MQ=255)                                                                                                                                               aCGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCGGCAGCCAAg      >  5:251434/1‑149 (MQ=255)                                                                                                                                                  cATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCTTGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCGGCAGCCAAGGCt   <  1:249316/149‑1 (MQ=255)                                                                                                                                                   aTTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCGGCAGCCAAGGCTg  <  4:191014/149‑1 (MQ=255)                                                                                                                                                    ttGAAACCATCAACTTCGACACCAGCGAAATCGATGCGCGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTc                >  1:57927/1‑135 (MQ=255)                                                                                                                                                    ttGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTATTGTCTGCGATCGCCTGCTTGCGGAGGAGCTGGTGGTGATCTTGCCGGAGGAGGCTGGCGTc                <  2:57927/134‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCACCTTGCAGGCACTGGCCTATCAATCGGGAGTAGGCACGTTGCGCCTGGCCACCCTGCCAACCTTCGGTTCGAAATGGCTGTTGCCGAGGCTGCATGCGTTCTACAGCGCTCACCCGGGCATGCTGGTGCACATTCATTCACGCATTGAAACCATCAACTTCGACACCAGCGAAATCGATGC‑CGCCATCGGCGTGGCCAGCCACGACTTGCCGGGTCTTGTCTGCCATCGCCTGCATGCCGAGGAGCTGGTGGTGATCTTGCCGCCGGAGGCTGGCGTCGGCAGCCAAGGCTG  >  NC_002947/5108517‑5108811

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