Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,463,400	T→A	7.9%	intergenic (‑51/+48)	PP_4799 ← / ← lipA	muramoyltetrapeptide carboxypeptidase/lipoyl synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,463,400	0	T	A	7.9%	178.9 / 5.9	63	intergenic (‑51/+48)	PP_4799/lipA	muramoyltetrapeptide carboxypeptidase/lipoyl synthase
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base A (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCGCGCCGCGCCCGCCGGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCATCTTGTAA  >  NC_002947/5463257‑5463544                                                                                                                                                 |                                                                                                                                                 ccgcgccgcgCCCGCCGGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGgc                                                                                                                                             >  3:13196/1‑149 (MQ=255) ccgcgccgcgCCCGCCGGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGgc                                                                                                                                             <  5:191782/149‑1 (MQ=255)       cgcgccCGCCGGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAGATGCCGATGCGCGCCCAAGGgcggcat                                                                                                                                        <  6:181130/149‑5 (MQ=255)       cgcgccCGCCGGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCAt                                                                                                                                       <  4:175863/149‑1 (MQ=255)              gccgGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATc                                                                                                                                >  8:44466/1‑149 (MQ=255)              gccgGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATc                                                                                                                                >  4:251342/1‑149 (MQ=255)                    gCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGc                                                                                                                          >  4:150047/1‑149 (MQ=255)                           tGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGctc                                                                                                                   >  1:74976/1‑149 (MQ=255)                           tGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGctc                                                                                                                   >  2:79098/1‑149 (MQ=255)                             gCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTa                                                                                                                 <  7:86267/149‑1 (MQ=255)                                     tGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTagag                                                                                                         <  2:118806/149‑1 (MQ=255)                                        gTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCAt                                                                                                      <  5:50240/149‑1 (MQ=255)                                            gCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTGGGGCGCGCAGCGGCAGCTGTGGCTCCGCTCTAGGCTAGAGCATAGCt                                                                                                  >  5:14255/1‑149 (MQ=255)                                                   gcgcGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCt                                                                                           <  3:89011/149‑1 (MQ=255)                                                          tGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATa                                                                                                     <  8:119160/132‑1 (MQ=255)                                                          tGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATa                                                                                                     >  7:119160/1‑132 (MQ=255)                                                          tGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCtt                                                                                    >  8:33628/1‑149 (MQ=255)                                                           ggCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTg                                                                                   <  3:251342/149‑1 (MQ=255)                                                            gCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTg                                                                                   <  8:266693/148‑1 (MQ=255)                                                              aGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTc                                                                                                                        <  7:27482/109‑1 (MQ=255)                                                              aGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTc                                                                                                                        >  8:27482/1‑109 (MQ=255)                                                                ggtgggtTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTgg                                                                                                                                               >  7:58576/1‑84 (MQ=255)                                                                ggtgggtTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTgg                                                                                                                                               <  8:58576/84‑1 (MQ=255)                                                                 gtgggtTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCg                                                                             >  7:35251/1‑149 (MQ=255)                                                                    ggtTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTgg                                                                          <  7:33628/149‑1 (MQ=255)                                                                         gAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCctg                                                                     >  6:75997/1‑149 (MQ=255)                                                                          aaGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCctgc                                                                    <  8:150774/149‑1 (MQ=255)                                                                          aaGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCctgc                                                                    >  5:133330/1‑149 (MQ=255)                                                                                ttGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCtgt                                                                                                                             >  8:143428/1‑86 (MQ=255)                                                                                ttGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCtgt                                                                                                                             <  7:143428/86‑1 (MQ=255)                                                                                    gcgaGTACATCGCCGCTCCCTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGt                                                          <  1:208175/146‑1 (MQ=255)                                                                                    cgcAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGt                                                          <  1:164497/149‑1 (MQ=255)                                                                                        gTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTa                                                      <  5:75997/149‑1 (MQ=255)                                                                                        gTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTa                                                      <  8:35251/149‑1 (MQ=255)                                                                                        gTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTa                                                      <  4:13196/149‑1 (MQ=255)                                                                                          aCATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTAc                                                     >  2:163804/1‑148 (MQ=255)                                                                                           cATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACg                                                    >  5:138999/1‑148 (MQ=255)                                                                                                    tCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCAc                                          >  8:59531/1‑149 (MQ=255)                                                                                                           aGCGCATTTGACACAACCAAAAAGGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCGCAGAGCTTGAGCTTGGCTTCGTGGGCCTGCTGGGCGGCGGGGTACGAAGAACGCACCAGCggg                                   <  6:14255/149‑1 (MQ=255)                                                                                                             cgcATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGc                                  >  7:229489/1‑148 (MQ=255)                                                                                                                 ttGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGaa                             <  7:44466/149‑1 (MQ=255)                                                                                                                    gACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGcc                                 >  1:295440/1‑142 (MQ=255)                                                                                                                    gACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGcc                                 <  2:295440/142‑1 (MQ=255)                                                                                                                         taaaaaaaaTGGCGATTC‑CGCCCTTGGGCGCGCATGGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTTGTa                                                      <  1:288038/112‑1 (MQ=255)                                                                                                                           acaaaAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCt                   >  6:61584/1‑149 (MQ=255)                                                                                                                           acaaaAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCt                   >  7:20773/1‑149 (MQ=255)                                                                                                                           acaaaAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCt                   <  2:283444/149‑1 (MQ=255)                                                                                                                            caaaAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCtt                  <  1:79098/149‑1 (MQ=255)                                                                                                                             aaaaaTGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCtt                  >  7:289933/1‑148 (MQ=255)                                                                                                                               aaaTGCCGATGC‑CGCCCTTGTGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGaa               <  8:213745/149‑1 (MQ=255)                                                                                                                                    ccgatgc‑cgCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGCCATCATGGGCGGCGTGGTAGGAACAACGCACCAGCGGGCAAGAGGCAACGTGGTTGAGGccc           >  2:107460/1‑148 (MQ=255)                                                                                                                                     cgatgcgcGCCCAAGGGCG‑GCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGaa                             <  7:173178/129‑1 (MQ=255)                                                                                                                                     cgatgcgcGCCCAAGGGCG‑GCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGaa                             >  8:173178/1‑129 (MQ=255)                                                                                                                                     cgatgc‑cgCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCa          >  6:180659/1‑148 (MQ=255)                                                                                                                                     cgatgc‑cgCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCAt         >  2:11047/1‑149 (MQ=255)                                                                                                                                     agatgc‑cgCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCa          >  1:316042/2‑148 (MQ=255)                                                                                                                                      gatgc‑cgCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCATc        <  8:229489/149‑1 (MQ=255)                                                                                                                                         gcgcgccCAAGGGCG‑GCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGTCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGccc           >  3:97303/4‑143 (MQ=255)                                                                                                                                         gcgcgccCAAGGGCG‑GCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGTCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGccc           <  4:97303/140‑1 (MQ=255)                                                                                                                                            cgccCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGaa                             <  6:255444/123‑1 (MQ=255)                                                                                                                                            cgccCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGaa                             >  5:255444/1‑123 (MQ=255)                                                                                                                                            cgccCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCATCTTGt    >  2:226115/1‑148 (MQ=255)                                                                                                                                             gccCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCATCTTGTaa  >  3:166328/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                 CCGCGCCGCGCCCGCCGGGGCGACAATGGCAAAACATGCGTTGGCCGGCAGCGCGGCTGGCAGGGTGGGTTGGAAGGTTTTGGCGCAGTACATCGCGGCTCCTTGCAGCGCATTGGACACAAACAAAAATGCCGATGC‑CGCCCTTGGGCGCGCATCGGCATCTGTGGCTCCGCTCTAGGCTAGAGCATAGCTCAGAGCTTGATCTTGGCTTCGTGGGCCTGCTGGTCGGCGTGGTACGAAGAACGCACCAGCGGGCCAGAAGCAACGTTCTTGAAGCCCATCTTGTAA  >  NC_002947/5463257‑5463544

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