Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,904,888	T→A	6.5%	intergenic (‑55/+180)	PP_4312 ← / ← PP_4313	hypothetical protein/peptidylprolyl isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,904,888	0	T	A	6.5%	233.7 / 4.0	77	intergenic (‑55/+180)	PP_4312/PP_4313	hypothetical protein/peptidylprolyl isomerase
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base A (2/3);  total (38/39)
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.64e-01

CCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAA  >  NC_002947/4904746‑4905035                                                                                                                                                 |                                                                                                                                                         ccTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGc                                                                                                                                                    <  2:120301/149‑1 (MQ=255) ccTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGc                                                                                                                                                    >  5:186951/1‑149 (MQ=255)     cAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAgg                                                                                                                                                 <  3:276023/148‑1 (MQ=255)       ctgcCATCGCCGCGCCGTTGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGgcg                                                                                                                                              >  5:184887/3‑149 (MQ=255)         gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCTGCAACTAGCCGGgcgat                                                                                                                                            >  4:207077/1‑147 (MQ=255)         gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCTGCAACTAGCCGGgc                                                                                                                                               <  2:318733/146‑1 (MQ=255)         gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCTGCAACTAGCCGGgc                                                                                                                                               >  1:318733/1‑146 (MQ=255)         gcCAACGGCGCGGCGATGGCAGTGAGCACCAGGCCCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGcc                                                                                                                                            >  2:68873/1‑149 (MQ=255)          ccAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAAAGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGTGTTACCCGCCACATCGCCCTGCAACTAGCCGGgcgat                                                                                                                                            <  1:61352/149‑3 (MQ=255)           cAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCTGCAACTAGCCGGgcgatgt                                                                                                                                          <  3:207077/149‑5 (MQ=255)           cAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCg                                                                                                                                          >  4:113171/1‑149 (MQ=255)             aCGGCGCGGCGATGGCAGTGAGCACCAGGCCCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTg                                                                                                                                        <  1:68873/149‑1 (MQ=255)               ggcgcggcgATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCa                                                                                                                                      <  3:48609/149‑1 (MQ=255)                      cgATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATg                                                                                                                               <  4:179055/149‑1 (MQ=255)                               tGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTg                                                                                                                      <  6:184887/149‑1 (MQ=255)                                 aGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTgg                                                                                                                     <  6:203862/148‑1 (MQ=255)                                  gCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGc                                                                                                                    >  5:294184/1‑148 (MQ=255)                                           tCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGtt                                                                                                                                                      <  4:21365/105‑1 (MQ=255)                                           tCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGtt                                                                                                                                                      >  3:21365/1‑105 (MQ=255)                                            cagcagCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑Agg                                                                                                       <  2:124345/149‑1 (MQ=255)                                                  caAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCaga                                                                                                 <  6:186951/149‑1 (MQ=255)                                                         cgacgacgGCGAAAGAATGTTTGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGcccc                                                                                          >  4:6808/1‑149 (MQ=255)                                                          gacgacgGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGcccc                                                                                          >  3:94324/1‑148 (MQ=255)                                                          gacgacgGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGcccc                                                                                          >  2:354553/1‑148 (MQ=255)                                                             gacgGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAAc                                                                                   <  3:218979/149‑1 (MQ=255)                                                             gacgGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAAc                                                                                   <  7:269635/149‑1 (MQ=255)                                                                ggCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCAGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTGAGAGTGCCCC‑‑‑GAACACg                                                                                >  8:329365/1‑149 (MQ=255)                                                                  cGAAAGAATGTTTGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGtt                                                                              <  3:6808/149‑1 (MQ=255)                                                                  cGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGtt                                                                              <  6:200314/149‑1 (MQ=255)                                                                   gAAAGAATGTTGGCAACTTCATCAAGGCAGGGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGAAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGTTCAGAGTGCCCC‑‑‑GAACCCGtt                                                                              >  1:215388/1‑149 (MQ=255)                                                                   gAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGttt                                                                             >  2:299767/1‑149 (MQ=255)                                                                     aaGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTgg                                                                           <  1:354553/149‑1 (MQ=255)                                                                      aGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGt                                                                          <  6:294184/149‑1 (MQ=255)                                                                              ggCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCt                                                                  >  3:98668/1‑149 (MQ=255)                                                                              ggCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGCGCCCC‑‑‑GAACCCGTTTGGTtcct                                                                      >  1:279876/1‑145 (MQ=255)                                                                              ggCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGCGCCCC‑‑‑GAACCCGTTTGGTtcct                                                                      <  2:279876/145‑1 (MQ=255)                                                                                 aaCTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTc                                                                                                    <  3:274953/115‑1 (MQ=255)                                                                                 aaCTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTc                                                                                                    >  4:274953/1‑115 (MQ=255)                                                                                  aCTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAg                                                              >  2:285321/1‑149 (MQ=255)                                                                                       atcaAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGaa                                                         >  1:156949/1‑149 (MQ=255)                                                                                        tcaAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAAt                                                        <  2:215388/149‑1 (MQ=255)                                                                                          aaGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGc                                                      <  3:113171/149‑1 (MQ=255)                                                                                           aggca‑ggcTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGcc                                                     >  7:246954/1‑149 (MQ=255)                                                                                           aggca‑ggcTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGcc                                                     <  4:98668/149‑1 (MQ=255)                                                                                            ggca‑ggcTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCg                                                    >  2:81166/1‑149 (MQ=255)                                                                                            ggca‑ggcTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCg                                                    >  5:156876/1‑149 (MQ=255)                                                                                               a‑ggcTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGtta                                                 >  5:302030/1‑148 (MQ=255)                                                                                                     cggggggCAGAGAATATCGGTG‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAg                                            <  4:203576/144‑1 (MQ=255)                                                                                                      gtgtGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAgg                                           >  1:218316/1‑149 (MQ=255)                                                                                                      gtgtGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAg                                            <  6:156876/148‑1 (MQ=255)                                                                                                         tgGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGc                                          >  7:97753/1‑147 (MQ=255)                                                                                                           gCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAg                                       <  6:302030/148‑1 (MQ=255)                                                                                                            cAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGc                                     >  8:221186/1‑149 (MQ=255)                                                                                                             agagAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGt                                                                                                                                    >  6:40376/1‑58 (MQ=255)                                                                                                             agagAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGt                                                                                                                                    <  5:40376/58‑1 (MQ=255)                                                                                                              gagaATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGcc                                    >  6:93793/1‑148 (MQ=255)                                                                                                              gagaATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGcc                                    >  2:231455/1‑148 (MQ=255)                                                                                                                gaATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAg                                 <  8:97753/149‑1 (MQ=255)                                                                                                                 aaTATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAgg                                >  5:338336/1‑149 (MQ=255)                                                                                                                  atatCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgg                                <  1:118651/148‑1 (MQ=255)                                                                                                                   tatCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGtt                              >  2:277399/1‑149 (MQ=255)                                                                                                                   tatCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                               >  1:309361/1‑148 (MQ=255)                                                                                                                   tatCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGt                               >  6:346006/1‑148 (MQ=255)                                                                                                                     tCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTTCAGTTAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTa                             <  8:246954/148‑1 (MQ=255)                                                                                                                     tCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGt                                                                          <  4:221840/103‑1 (MQ=255)                                                                                                                     tCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGt                                                                          >  3:221840/1‑103 (MQ=255)                                                                                                                       gggT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATTCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCCTTGCCTAgggaacg                           <  7:329365/148‑6 (MQ=255)                                                                                                                        ggT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAgggg                         <  2:346614/149‑1 (MQ=255)                                                                                                                         gT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGc                        >  7:113355/1‑149 (MQ=255)                                                                                                                              aCCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTgg                     <  8:113355/148‑1 (MQ=255)                                                                                                                                 cGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGgtgt                 >  6:369211/1‑149 (MQ=255)                                                                                                                                     acaTCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTgcg             >  6:321394/1‑149 (MQ=255)                                                                                                                                       aTCGCCCGGCTAGTTTCAGGGCGCCCGTGCAGGGCATGCCATCCTTGGCAAAGCGTGTCTGGGTCAGAGTGCCCCCCAGAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGtt                <  1:299767/149‑1 (MQ=255)                                                                                                                                          gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTg               >  7:101157/1‑142 (MQ=255)                                                                                                                                          gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTg               <  8:101157/142‑1 (MQ=255)                                                                                                                                          gCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAAtt        <  4:94324/149‑1 (MQ=255)                                                                                                                                                cTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGATTGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTaa  <  1:285321/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                         CCTCCAGCGCCAACGGCGCGGCGATGGCAGTGAGCACCAGGCTCAGCAGCAAGTGACGACGACGGCGAAAGAATGTTGGCAACTTCATCAAGGCA‑GGCTCGTGTGGCAGAGAATATCGGGT‑TTACCCGCCACATCGCCCGGCTAGTTGCAGGGCGCCCGTGCAGTGAATGCCATCCTTGGCAACGCGTG‑‑AGGGTCAGAGTGCCCC‑‑‑GAACCCGTTTGGTTCCTCCCTGGAGGTGAATGCCGTTCCCTAGGCAAGGCCTAGGTTAGGGGCTGGGTGTTGCGCAATTGTGTAA  >  NC_002947/4904746‑4905035

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