Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,710,169	G→T	5.9%	A104D (GCC→GAC)	PP_4168 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,710,169	0	G	T	5.9%	318.8 / 6.0	98	A104D (GCC→GAC)	PP_4168	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (48/44);  new base T (4/2);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTAAGCCGCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCA  >  NC_002947/4710023‑4710316                                                                                                                                                      |                                                                                                                                                     aGGTAAGCCGCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGcc                                                                                                                                                    >  3:317188/1‑149 (MQ=255)    tgaGCCGCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGt                                                                                                                                                 <  5:210946/147‑1 (MQ=255)      aGCCGCACTGTGCAGCACCTGCTCGTGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTg                                                                                                                                                <  2:65897/148‑1 (MQ=255)       gCCGCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGGTGATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGc                                                                                                                                               >  1:300105/1‑149 (MQ=255)              tgtgCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGt                                                                                                                                       >  5:320467/1‑149 (MQ=255)                tgCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCATCACCCACAGCTTCAAATCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGTCCGGTGCCACgtcattc                                                                                                                                     >  4:289811/1‑143 (MQ=255)                tgCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTgc                                                                                                                                     <  6:400320/149‑1 (MQ=255)                  cagcaCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTgcgc                                                                                                                                   >  7:325120/1‑149 (MQ=255)                   agcaCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGTCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCgct                                                                                                                                  >  2:15402/1‑149 (MQ=255)                     caCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCgctgc                                                                                                                                >  4:27962/1‑149 (MQ=255)                     caCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCgctgc                                                                                                                                <  2:311512/149‑1 (MQ=255)                     caCCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCgctgc                                                                                                                                >  1:78063/1‑149 (MQ=255)                        cTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCtt                                                                                                                              >  3:205136/1‑148 (MQ=255)                            tcgCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCt                                                                                                                         >  6:206591‑M2/4‑149 (MQ=255)                                gACCATGGGCAGTCAACCGGCACTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTTCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCg                                                                                                                     <  1:303940/149‑1 (MQ=255)                                   cATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGgcg                                                                                                                  >  6:395579/1‑149 (MQ=255)                                    aTGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATGCCGCGGCCGGTGCCACGGCGTGCGCTGCTTCAcc                                                                                                                          <  1:52641/141‑1 (MQ=255)                                    aTGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCAcc                                                                                                                          >  2:52641/1‑140 (MQ=255)                                        gCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCAcgc                                                                                                             >  1:328042/1‑149 (MQ=255)                                          aGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCgg                                                                                                           <  8:325120/149‑1 (MQ=255)                                          aGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCgg                                                                                                           <  2:328042/149‑1 (MQ=255)                                           gTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGCCG‑CGCATCCGGTGcgt                                                                                                                                             <  8:32237/114‑3 (MQ=255)                                           gTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGCCG‑CGCATCCGGTGcgt                                                                                                                                             >  7:26266/1‑112 (MQ=255)                                           gTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGCCG‑CGCATCCGGTGcgt                                                                                                                                             >  7:32237/1‑112 (MQ=255)                                           gTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGCCG‑CGCATCCGGTGcgt                                                                                                                                             <  8:26266/114‑3 (MQ=255)                                           gTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCggg                                                                                                          <  7:357021/149‑1 (MQ=255)                                              aaCCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGggcc                                                                                                       >  4:200831/1‑149 (MQ=255)                                              aaCCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGggcc                                                                                                       >  8:1455/1‑149 (MQ=255)                                               aCCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGggccg                                                                                                      <  3:196228/149‑1 (MQ=255)                                                ccGGCGCTCGG‑CATCCGTCTTGGCCCGGCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGggccg                                                                                                      >  1:217431/1‑149 (MQ=255)                                                 cGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGTCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCCCAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGggccggc                                                                                                    <  3:23307/149‑1 (MQ=255)                                                 cGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGggccggc                                                                                                    <  6:103215/149‑1 (MQ=255)                                                      cTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTc                                                                                               >  7:95637/1‑149 (MQ=255)                                                         gg‑CATCCGTCTTGGTCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGc                                                                                            <  1:15402/149‑1 (MQ=255)                                                         gg‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGc                                                                                            <  8:95637/149‑1 (MQ=255)                                                         gg‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCACTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGc                                                                                            >  2:159596/1‑149 (MQ=255)                                                             aTCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGaa                                                                                         >  2:212492/1‑149 (MQ=255)                                                              tCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAAc                                                                                        <  2:78063/149‑1 (MQ=255)                                                                cGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACgg                                                                                      <  1:95845/149‑1 (MQ=255)                                                                cGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACg                                                                                       <  4:353431/148‑1 (MQ=255)                                                                      ggCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACgg                                                                                >  3:289616/1‑149 (MQ=255)                                                                       gCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACgg                                                                                <  2:300105/148‑1 (MQ=255)                                                                        ccc‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCt                                                                              >  6:251128/1‑149 (MQ=255)                                                                         cc‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCtt                                                                             <  2:91349/149‑1 (MQ=255)                                                                         cc‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCtt                                                                             <  4:317188/149‑1 (MQ=255)                                                                             cgGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTg                                                                          <  2:40567/149‑1 (MQ=255)                                                                              ggCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTgg                                                                         <  5:395579/149‑1 (MQ=255)                                                                                cTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTc                                                                       >  1:159721/1‑149 (MQ=255)                                                                                 tCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGgc                                                                                                                   >  6:168362/1‑104 (MQ=255)                                                                                 tCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGgc                                                                                                                   <  5:168362/104‑1 (MQ=255)                                                                                 tCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCtcgtcg                                                                      >  4:222230/1‑149 (MQ=255)                                                                                 tCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCg                                                                      >  2:63721/1‑149 (MQ=255)                                                                                   aGCTTCACCGTGGCGCAGCACCCACAGCTTCCCAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCt                                                                    <  8:190825/149‑1 (MQ=255)                                                                                   aGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGc                                                                     <  7:1455/148‑1 (MQ=255)                                                                                     cTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCtgt                                                                  >  7:72864/1‑149 (MQ=255)                                                                                        cACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTgg                                                                <  4:205136/148‑1 (MQ=255)                                                                                        cACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTgg                                                                >  2:380012/1‑148 (MQ=255)                                                                                           cGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAAc                                                            <  3:200831/149‑1 (MQ=255)                                                                                            gTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACt                                                           <  3:27962/149‑1 (MQ=255)                                                                                             tGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCTTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTg                                                          <  3:289811/149‑1 (MQ=255)                                                                                             tGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTg                                                          >  2:357026/1‑149 (MQ=255)                                                                                                cgcAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCg                                                       <  2:217431/149‑1 (MQ=255)                                                                                                 gcagcaCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCg                                                       >  6:19736/1‑148 (MQ=255)                                                                                                   agcaCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCACTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAg                                                    <  1:159596/149‑1 (MQ=255)                                                                                                     caCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAAc                                                            <  5:40069/139‑1 (MQ=255)                                                                                                     caCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAAc                                                            >  6:40069/1‑139 (MQ=255)                                                                                                      aCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCgg                                                 <  8:72864/149‑1 (MQ=255)                                                                                                      aCCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGACACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAACCgg                                                 >  4:41594/1‑149 (MQ=255)                                                                                                        ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑cc                                              >  5:80812/1‑149 (MQ=255)                                                                                                        ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑cc                                              >  7:28482/1‑149 (MQ=255)                                                                                                        ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑cc                                              >  8:239227/1‑149 (MQ=255)                                                                                                        ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑cc                                              >  4:145807/1‑149 (MQ=255)                                                                                                        ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑cc                                              >  3:333244/1‑149 (MQ=255)                                                                                                        ccACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑cc                                              >  2:100090/1‑149 (MQ=255)                                                                                                          acagcttcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGa                                            >  6:82674/1‑149 (MQ=255)                                                                                                          acagcttcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGa                                            >  1:147982/1‑149 (MQ=255)                                                                                                          acagcttcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGa                                            >  2:298370/1‑149 (MQ=255)                                                                                                          acagcttcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGa                                            <  3:17713/149‑1 (MQ=255)                                                                                                           cagcttcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGAt                                           <  5:251128/149‑1 (MQ=255)                                                                                                                tcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCtgtg                                      >  2:75177/1‑149 (MQ=255)                                                                                                                tcacagcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCtgtg                                      >  3:146136/1‑149 (MQ=255)                                                                                                                    agcttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAg                                  >  6:381637/1‑149 (MQ=255)                                                                                                                      cttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAgg                                 >  8:8686/1‑148 (MQ=255)                                                                                                                      cttGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAgg                                 >  5:42932/1‑148 (MQ=255)                                                                                                                          ggTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACt                             <  5:206591/148‑1 (MQ=255)                                                                                                                          ggTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTg                            <  7:297522/149‑1 (MQ=255)                                                                                                                          ggTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTg                            >  2:262770/1‑149 (MQ=255)                                                                                                                               cTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCa                       <  2:159721/149‑1 (MQ=255)                                                                                                                                tCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCt                         >  1:4617/1‑146 (MQ=255)                                                                                                                                  gTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGa                     <  1:357026/148‑1 (MQ=255)                                                                                                                                      ggCACCGGCCG‑CGCATCCGGTGCCACGGCGTGCGCTGCGTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGCCTGTc                >  7:26233/2‑149 (MQ=255)                                                                                                                                      cGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTc                <  1:100090/149‑1 (MQ=255)                                                                                                                                         aCCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGcc             >  4:156944/1‑149 (MQ=255)                                                                                                                                         aCCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGcc             >  8:211495/1‑149 (MQ=255)                                                                                                                                            ggATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGaa          <  6:320467/149‑1 (MQ=255)                                                                                                                                             gATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAg         <  3:145807/149‑1 (MQ=255)                                                                                                                                             gATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAg         <  5:381637/149‑1 (MQ=255)                                                                                                                                                    cgGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCATGCTGCGAACGGCCACGGCTTCTGGGCGCTGTGGAAACTGCCGAAGCGGCCCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTtcat    <  1:75177/149‑1 (MQ=255)                                                                                                                                                     ggCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTtcatca  <  1:298370/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGGTAAGCCGCACTGTGCAGCACCTGCTCGCGACCATGGGCAGTCAACCGGCGCTCGG‑CATCCGTCTTGGCCC‑GCGGCTCAGCTTCACCGTGGCGCAGCACCCACAGCTTCACAGCTTGGGTTCCTCGTCACGCACCGGATG‑CGCGGCCGGTGCCACGGCGTGCGCTGCTTCACCTTCCGGCGCACGCGGGGCCGGCCAGTCTGCGAACGGCCACGGCTTCTGGTCGCTGTGGAAACTGCCGAAGCGG‑CCGATCTGTGCCAGGTACTGGCTCAGACTGTCGCCGAAGTTCATCA  >  NC_002947/4710023‑4710316

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