Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,749,873	T→A	6.1%	K206M (AAG→ATG)	PP_2406 ←	quinate/shikimate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,749,873	0	T	A	6.1%	199.9 / 2.6	65	K206M (AAG→ATG)	PP_2406	quinate/shikimate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base A (2/2);  total (35/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 = 9.30e-01

AACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCG  >  NC_002947/2749726‑2750019                                                                                                                                                    |                                                                                                                                                    aaCACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCtt                                                                                                                                                    <  8:43761/149‑1 (MQ=255)      cGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCa                                                                                                                                               >  3:45909/1‑149 (MQ=255)         cATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCATGGCCATGc                                                                                                                                            >  5:155468/1‑149 (MQ=255)         cATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCATGGCCATGc                                                                                                                                            >  5:149934/1‑149 (MQ=255)                  gcCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCg                                                                                                                                   >  5:211364/1‑149 (MQ=255)                   ccATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGc                                                                                                                                    <  7:67873/147‑1 (MQ=255)                          gCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGtt                                                                                                                            >  5:6809/1‑148 (MQ=255)                          gCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTg                                                                                                                           <  6:303647/149‑1 (MQ=255)                            gTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGGCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGAc                                                                                                                         <  2:14174/149‑1 (MQ=255)                                ggCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGc                                                                                                                     >  5:115613/1‑149 (MQ=255)                                        aGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGc                                                                                                             <  5:13509/149‑1 (MQ=255)                                           gCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTc                                                                                                          >  3:73616/1‑149 (MQ=255)                                            cccGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTc                                                                                                          <  4:45909/148‑1 (MQ=255)                                            cccGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCg                                                                                                         <  6:115613/149‑1 (MQ=255)                                             ccGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCgg                                                                                                        <  4:73616/149‑1 (MQ=255)                                              cGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTg                                                                                                                                                   <  7:332441/105‑1 (MQ=255)                                              cGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTg                                                                                                                                                   >  8:332441/1‑105 (MQ=255)                                               gggCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGc                                                                                                       >  4:194431/1‑148 (MQ=255)                                                 gCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCATGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCaa                                                                                                    <  6:149934/149‑1 (MQ=255)                                                 gCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCATGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCaa                                                                                                    <  6:155468/149‑1 (MQ=255)                                                    gcgcgcAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGc                                                                                                 <  6:6809/149‑1 (MQ=255)                                                      gcgcAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCtg                                                                                               >  3:325351/1‑149 (MQ=255)                                                                  ggTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCtgt                                                                                              <  1:87091/138‑1 (MQ=255)                                                                  ggTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCtgt                                                                                              >  2:87091/1‑138 (MQ=255)                                                                      tcCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGc                                                                               >  2:250475/1‑148 (MQ=255)                                                                                    gatgatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCt                                                                >  2:301734/1‑149 (MQ=255)                                                                                       gatCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTACGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACCCCGCGCTCcg                                                             >  6:50746/1‑149 (MQ=255)                                                                                         tctcTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGgcgcggagcgcg                                                           >  6:1485/1‑142 (MQ=255)                                                                                          ctctGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGC‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCgcgc                                                          >  6:178727/1‑149 (MQ=255)                                                                                           tctGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGGAGCGCGcc                                                         >  6:314876/1‑149 (MQ=255)                                                                                            ctGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGcc                                                         >  5:232958/1‑148 (MQ=255)                                                                                             tGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCcgcg                                                           <  6:211364/145‑1 (MQ=255)                                                                                               ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGGAGCGCGCCgtg                                                      <  5:314876/148‑3 (MQ=255)                                                                                               ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGGAGCGCGCCgt                                                       <  5:1485/147‑2 (MQ=255)                                                                                               ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGaaa                                                     >  5:233072/1‑149 (MQ=255)                                                                                               ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGaaa                                                     >  5:227785/1‑149 (MQ=255)                                                                                               ccacccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGCTAGCTTGGCCAAGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGGAGCGCGCCgt                                                       <  1:161157/147‑2 (MQ=255)                                                                                                  cccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGC‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATgg                                                  <  5:178727/149‑1 (MQ=255)                                                                                                   ccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCaa                                                                            <  5:224599/122‑1 (MQ=255)                                                                                                   ccacAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATggg                                                 >  2:23010/1‑149 (MQ=255)                                                                                                   caacaGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCaa                                                                            >  6:224599/3‑122 (MQ=255)                                                                                                    cacaGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGgt                                                <  1:23010/149‑1 (MQ=255)                                                                                                    cacaGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGgt                                                >  5:81930/1‑149 (MQ=255)                                                                                                    cacaGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGgt                                                <  3:194431/149‑1 (MQ=255)                                                                                                     acaGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGGCCCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGgt                                                <  1:250475/149‑1 (MQ=255)                                                                                                           ggggATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCTAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAgg                                         <  5:50746/149‑1 (MQ=255)                                                                                                               aTGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGt                                     <  4:325351/149‑1 (MQ=255)                                                                                                                tGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTc                                    >  7:329072/1‑149 (MQ=255)                                                                                                                tGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTc                                    <  2:342930/149‑1 (MQ=255)                                                                                                                   aacaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccac                                 <  6:227785/149‑1 (MQ=255)                                                                                                                   aacaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTccac                                 <  6:233072/149‑1 (MQ=255)                                                                                                                   aacaacTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCACGCCGAAATGGGTGTTCAGGTTGTccac                                 >  3:307856/1‑149 (MQ=255)                                                                                                                        cTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAgg                                         <  1:33790/136‑1 (MQ=255)                                                                                                                        cTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAgg                                         >  2:33790/1‑136 (MQ=255)                                                                                                                        cTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACg                            >  8:216270/1‑149 (MQ=255)                                                                                                                        cTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACg                            >  3:128211/1‑149 (MQ=255)                                                                                                                           cACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGTTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGt                                     <  6:216394/137‑1 (MQ=255)                                                                                                                           cACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGt                                     >  5:216394/1‑137 (MQ=255)                                                                                                                           cACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCt                         <  1:301734/149‑1 (MQ=255)                                                                                                                           cACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCt                         >  1:13758/1‑149 (MQ=255)                                                                                                                               ggCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTgcgc                     >  5:93783/1‑149 (MQ=255)                                                                                                                                      gCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGCCCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTc              >  4:28637/1‑149 (MQ=255)                                                                                                                                               gTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCtgt                                                                                              <  7:59234/61‑1 (MQ=255)                                                                                                                                               gTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCtgt                                                                                              >  8:59234/1‑61 (MQ=255)                                                                                                                                                   cTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCg  >  3:57932/1‑148 (MQ=255)                                                                                                                                                   cTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCg  >  7:141516/1‑148 (MQ=255)                                                                                                                                                    ttGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCAc      <  3:88325/143‑1 (MQ=255)                                                                                                                                                    ttGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCAc      >  4:88325/1‑143 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AACACCGCCATGTTGCTGCCATCCAGCGTGCGGCAGCCCAGTGCCCGGGCCGCGCGCAGCAGTTCGGTCTCCAGCGGGAAGTAGATGATCTCTGCCACCCACAGGCGGGGATGCAACAACTCCACCGGCAGTGGCGTGCCCGGTAGCTTGGCCATGCCCACCGGCGTGGTGTTGACCAGCCCGTCCGCTTCGGCCAAGGCTGTGG‑CCAGGTCGGTGCCAAGCACGGCGCGCTCCGCGCCGAAATGGGTGTTCAGGTTGTCCACCAACGCCTGCGCACGGGTCGCATCCACTTCG  >  NC_002947/2749726‑2750019

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑