Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,625,095	C→A	6.0%	P222P (CCG→CCT)	PP_4943 ←	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,625,095	0	C	A	6.0%	211.7 / 2.7	67	P222P (CCG→CCT)	PP_4943	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base A (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GAACGCCTGGAAATAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGGCTGCTTGCAGGGC  >  NC_002947/5624948‑5625241                                                                                                                                                    |                                                                                                                                                   gAACGCCTGGATATAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCg                                                                                                                                                   >  4:118064/1‑149 (MQ=255) gAACGCCTGGAAATAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCg                                                                                                                                                   >  1:126204/1‑149 (MQ=255)       cTGGAAATAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTg                                                                                                                                             >  6:299964/1‑149 (MQ=255)       cTGGAAATAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTg                                                                                                                                             >  4:208955/1‑149 (MQ=255)           aaaTAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTg                                                                                                                                         >  8:53948/1‑149 (MQ=255)            aaTAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTgg                                                                                                                                        <  4:142029/149‑1 (MQ=255)             ttAACGGCGCGCATCCGGTCCCTGGCCAAGAAAGTCGTCTTGCCCGGCATTTCCAATTTCGTCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTggg                                                                                                                                       <  4:193494/148‑1 (MQ=255)             aTAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTggg                                                                                                                                       <  3:124526/149‑1 (MQ=255)                  ggCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCaaaa                                                                                                                                  >  2:124972/1‑149 (MQ=255)                     gcgcATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACcg                                                                                                                               >  5:32162/1‑149 (MQ=255)                      cgcATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACcgc                                                                                                                              <  6:32162/149‑1 (MQ=255)                              gTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAgg                                                                                                                      <  5:123845/149‑1 (MQ=255)                                  cTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTgg                                                                                                                  >  8:244228/1‑149 (MQ=255)                                     gCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGc                                                                                                                 <  7:53948/147‑1 (MQ=255)                                       cAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCGCCGTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGa                                                                                                             >  5:235309/1‑149 (MQ=255)                                       cAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGa                                                                                                             >  4:46634/1‑149 (MQ=255)                                       cAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGa                                                                                                             >  4:161007/1‑149 (MQ=255)                                        aaGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGAt                                                                                                            >  1:273497/1‑149 (MQ=255)                                         aGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATc                                                                                                           >  2:52060/1‑149 (MQ=255)                                         aGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATc                                                                                                           <  6:135193/149‑1 (MQ=255)                                          gAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATcc                                                                                                          <  5:3937/149‑1 (MQ=255)                                              gTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGt                                                                                                      <  3:208955/149‑1 (MQ=255)                                               tCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCAGGCACTGCCGGGGctgc                                                                                                                                   >  6:194800/1‑116 (MQ=255)                                               tCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCAGGCACTGCCGGGGctgc                                                                                                                                   <  5:194800/119‑4 (MQ=255)                                                      gCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGt                                                                                                <  1:152499/147‑1 (MQ=255)                                                        ccGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAg                                                                                            >  4:217656/1‑149 (MQ=255)                                                         cGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTgg                                                                                                                                        >  3:217457/1‑104 (MQ=255)                                                         cGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTgg                                                                                                                                        <  4:217457/104‑1 (MQ=255)                                                         cGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAgg                                                                                           >  3:210086/1‑149 (MQ=255)                                                                 ccAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCaa                                                                                    <  5:308115/148‑1 (MQ=255)                                                                 ccAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCaa                                                                                    <  1:162851/148‑1 (MQ=255)                                                                 ccAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCaa                                                                                    <  7:245945/148‑1 (MQ=255)                                                                 ccAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCaa                                                                                    >  8:245945/1‑148 (MQ=255)                                                                 ccAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAAc                                                                                   <  2:273497/149‑1 (MQ=255)                                                                  cAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACa                                                                                  >  6:83094/1‑149 (MQ=255)                                                                    aGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACAtt                                                                                <  3:161007/149‑1 (MQ=255)                                                                     gTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTg                                                                               <  1:124972/149‑1 (MQ=255)                                                                         ggcggcCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCCCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCcga                                                                           <  3:201481/149‑1 (MQ=255)                                                                              ccAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATc                                                                                                           >  4:103975/1‑112 (MQ=255)                                                                              ccAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATc                                                                                                           <  3:103975/112‑1 (MQ=255)                                                                               cAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCaaaa                                                                                                                                  >  6:31452/1‑88 (MQ=255)                                                                               cAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCaaaa                                                                                                                                  <  5:31452/88‑1 (MQ=255)                                                                                aaCGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCa                                                                    >  7:241449/1‑149 (MQ=255)                                                                                 aCGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCAc                                                                   >  5:316386/1‑149 (MQ=255)                                                                                 aCGCCTGGAGCCGGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCAc                                                                   <  6:235309/149‑1 (MQ=255)                                                                                      ttGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCt                                                               >  3:154801/1‑148 (MQ=255)                                                                                      ttGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTg                                                              >  4:178098/1‑149 (MQ=255)                                                                                        gAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGa                                                             <  5:83094/148‑1 (MQ=255)                                                                                             tGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACAATGCCGACGATCCACGCCTGAGCATc                                                        >  1:159921/1‑148 (MQ=255)                                                                                                ccTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAgg                                                    <  3:288037/149‑1 (MQ=255)                                                                                                     aGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCaa                                               <  7:244228/149‑1 (MQ=255)                                                                                                      gACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAg                                              <  3:46634/149‑1 (MQ=255)                                                                                                         gCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTAGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGcc                                           >  1:88247/1‑149 (MQ=255)                                                                                                          ccccTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCt                                          <  2:88247/149‑1 (MQ=255)                                                                                                             cTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGgc                                       <  1:52060/149‑1 (MQ=255)                                                                                                                 ccGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGgcgtg                                    <  2:159921/148‑1 (MQ=255)                                                                                                                 ccGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGCCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGgcgtg                                    <  6:48770/148‑1 (MQ=255)                                                                                                                   gAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGgcgtgcgt                                 >  8:168228/1‑149 (MQ=255)                                                                                                                           ccAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAgcgcgc                         >  2:253397/1‑149 (MQ=255)                                                                                                                           ccAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAgcgcgc                         <  2:126204/149‑1 (MQ=255)                                                                                                                            cAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGcgcgcg                        >  8:22845/1‑149 (MQ=255)                                                                                                                            cAGGCGGGCGCCAGCCGGCAGCCCAGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGcgcgcg                        >  6:77872/1‑149 (MQ=255)                                                                                                                                cgggcgCCAGCCGGCAGCCCAGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCAc                    <  5:77872/149‑1 (MQ=255)                                                                                                                                 gggcgCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCAc                    >  4:96738/1‑148 (MQ=255)                                                                                                                                   gcgcCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAg                 <  7:168228/149‑1 (MQ=255)                                                                                                                                               aGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGGCTGCTTGCAg     <  1:253397/149‑1 (MQ=255)                                                                                                                                                gCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGGCTGCTTGCAgg    <  4:210086/149‑1 (MQ=255)                                                                                                                                                 ccccGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGGCTGCTTGCAggg   <  3:217656/149‑1 (MQ=255)                                                                                                                                                  cccGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGgctgc           >  8:123879/1‑140 (MQ=255)                                                                                                                                                  cccGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGgctgc           <  7:123879/140‑1 (MQ=255)                                                                                                                                                  cccGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCATGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGGCTGCTTGCAGGGc  >  6:112159/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAACGCCTGGAAATAACGGCGCGCATCCGGTACCTGGCCAAGAAAGTCGACTTGCCCGGCAATGCCAAGTTCGGCGGCCAACGCCTTGAGCCTGGCCTCCAGACGCCCCTTGCCGAGGATTGCCAGGCGGGCGCCAGCCGGCAGCCCCGGCAGTGCCTGGGCAAAACCGCGCAACAGGGTGGCCTGATCCTTGTCCGGGTGCAGGCGCCCAACATTGCCGACGATCCACGCCTGAGCATCCAGGCCCAAGGCCTGGCGTGCGTCAGCGCGCGGCACCAGGGCTGCTTGCAGGGC  >  NC_002947/5624948‑5625241

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