Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,299,365	G→A	7.9%	F271F (TTC→TTT)	PP_1134 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,299,365	0	G	A	7.9%	220.2 / 6.5	76	F271F (TTC→TTT)	PP_1134	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/38);  new base A (3/3);  total (35/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.29e-01

AGCTTGTAGTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGCCCCCAGGG  >  NC_002947/1299219‑1299513                                                                                                                                                   |                                                                                                                                                       aGCTTGTAGTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCgg                                                                                                                                                        >  3:225977/1‑147 (MQ=255) aGCTTGTAGTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGaa                                                                                                                                                      <  2:181620/149‑1 (MQ=255)  gCTTGTAGTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑g                                                                                                                                                    <  7:199327/149‑1 (MQ=255)      gTAGTCGGCGGCACCAAGTTTTGGGTGGCCGGCATACTTGACCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCTATACCGAAGCCCACCGATAACGTTTTGGCTTCCATGTAGACACCTTAGGCGTGCAGCATGCGGAA‑Ggt                                                                                                                                                  <  2:48504/147‑1 (MQ=255)         gTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTc‑gc                                                                                                                                            <  3:321693/149‑1 (MQ=255)         gTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTc‑gc                                                                                                                                            >  8:48924/1‑149 (MQ=255)          tCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTc‑gcg                                                                                                                                           >  1:313630/1‑149 (MQ=255)          tCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTc‑gcg                                                                                                                                           >  7:121214/1‑149 (MQ=255)            ggcggcACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCgg                                                                                                                                          <  2:13627/148‑1 (MQ=255)            ggcggcACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACCGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTCCGCATGCTGCACGCCAAA‑GGTGTctaca                                                                                                                                           <  3:68554/148‑5 (MQ=38)                  aCCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTCCGCATGCTGCACGCCAAA‑GGTGTctac                                                                                                                                            <  1:346122/141‑4 (MQ=38)                  aCCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTCCGCATGCTGCACGCCAAA‑GGTGTctac                                                                                                                                            >  2:346122/1‑138 (MQ=38)                    ccAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTgg                                                                                                                                 >  5:84845/1‑149 (MQ=255)                     cAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGc                                                                                                                                >  3:310146/1‑149 (MQ=255)                     cAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTCCGCATGCTGCACGCCAAA‑GGTGTc                                                                                                                                               <  5:294002/135‑1 (MQ=38)                     cAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTCCGCATGCTGCACGCCAAA‑GGTGTc                                                                                                                                               >  6:294002/1‑135 (MQ=38)                      aGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGcc                                                                                                                               >  7:142482/1‑149 (MQ=255)                           tGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAgg                                                                                                                           <  5:4779/148‑1 (MQ=255)                                 gCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATa                                                                                                                    >  1:307184/1‑149 (MQ=255)                                 gCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATa                                                                                                                    >  2:113172/1‑149 (MQ=255)                                     gCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACca                                                                                                                 >  1:353488/1‑148 (MQ=255)                                         aCTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATgg                                                                                                            >  4:56412/1‑149 (MQ=255)                                          cTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATgg                                                                                                            <  3:56412/148‑1 (MQ=255)                                          cTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGgc                                                                                                           <  6:84845/149‑1 (MQ=255)                                           ttGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGgcg                                                                                                          >  4:36482/1‑149 (MQ=255)                                              cccAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGa                                                                                                       <  5:252303/149‑1 (MQ=255)                                                cAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGt                                                                                                     <  3:36482/149‑1 (MQ=255)                                                cAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGt                                                                                                     <  8:121214/149‑1 (MQ=255)                                                 aGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTg                                                                                                    <  8:142482/149‑1 (MQ=255)                                                 aGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTg                                                                                                    <  6:74212/149‑1 (MQ=255)                                                   gCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAACGGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGc                                                                                                   <  1:91368/149‑1 (MQ=255)                                                   gCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGc                                                                                                   >  4:364245/1‑148 (MQ=255)                                                    cGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCaa                                                                                                 <  3:364245/149‑1 (MQ=255)                                                    cGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCaa                                                                                                 <  7:48924/149‑1 (MQ=255)                                                     gTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAAc                                                                                                >  7:68526/1‑149 (MQ=255)                                                            gTCGATCAGCTGTTGCGGGTGCTCGATTCGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTcg                                                                                          <  7:151298/148‑1 (MQ=255)                                                               gATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACGTGGCGCGAGTGCAACTGCTcgccgc                                                                                      <  4:310146/149‑1 (MQ=255)                                                                 tCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCtt                                                                                    <  1:266727/149‑1 (MQ=255)                                                                      gTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGGTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAg                                                                               >  8:244166/1‑149 (MQ=255)                                                                      gTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAg                                                                               <  1:113172/149‑1 (MQ=255)                                                                           cGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATg                                                                                                             >  2:340829/1‑114 (MQ=255)                                                                           cGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATg                                                                                                             <  1:340829/114‑1 (MQ=255)                                                                                   cGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGt                                                                   >  7:338951/1‑148 (MQ=255)                                                                                    gATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGt                                                                   >  4:327178/1‑147 (MQ=255)                                                                                        cGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTg                                                             <  4:225977/149‑1 (MQ=255)                                                                                         ggAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTgg                                                            >  5:168450/1‑149 (MQ=255)                                                                                          gAAGACCACCGAGAACGGTTTGGCTTCCATGTAGACAACTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTTGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCTGTCAGTTGGc                                                           <  6:303559/149‑1 (MQ=255)                                                                                               ccACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCt                                                      <  8:68526/149‑1 (MQ=255)                                                                                                cACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTCTGCGGGCGCTGTGGCCCTGGGCCacagc                                                                                                                   <  1:214197/88‑5 (MQ=255)                                                                                                 aCCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCa                                                                 >  8:151134/1‑136 (MQ=255)                                                                                                 aCCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCa                                                                 <  7:151134/136‑1 (MQ=255)                                                                                                 aCCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTc                                                    >  5:111512/1‑149 (MQ=255)                                                                                                   cGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAgg                                                                                                                           >  8:101446/1‑76 (MQ=255)                                                                                                   cGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGa                                                  <  6:168450/149‑1 (MQ=255)                                                                                                   cGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGa                                                  <  3:327178/149‑1 (MQ=255)                                                                                                     agaACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGAt                                                 >  5:583/1‑148 (MQ=255)                                                                                                     agaACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATc                                                <  2:208225/149‑1 (MQ=255)                                                                                                      gaACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATc                                                >  8:339125/1‑148 (MQ=255)                                                                                                       aaCGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCaa                                              <  2:307184/149‑1 (MQ=255)                                                                                                          ggTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAg                                           >  6:205686/1‑149 (MQ=255)                                                                                                           gTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAgg                                          <  3:13503/149‑1 (MQ=255)                                                                                                                gCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGt                                     <  2:340616/149‑1 (MQ=255)                                                                                                                      aTGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCt                                >  3:32427/1‑148 (MQ=255)                                                                                                                           gACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGc                                                                                                                         <  6:249246/54‑1 (MQ=255)                                                                                                                           gACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGc                                                                                                                         >  5:249246/1‑54 (MQ=255)                                                                                                                           gACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCt                          >  3:33119/1‑149 (MQ=255)                                                                                                                           gACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCt                          <  6:583/149‑1 (MQ=255)                                                                                                                           gACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCt                          <  8:338951/149‑1 (MQ=255)                                                                                                                           gACACCTTCCGCATGCTGCACGCCAAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCt                          >  6:108612/1‑149 (MQ=255)                                                                                                                            acacCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTg                         >  2:328964/1‑149 (MQ=255)                                                                                                                              acCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTgg                        >  6:146460/1‑148 (MQ=255)                                                                                                                              acCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGa                       <  4:223345/149‑1 (MQ=255)                                                                                                                                   tGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAgc                    <  2:353488/147‑1 (MQ=255)                                                                                                                                            gcaTGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGcc         <  4:109366/149‑1 (MQ=255)                                                                                                                                               tGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGACGCTTTCAAGCACTACACCGGTCAGTTTGCCATCTTCGATAAACAGGTCGGTGACCTTCTCCTGGAAGCTCACTTCGTcccc       >  5:280972/1‑148 (MQ=255)                                                                                                                                                gCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGcccc       >  3:201957/1‑147 (MQ=255)                                                                                                                                                gCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGcccc       >  2:7763/1‑147 (MQ=255)                                                                                                                                                   gAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGCCCCCAggg  <  2:313630/149‑1 (MQ=255)                                                                                                                                                   gAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGCCCCCAggg  >  4:372723/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                       AGCTTGTAGTCGGCGGCACCCAGTTTTGGGTGGCCGGCATACTTGCCCAGGCGTGCCTTGTCGATCAGCGTTTGCGGGTGCTCGATACGGAAGCCCACCGAGAACGGTTTGGCTTCCATGTAGACACCTTTGGCGTGCAGCATGCGGAA‑GGTGTC‑GCGGGCGCTGTGGCCCAGGGCCAATACCACATGGCGCGAGTGCAACTGCTCGCCGCTTTCAAGCACTACACCGGTCAGTTGGCCATCTTCGATCAACAGGTCGGTGACCTTCTCCTGGAAGCGCACTTCGGCCCCCAGGG  >  NC_002947/1299219‑1299513

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