Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,866	T→C	5.3%	Y37H (TAC→CAC)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,866	0	T	C	5.3%	222.5 / 2.2	75	Y37H (TAC→CAC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/38);  new base C (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GTGTAAGCCTGACCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGG  >  NC_002947/3762718‑3763010                                                                                                                                                     |                                                                                                                                                   gtgtAAGCCTGACCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCgtgt                                                                                                                                                    <  8:2248/149‑1 (MQ=255)  ggtAAGCCTGACCCTGAGCCTTTGAACCGAGGAACCATCATGCATGCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCGCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTa                                                                                                                                                   <  4:75405/148‑1 (MQ=255)     aaGCCTGACCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACgg                                                                                                                                                >  7:290832/1‑149 (MQ=255)            aCCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACgg                                                                                                                                          <  4:132987/148‑1 (MQ=255)            aCCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGgc                                                                                                                                         >  3:195374/1‑149 (MQ=255)            aCCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTACACGGCCTCga                                                                                                                                          <  5:70386/148‑2 (MQ=255)               cTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGcc                                                                                                                                      >  5:101007/1‑149 (MQ=255)               cTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTACACGGCCTCgat                                                                                                                                         <  6:228742/146‑3 (MQ=255)               cTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTACACGGCCTCgat                                                                                                                                         >  5:228742/1‑144 (MQ=255)                tGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGcc                                                                                                                                      >  2:241160/1‑148 (MQ=255)                 gAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGcc                                                                                                                                      >  8:3884/1‑147 (MQ=255)                      tttGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGAt                                                                                                                                >  5:342236/1‑148 (MQ=255)                       ttGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATaa                                                                                                                              >  1:308663/1‑149 (MQ=255)                        tGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGCCCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAAc                                                                                                                             >  6:280056/1‑149 (MQ=255)                            ccGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGcc                                                                                                                          >  2:229105/1‑148 (MQ=255)                             cGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTa                                                                                                                        >  7:166479/1‑149 (MQ=255)                              gAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTa                                                                                                                        >  4:317455/1‑148 (MQ=255)                               aGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTAc                                                                                                                       >  8:182239/1‑148 (MQ=255)                               aGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTAc                                                                                                                       >  1:114507/1‑148 (MQ=255)                                     catcatGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCgg                                                                                                               <  5:170679/149‑1 (MQ=255)                                     catcatGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCg                                                                                                                >  2:132275/1‑148 (MQ=255)                                     catcatGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCg                                                                                                                >  1:120121/1‑148 (MQ=255)                                         atGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGaa                                                                                                             >  8:100931/1‑147 (MQ=255)                                            catccatcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGa                                                                                                         >  8:167894/1‑148 (MQ=255)                                             atccatcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATc                                                                                                       <  5:44157/149‑1 (MQ=255)                                             atccatcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATc                                                                                                       <  4:195374/149‑1 (MQ=255)                                             atccatcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATc                                                                                                       <  7:182239/149‑1 (MQ=255)                                                catcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCaa                                                                                                     <  5:11980/148‑1 (MQ=255)                                                catcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCaa                                                                                                     <  7:167894/148‑1 (MQ=255)                                                 atcGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAgc                                                                                                   <  8:229882/149‑1 (MQ=255)                                                     cTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAgcgcgc                                                                                               <  3:317455/149‑1 (MQ=255)                                                         aaGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGTGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTgg                                                                                           <  2:308663/149‑1 (MQ=255)                                                         aaGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTgg                                                                                           <  8:290832/149‑1 (MQ=255)                                                           gTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGc                                                                                         <  6:101007/149‑1 (MQ=255)                                                               cTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                     <  1:11054/149‑1 (MQ=255)                                                                ttCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTc                                                                                     >  1:243540/1‑148 (MQ=255)                                                                  cccTGTTGGCCCTGATGGTCAGTGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑c                                                                                 >  4:163378/1‑149 (MQ=255)                                                                  cccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACGGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGa                                                                                   <  1:282051/148‑1 (MQ=255)                                                                  cccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCg                                                                                    <  2:239928/147‑1 (MQ=255)                                                                   ccTGTTGGCCCTGATGGTCAGCGTCAACGCGGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑Cg                                                                                >  5:48784/1‑149 (MQ=255)                                                                   ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGaaaga                                                                                                          <  6:187487/124‑3 (MQ=255)                                                                   ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGaaaga                                                                                                          >  5:187487/1‑122 (MQ=255)                                                                   ccTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑c                                                                                 <  2:120121/148‑1 (MQ=255)                                                                    cTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGt                                                                               >  5:126100/1‑149 (MQ=255)                                                                       ttGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGa                                                                             <  6:342236/148‑1 (MQ=255)                                                                           cccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTgg                                                                                           <  5:128326/131‑1 (MQ=255)                                                                           cccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTgg                                                                                           >  6:128326/1‑131 (MQ=255)                                                                              tGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGaa                                                                      >  3:126767/1‑148 (MQ=255)                                                                                  ggTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATg                                                                  >  1:251244/1‑149 (MQ=255)                                                                                    tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                               <  8:166479/149‑1 (MQ=255)                                                                                    tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                               <  5:199009/149‑1 (MQ=255)                                                                                    tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                               >  6:332851/1‑149 (MQ=255)                                                                                    tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                               <  6:126100/149‑1 (MQ=255)                                                                                        cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCgg                                                           <  7:3884/149‑1 (MQ=255)                                                                                            aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCa                                                        <  1:241160/148‑1 (MQ=38)                                                                                                  ggCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGa                                                 <  2:114507/149‑1 (MQ=255)                                                                                                     gcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCACgtcgaacg                                               >  4:75965/1‑142 (MQ=38)                                                                                                      cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAacgc                                             <  2:251244/149‑1 (MQ=38)                                                                                                          cGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑GTCGAACGcccag                                         <  3:75965/149‑4 (MQ=37)                                                                                                           ggCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGc                                        <  2:203370/149‑1 (MQ=255)                                                                                                              cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                     <  4:126767/149‑1 (MQ=255)                                                                                                                 tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGAt                                    <  5:280056/147‑1 (MQ=255)                                                                                                                   gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                <  5:332851/149‑1 (MQ=255)                                                                                                                   gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                >  7:34612/1‑149 (MQ=255)                                                                                                                        aggagCGACGTGGCAATATCGAGGCCGTGGACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                           <  3:280017/149‑1 (MQ=255)                                                                                                                         ggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                           <  1:295843/148‑1 (MQ=255)                                                                                                                          gatcgtcGTTGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGTGAACATGATCCCGGTCACGGTCGAGCGCGTGGCGGATCATGCCGtgct                         <  5:192089/143‑1 (MQ=255)                                                                                                                          gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgct                         >  2:253640/1‑149 (MQ=255)                                                                                                                          gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgct                         <  6:48784/149‑1 (MQ=255)                                                                                                                           agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgctg                        <  1:253640/149‑1 (MQ=255)                                                                                                                           agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                           >  3:81548/1‑146 (MQ=255)                                                                                                                           agCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                           <  4:81548/146‑1 (MQ=255)                                                                                                                                cGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAcc                   <  1:132275/149‑1 (MQ=255)                                                                                                                                     cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAGCATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCACCCACTTa              >  1:153195/1‑149 (MQ=255)                                                                                                                                          tCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCg         >  6:2324/1‑149 (MQ=255)                                                                                                                                           cgtggCCGTACACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc        >  8:136470/4‑149 (MQ=255)                                                                                                                                                  gtgtACGGCCACGGCGCCGAGGATAACGCCTACAGAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGATCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg  >  1:47095/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GTGTAAGCCTGACCCTGAGCCTTTGAACCGAGGAACCATCATGCATCCATCGCTGAAAGTCGCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACA‑AGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGG  >  NC_002947/3762718‑3763010

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