Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,054,368	G→A	6.0%	N86N (AAC→AAT)	PP_1829 ←	alpha/beta family hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,054,368	0	G	A	6.0%	201.7 / 2.2	66	N86N (AAC→AAT)	PP_1829	alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (33/33)
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.21e-01

CGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAGGCTGCCGTGC  >  NC_002947/2054223‑2054515                                                                                                                                                  |                                                                                                                                                    cGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTg                                                                                                                                                  <  1:196932/149‑1 (MQ=255) cGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTg                                                                                                                                                  <  2:103431/149‑1 (MQ=255) cGAAGCGCGGGTCGGCGTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTg                                                                                                                                                  <  8:15005/149‑1 (MQ=255)      cgcgGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCcag                                                                                                                                             <  3:322904/149‑1 (MQ=255)            tCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTTGTAACGCACGCTGGCAATGTGTTCGTTGGCcagcagcag                                                                                                                                       >  7:210292/1‑149 (MQ=255)             cGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAgg                                                                                                                                      >  1:175195/1‑149 (MQ=255)                    gAGCTTGCGGCTCCAGGCGACCACATCGGCGACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGAc                                                                                                                               >  1:30281/1‑149 (MQ=255)                    gAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGa                                                                                                                                <  8:51139/148‑1 (MQ=255)                     aGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACg                                                                                                                              <  3:131794/149‑1 (MQ=255)                      gCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCTACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGc                                                                                                                             <  2:288956/149‑1 (MQ=255)                      gCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGc                                                                                                                             <  2:30281/149‑1 (MQ=255)                        ttGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCtt                                                                                                                           >  2:300559/1‑149 (MQ=255)                          gCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCtt                                                                                                                           <  1:300559/147‑1 (MQ=255)                           cGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCa                                                                                                                         <  5:294552/148‑1 (MQ=255)                                 cAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTg                                                                                                                    <  1:7468/147‑1 (MQ=255)                                 cAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTg                                                                                                                    <  1:144716/147‑1 (MQ=255)                                 cAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTc                                                                                                                  <  6:334004/149‑1 (MQ=255)                                  aGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTc                                                                                                                  >  3:102672/1‑148 (MQ=255)                                  aGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCg                                                                                                                 >  7:320791/1‑149 (MQ=255)                                    gCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGAc                                                                                                               >  4:12078/1‑149 (MQ=255)                                         cacaTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGc                                                                                                          >  4:332913/1‑149 (MQ=255)                                          acaTCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGcc                                                                                                         <  3:332913/149‑1 (MQ=255)                                             tCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTg                                                                                                                                                  <  8:78842/105‑1 (MQ=255)                                             tCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTg                                                                                                                                                  >  7:78842/1‑105 (MQ=255)                                              cGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGAc                                                                                                     >  3:34734/1‑149 (MQ=255)                                                  cacaTAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCAGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCg                                                                                                 >  7:172931/1‑149 (MQ=255)                                                           cTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCa                                                                                                                         <  4:12646/116‑1 (MQ=255)                                                           cTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCa                                                                                                                         >  3:12646/1‑116 (MQ=255)                                                            tCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTAcc                                                                                       >  6:154048/1‑149 (MQ=255)                                                                       gTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGcggtcgg                                                                            >  3:212330/1‑149 (MQ=255)                                                                          aCGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTgg                                                                         >  1:35132/1‑149 (MQ=255)                                                                              tCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGcc                                                                     >  7:282347/1‑149 (MQ=255)                                                                                gTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCg                                                                    >  6:156776/1‑148 (MQ=255)                                                                                      cGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCgg                                                             <  4:34734/149‑1 (MQ=255)                                                                                         ggCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCga                                                          <  5:154048/149‑1 (MQ=255)                                                                                             ggctggctGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCgatgat                                                      <  2:175195/149‑1 (MQ=255)                                                                                              gctggctgGCGGCGACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCATCAGCATGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATc                                                     <  6:232141/149‑1 (MQ=255)                                                                                               ctggctggCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCACGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCa                                                    <  4:102672/149‑1 (MQ=255)                                                                                                 ggctggcGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGc                                                  <  3:12078/149‑1 (MQ=255)                                                                                                  gctggcGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAg                                                   <  6:286895/147‑1 (MQ=255)                                                                                                    tggcGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGcacc                                               >  1:302387/1‑149 (MQ=255)                                                                                                    tggcGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGcac                                                >  3:193573/1‑148 (MQ=255)                                                                                                          ccACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCATTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTggc                                         >  8:313616/1‑149 (MQ=255)                                                                                                           cACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTggcg                                        >  2:20183/1‑149 (MQ=255)                                                                                                            aCCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGACGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTggcgg                                       >  6:137437/1‑149 (MQ=255)                                                                                                            aCCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCATTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTggcgg                                       >  5:41242/1‑149 (MQ=255)                                                                                                              cccGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGa                                                                  >  1:37397/1‑120 (MQ=255)                                                                                                              cccGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGa                                                                  <  2:37397/120‑1 (MQ=255)                                                                                                              cccGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGgcggcg                                     >  3:297675/1‑149 (MQ=255)                                                                                                               ccGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGcggcgg                                    <  4:297675/149‑1 (MQ=255)                                                                                                                cgcgTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCggg                                   <  8:320791/149‑1 (MQ=255)                                                                                                                        tCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTg                           >  7:172007/1‑149 (MQ=255)                                                                                                                        tCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTg                           >  1:319286/1‑149 (MQ=255)                                                                                                                        tCGTAACGCACGCTGGCAATGTGTTCATTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTg                           <  6:41242/149‑1 (MQ=255)                                                                                                                            aaCGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGc                       >  7:8600/1‑149 (MQ=255)                                                                                                                              cgcacgcTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCgg                     <  2:35132/149‑1 (MQ=255)                                                                                                                               gcacgcTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCgg                     >  1:218870/1‑148 (MQ=255)                                                                                                                               gcacgcTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCAGGCCCGACGCCGGGTTGTTAACGTCGCGGTCGGTGGGGCCGGAGCCGGCCATGATCAGCACCACTGGCGGCGGGGTGGCCTgcagatgc                    >  1:336281/1‑142 (MQ=255)                                                                                                                                       gCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGcagcagcag            <  8:210292/149‑1 (MQ=255)                                                                                                                                       gCAATGTGTTCATTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTAGCCTGCTGCTGccgcaacag            <  7:313616/149‑1 (MQ=255)                                                                                                                                        cAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGCGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAgg           >  2:247523/1‑149 (MQ=255)                                                                                                                                        cAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAgg           <  4:212330/149‑1 (MQ=255)                                                                                                                                          aTGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAGGCt         <  4:193573/149‑1 (MQ=255)                                                                                                                                              gtTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTCGCCTGTTGCGGCAGCAGCAGGCTGCCg     >  8:33878/1‑149 (MQ=255)                                                                                                                                                tCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAGGCTGCCGTg   >  8:81801/1‑149 (MQ=255)                                                                                                                                                 cGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAGGCTGCCGTg   <  8:282347/148‑1 (MQ=255)                                                                                                                                                  gTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAGGCTGCCGTgc  <  7:81801/148‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CGAAGCGCGGGTCGGCTTTGAGCTTGCGGCTCCAGGCGACCACATCGGCCACATAGCGCTCCACGCTGAGGTCACGCTCGTCAGGCGTGGCCGGCTGGCTGGCGGCCACCCCGCGTTTGTCGTAACGCACGCTGGCAATGTGTTCGTTGGCCAGCAGCAGGGCCAGACGCTTCAGGTTGTCGACCCGGCCCGACGCCGGGTTGTTACCGTCGCGGTCGGTGGGGCCGGAGCCGGCGATGATCAGCACCACTGGCGGCGGGGTGGCCTGTTGCGGCAGCAGCAGGCTGCCGTGC  >  NC_002947/2054223‑2054515

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