Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,238,149	C→A	6.6%	S196S (TCC→TCA)	PP_2836 →	2‑keto‑3‑deoxyxylonate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,238,149	0	C	A	6.6%	190.0 / 3.8	61	S196S (TCC→TCA)	PP_2836	2‑keto‑3‑deoxyxylonate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (31/26);  new base A (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

GACCTGTCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGACTTCGAAGG  >  NC_002947/3238009‑3238295                                                                                                                                             |                                                                                                                                                    gACCTGTCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTgg                                                                                                                                             <  7:303984/149‑1 (MQ=255) gACCTGTCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCGGACAGCGGTTgg                                                                                                                                             <  3:295565/149‑1 (MQ=255)  aCCTGTCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGc                                                                                                                                            >  2:265683/1‑149 (MQ=255)    cTGTCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCa                                                                                                                                          >  5:298744/1‑149 (MQ=255)       tCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCAc                                                                                                                                         <  5:55230/147‑1 (MQ=255)         gCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGgcg                                                                                                                                     >  2:47876/1‑149 (MQ=255)          cAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGgcgc                                                                                                                                    <  3:270371/149‑1 (MQ=255)           aGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGcc                                                                                                                                   <  4:145753/149‑1 (MQ=255)               cGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACa                                                                                                                               >  8:137834/1‑149 (MQ=255)                gTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACAt                                                                                                                              >  2:203645/1‑149 (MQ=255)                gTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACAt                                                                                                                              <  6:23868/149‑1 (MQ=255)                    ccGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCgg                                                                                                                           >  7:41873/1‑148 (MQ=255)                           ggctcAGGCCGAAGTGCTGCGCTAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCATTAGCTGTAGGTGGGTATCGGCCCGGATTCCGATGTGTTCACCATGGCCCAACAGCTGTGCGCGGTTGGGCTCGGCGCCGACATCGG‑CATCCAc                                                                                                                   <  2:264271/148‑1 (MQ=255)                             ctcAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCAcccc                                                                                                                <  7:137834/149‑1 (MQ=255)                              tcAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCAcccc                                                                                                                >  6:113660/1‑148 (MQ=255)                              tcAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCGGACAGCGGTTgg                                                                                                                                             >  5:50613/1‑120 (MQ=255)                              tcAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCGGACAGCGGTTgg                                                                                                                                             <  6:50613/120‑1 (MQ=255)                                  gCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTc                                                                                                           >  3:246314/1‑149 (MQ=255)                                    cGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCACAAGTc                                                                                                           >  8:317017/1‑147 (MQ=255)                                     gAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAg                                                                                                        >  7:246401/1‑149 (MQ=255)                                      aaGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGc                                                                                                       <  5:87241/149‑1 (MQ=255)                                              cgcAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGaacaa                                                                                               <  1:47876/149‑1 (MQ=255)                                              cgcAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGaaca                                                                                                <  1:203645/148‑1 (MQ=255)                                              cgcAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGCCCAGCTGGaacaa                                                                                               <  6:276464/149‑1 (MQ=255)                                                    gTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCaga                                                                                         <  6:298744/149‑1 (MQ=255)                                                         ggTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTg                                                                                                                                              >  4:126589/1‑92 (MQ=255)                                                         ggTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTg                                                                                                                                              <  3:126589/92‑1 (MQ=255)                                                          gTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCgg                                                                                    >  7:1304/1‑148 (MQ=255)                                                               gCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAgtgg                                                                              >  2:205493/1‑149 (MQ=255)                                                                cAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGaa                                                                                  >  6:222878/1‑144 (MQ=255)                                                                cAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGaa                                                                                  <  5:222878/144‑1 (MQ=255)                                                                  gggCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGc                                                                           <  4:246314/149‑1 (MQ=255)                                                                         tgGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCgt                                                                    <  1:205493/149‑1 (MQ=255)                                                                           gTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCgtgt                                                                  >  6:133956/1‑149 (MQ=255)                                                                               cAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCCCCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAg                                                              >  5:218522/1‑149 (MQ=255)                                                                                        gAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACAGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACCACCCAGAGCCGGAAGTGGTGATGGCCGTGTCGAGCAATGgcag                                                     >  6:351388/1‑149 (MQ=255)                                                                                            tGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATc                                                 <  7:317017/149‑1 (MQ=255)                                                                                              ggTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTg                                                                         <  6:156079/123‑1 (MQ=255)                                                                                              ggTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTg                                                                         >  5:156079/1‑123 (MQ=255)                                                                                              ggTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCaa                                               <  5:178185/149‑1 (MQ=255)                                                                                                     gCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCgtgt                                                                  <  8:176102/123‑1 (MQ=255)                                                                                                     gCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCgtgt                                                                  >  7:176102/1‑123 (MQ=255)                                                                                                     gCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGgcgc                                         >  1:68145/1‑148 (MQ=255)                                                                                                     gCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGgcgc                                         >  5:72612/1‑148 (MQ=255)                                                                                                            tGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCg                                 <  8:246401/149‑1 (MQ=255)                                                                                                               cGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATgg                                                        >  4:16703/1‑123 (MQ=255)                                                                                                               cGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATgg                                                        <  3:16703/123‑1 (MQ=255)                                                                                                                      ttCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGACGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGt                       <  1:265683/149‑1 (MQ=255)                                                                                                                         aCCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCGGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCaa                    >  7:72334/1‑149 (MQ=255)                                                                                                                         aCCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCGGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCaa                    >  7:67483/1‑149 (MQ=255)                                                                                                                         aCCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCaa                    >  3:57565/1‑149 (MQ=255)                                                                                                                           cAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAAcc                  >  7:22212/1‑149 (MQ=255)                                                                                                                              ggCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTgc               <  6:218522/149‑1 (MQ=255)                                                                                                                                 ccAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCaa                                               >  7:74193/1‑114 (MQ=255)                                                                                                                                 ccAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCaa                                               <  8:74193/114‑1 (MQ=255)                                                                                                                                  cAACCGCGGACAGCGGTTGGCCACGGCGCCGACATCGGCCATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGACGAGCAATGGCAGCATCAAggcc                                           >  1:240037/1‑116 (MQ=255)                                                                                                                                   aaCCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGTCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGAc          >  7:91485/1‑149 (MQ=255)                                                                                                                                        cTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGACTACgc     >  2:340452/1‑148 (MQ=255)                                                                                                                                          gTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGTTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGACTTCGAAg   >  1:341366/1‑149 (MQ=255)                                                                                                                                           tCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGACTTCGaa    <  5:193578/147‑1 (MQ=255)                                                                                                                                           tCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGACTTCGAAgg  >  7:197693/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                    GACCTGTCGCAGATCGTGCCCGGCTCGGCTCAGGCCGAAGCGCTGCGCAAGGTGCTGGTAGAGCAGGGCCTGTGGTCGCAGTACCTGGAGGTGGGTATCGGCCCGGATGCCGAAGTGTTCACCAAGGCCCAACCGCTGTCCGCGGTTGGCCACGGCGCCGACATCGG‑CATCCACCCCAAGTCCAGCTGGAACAACCCAGAGCCGGAAGTGGTGCTGGCCGTGTCGAGCAATGGCAGCATCAAGGGCGCCATGCTCGGCAACGACGTCAACCTGCGCGACTTCGAAGG  >  NC_002947/3238009‑3238295

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