Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,490,867	G→C	7.0%	P175R (CCG→CGG) ‡	clpV ←	protein ClpV1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,490,867	0	G	C	7.0%	173.0 / 4.0	57	P175R (CCG→CGG) ‡	clpV	protein ClpV1
Reads supporting (aligned to +/- strand):  ref base G (30/23);  new base C (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01

GCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGG  >  NC_002947/3490721‑3491013                                                                                                                                                   |                                                                                                                                                   gCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCggg                                                                                                                                                  >  4:76932/1‑149 (MQ=255)  ccGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGc                                                                                                                                                 >  8:138201/1‑149 (MQ=255)  ccGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCAATGGCGCCGCTGGGCTGGCCAGCGGGc                                                                                                                                                 >  2:291808/1‑149 (MQ=255)   cGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGCCAGCCCAGCggcgc                                                                                                                                                <  5:289487/149‑4 (MQ=255)    gTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               >  6:50003/1‑149 (MQ=255)    gTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTg                                                                                                                                               >  4:52495/1‑149 (MQ=255)       cTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcc                                                                                                                                            >  3:219488/1‑149 (MQ=255)       cTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcc                                                                                                                                            >  3:150587/1‑149 (MQ=255)       cTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgc                                                                                                                                             <  2:322816/148‑1 (MQ=255)           gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctgc                                                                                                                                         <  4:219488/148‑1 (MQ=255)           gAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCAGGCCAGCCCGCTGgccagcc                                                                                                                                        <  7:329605/149‑1 (MQ=255)                 aCGGCCCAGCACCGGGTCGATGCGGCCTTTGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACt                                                                                                                                  >  4:358287/1‑149 (MQ=255)                  cGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTc                                                                                                                                 <  3:76932/149‑1 (MQ=255)                  cGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTGGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTAGGCAGGCCAGCCCGCTGgccagcccagcggc                                                                                                                                 <  3:7598/149‑8 (MQ=255)                   ggCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCg                                                                                                                                >  1:47917/1‑149 (MQ=255)                       cAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGgcgc                                                                                                                            >  1:177414/1‑149 (MQ=255)                              gggTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctgcc                                                                                                                                        <  6:85839/130‑1 (MQ=255)                              gggTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGgcctgcc                                                                                                                                        >  5:85839/1‑130 (MQ=255)                               ggTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTcc                                                                                                                    <  3:52495/149‑1 (MQ=255)                                gTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGCCAGCCCAGCggcgc                                                                                                                                                <  6:115673/120‑4 (MQ=255)                                gTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGCCAGCCCAGCggcgc                                                                                                                                                >  5:115673/1‑117 (MQ=255)                                   gATGCGGCCTTCGCGGGCGCTTTGGGTACGGTTGACCGTGTACTGGTCCAGCGCCGGTGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCGGCCTGGACTCGGCGCTACCCTccccc                                                                                                                 <  2:316943/148‑1 (MQ=255)                                    aggCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAg                                                                                                               <  6:10232/147‑1 (MQ=255)                                      gCGGCCTTCGCGAGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAgcg                                                                                                             >  7:88204/1‑149 (MQ=255)                                                   gcgcTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAg                                                                                                      <  1:214399/143‑1 (MQ=255)                                                   gcgcTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAg                                                                                                      >  2:214399/1‑143 (MQ=255)                                                       tttGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGgcg                                                                                             >  5:260095/1‑148 (MQ=255)                                                          gggTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAgg                                                                                         >  6:26749/1‑149 (MQ=255)                                                            gTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGt                                                                                        >  1:82069/1‑148 (MQ=255)                                                             tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGcgcg                                                                                                           <  1:240464/128‑1 (MQ=255)                                                             tAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGcgcg                                                                                                           >  2:240464/1‑128 (MQ=255)                                                                ggTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCg                                                                                   <  4:150587/149‑1 (MQ=255)                                                                    gACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCg                                                                                <  8:299691/148‑1 (MQ=255)                                                                         tgtACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtt                                                                          <  5:26749/149‑1 (MQ=255)                                                                                                  gCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAgg                                                    <  2:9761/146‑1 (MQ=255)                                                                                                           cTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATcac                                             <  6:299583/144‑1 (MQ=255)                                                                                                           cTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATcac                                             >  5:299583/1‑144 (MQ=255)                                                                                                              gCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGTGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtt                                                                          >  2:354549/1‑112 (MQ=255)                                                                                                              gCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGTGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGAtt                                                                          <  1:354549/112‑1 (MQ=255)                                                                                                              gCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGAc                                     >  4:175273/1‑149 (MQ=255)                                                                                                               cTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACg                                    >  1:106722/1‑149 (MQ=255)                                                                                                                       gCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAg                            >  1:219063/1‑149 (MQ=255)                                                                                                                         tCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTATTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCa                          >  5:319788/1‑149 (MQ=255)                                                                                                                          cACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCAt                         >  2:355428/1‑149 (MQ=255)                                                                                                                           aCTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGATCAGGCGCAGGTCATCGGCGTTGATTATTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCa                          >  7:336703/1‑147 (MQ=255)                                                                                                                            cTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCg                       >  8:308524/1‑149 (MQ=255)                                                                                                                             tGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGt                      >  3:255475/1‑149 (MQ=255)                                                                                                                               gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTTCCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc                     <  1:291808/148‑1 (MQ=255)                                                                                                                               gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGt                      <  2:177414/147‑1 (MQ=255)                                                                                                                               gcgcCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc                     <  8:88204/148‑1 (MQ=255)                                                                                                                                cgccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTc                     >  6:14737/1‑147 (MQ=255)                                                                                                                                 gccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg                  >  4:162676/1‑149 (MQ=255)                                                                                                                                 gccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg                  <  3:162676/149‑1 (MQ=255)                                                                                                                                 gccgcTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAg                  <  2:47917/149‑1 (MQ=255)                                                                                                                                       aggcaggCCAGCCCGCTGGCCAGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGc            >  4:274028/2‑149 (MQ=255)                                                                                                                                        ggcaggCCAGCCCGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGcc           >  3:134713/1‑149 (MQ=255)                                                                                                                                          ctggCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTg         <  7:308524/149‑1 (MQ=255)                                                                                                                                             gCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTgcag      >  1:245126/1‑149 (MQ=255)                                                                                                                                               cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCa                             >  6:332628/1‑124 (MQ=255)                                                                                                                                               cAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCa                             <  5:332628/124‑1 (MQ=255)                                                                                                                                                 gCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAgg  <  3:175273/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCCGTACTTCGAACTCACGGCCCAGCACCGGGTCGATGCGGCCTTCGCGGGCGCTTTGGGTAAGGTTGACCGTGTACTGGTCCAGCGCCGGGGTCTTGCCGCTGGCCTTGCTCGACGTGCTCACTGGCGCCGCTGGGCTGGCCAGCGGGCTGGCCTGCCTGGACTCGGCGCTACCCTCCACCAGCGCGGCAAGGTTCAGGCGCAGGTCATCGGCGTTGATTTTTTCCAGCTCCGGTGCCGAGGCGATCACCACGCGACGCAGTTCAGCATCGTCAAGCAGAGCCTGCAGCAGG  >  NC_002947/3490721‑3491013

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