Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,032,761	C→G	7.6%	R94G (CGA→GGA)	PP_4435 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,032,761	0	C	G	7.6%	249.1 / 8.6	79	R94G (CGA→GGA)	PP_4435	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base G (4/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCTAAA  >  NC_002947/5032615‑5032909                                                                                                                                                    |                                                                                                                                                     cGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTgcgc                                                                                                                                                      >  8:369680/1‑147 (MQ=255) cGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGa                                                                                                                                                    <  3:354724/149‑1 (MQ=255)  gAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAg                                                                                                                                                   >  3:231884/1‑149 (MQ=255)  gAACATTGACAGGCTTCTTGAGGAGGACGGTAGTCACAAGTAGAAAATCCTCACAGTGACTATCTATCTGAACGACATCAACACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCTG‑AAGGTCAGTCTCCTGCGCGAg                                                                                                                                                   >  1:417113/1‑149 (MQ=255)   aaCATTGACAGGCTTCTTGAGTAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGc                                                                                                                                                  <  6:141522/149‑1 (MQ=255)   aaCATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGc                                                                                                                                                  <  8:359498/149‑1 (MQ=255)       ttGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCGCGCAGGAGCCtga                                                                                                                                              <  8:358780/149‑2 (MQ=255)         gACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTgtgt                                                                                                                                            >  5:28624/1‑149 (MQ=255)           cAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTg                                                                                                                                          <  1:79101/149‑1 (MQ=255)            atggttcttGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGa                                                                                                                                         <  5:4582/145‑1 (MQ=255)            aGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCGCGCAGGAGCCtgac                                                                                                                                             >  2:248370/1‑143 (MQ=255)            aGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCGCGCAGGAGCCtgac                                                                                                                                             <  1:248370/145‑3 (MQ=255)            aGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCtgt                                                                                                                                              <  5:7242/144‑1 (MQ=255)            aGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCtgt                                                                                                                                              >  6:7242/1‑144 (MQ=255)              gCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAgg                                                                                                                                       >  8:261772/1‑149 (MQ=255)                   ttGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAggcggc                                                                                                                                   >  8:113887/1‑148 (MQ=255)                           gCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaa                                                                                                                          >  2:73880/1‑149 (MQ=255)                           gCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTc                                                                                                                            >  3:204425/1‑147 (MQ=255)                            ccGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaag                                                                                                                         <  2:186105/149‑1 (MQ=255)                                gaGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaagccaa                                                                                                                     <  4:324677/148‑1 (MQ=255)                                  gTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCTACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaagccaagc                                                                                                                   <  1:196323/149‑1 (MQ=255)                                   tCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCaagccaagcc                                                                                                                  <  3:15169/149‑1 (MQ=255)                                      caAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGc                                                                                                               >  5:285472/1‑149 (MQ=255)                                      caAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGc                                                                                                               >  4:151135/1‑149 (MQ=255)                                       aactCGAAAATCCTCACACGGACTATCTATCTGAACGACATCACCACCGTCTTCGGGGCCATGAACGAGGTTTGGGAACCGTGGGGCGCGGTAGGGTCAGGCTCCTGCGCGAGCCTGTGTTg                                                                                                                                          <  1:3918/119‑1 (MQ=255)                                       aaGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGAGTTCGCTGCCATGAACGAGGTTTGGGATGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCt                                                                                                              >  6:72397/1‑149 (MQ=255)                                       aaGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTg                                                                                                                                          >  2:3918/1‑121 (MQ=255)                                       aaGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaa                                                                                                                          >  3:247259/1‑137 (MQ=255)                                       aaGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTcaa                                                                                                                          <  4:247259/137‑1 (MQ=255)                                        aGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTc                                                                                                             <  6:196816/149‑1 (MQ=255)                                            gAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                                          >  4:237704/1‑148 (MQ=255)                                            gAAAATCCTCACAGTGACTATCGATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                                          >  6:70928/1‑148 (MQ=255)                                             aaaaTCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTga                                                                                                        <  7:261772/149‑1 (MQ=255)                                              aaaTCCTCACATTTAATATCTATCTGAACGACATCACCACCGACTTCGCTTCCATGAACGAGGTTTGGGAAGAGTTGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTgag                                                                                                       <  5:263982/149‑1 (MQ=255)                                              aaaTCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTgag                                                                                                       <  1:56517/149‑1 (MQ=255)                                               aaTCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTgaga                                                                                                      <  1:73880/149‑1 (MQ=255)                                                    tCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTAt                                                                                                                                <  8:210123/118‑1 (MQ=255)                                                    tCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTAt                                                                                                                                >  7:210123/1‑118 (MQ=255)                                                      acaGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTcc                                                                                                 <  6:285472/147‑1 (MQ=255)                                                      acaGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCg                                                                                                >  3:316514/1‑148 (MQ=255)                                                        aGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCgtag                                                                                             <  8:331656/149‑1 (MQ=255)                                                                tctatctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTagcagc                                                                                      >  3:99346/1‑148 (MQ=255)                                                                   atctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATc                                                                                  >  5:177887/1‑149 (MQ=255)                                                                    tctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCg                                                                                 >  1:15084/1‑149 (MQ=255)                                                                    tctGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCg                                                                                 >  2:281231/1‑149 (MQ=255)                                                                       gAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                                          >  8:24668/1‑121 (MQ=255)                                                                       gAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGtt                                                                                                          <  7:24668/121‑1 (MQ=255)                                                                        aaCGACATCACCCCCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                                              <  4:254442/148‑1 (MQ=255)                                                                        aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGc                                                                                <  6:177887/146‑1 (MQ=255)                                                                        aaCGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCt                                                                              >  6:183121/1‑148 (MQ=255)                                                                                 accaccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTaaa                                                                    <  7:369680/149‑1 (MQ=255)                                                                                 accaccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTaa                                                                     >  2:316459/1‑148 (MQ=255)                                                                                  ccaccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTaaaa                                                                   <  7:113887/149‑1 (MQ=255)                                                                                   caccGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAg                                                                  >  2:403250/1‑149 (MQ=255)                                                                                    accGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAgc                                                                 <  8:391208/149‑1 (MQ=255)                                                                                                 ccATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCgg                                                         <  6:424020/144‑1 (MQ=255)                                                                                                 ccATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCgg                                                         >  5:424020/1‑144 (MQ=255)                                                                                                   aTGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGgcgc                                                  >  4:29062/1‑149 (MQ=255)                                                                                                   aTGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGgcgc                                                  >  2:89892/1‑149 (MQ=255)                                                                                                   aTGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGgcgc                                                  >  7:207638/1‑149 (MQ=255)                                                                                                    tGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGcc                                                 >  6:274037/1‑149 (MQ=255)                                                                                                         gAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTa                                               <  8:244501/146‑1 (MQ=255)                                                                                                         gAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTa                                               >  7:244501/1‑146 (MQ=255)                                                                                                         gAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATc                                             >  3:320753/1‑148 (MQ=255)                                                                                                          aGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCa                                           <  6:28624/149‑1 (MQ=255)                                                                                                          aGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCa                                           <  4:231884/149‑1 (MQ=255)                                                                                                            gTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAg                                          <  4:204425/148‑1 (MQ=255)                                                                                                               tGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                      >  1:51882/1‑149 (MQ=255)                                                                                                                 ggAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTgg                                     >  1:216801/1‑148 (MQ=255)                                                                                                                     gagTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACa                                <  3:163446/149‑1 (MQ=255)                                                                                                                          ggTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGt                             >  6:66997/1‑147 (MQ=255)                                                                                                                            tCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCaga                         <  1:281231/149‑1 (MQ=255)                                                                                                                                      cacAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGc                >  1:414246/3‑149 (MQ=255)                                                                                                                                       tCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGGAATCGACTGACGCTTAAAAATCCACACTGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc               >  1:142550/1‑149 (MQ=255)                                                                                                                                       tCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc               >  6:98611/1‑149 (MQ=255)                                                                                                                                       acAGGCTCGCGCAGGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGcc               >  2:351053/2‑149 (MQ=255)                                                                                                                                        cAGGCTCGCGCAGGAGCCTGTGATGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGACCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGAGCACCCTGCAGGCGCCTATCCAGCCTGGTAACAGGTCACAGGCGCTCGCCg              >  7:387798/1‑149 (MQ=255)                                                                                                                                           gCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                      <  2:222591/122‑1 (MQ=255)                                                                                                                                           gCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTg                                      >  1:222591/1‑122 (MQ=255)                                                                                                                                           gCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGc           <  4:316514/149‑1 (MQ=255)                                                                                                                                                  cgcgAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCt     >  8:228352/1‑148 (MQ=255)                                                                                                                                                    cgAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCTaaa  >  6:219128/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGAACATTGACAGGCTTCTTGAGGAGGCCGGTAGTCACAAGTCGAAAATCCTCACAGTGACTATCTATCTGAACGACATCACCACCGACTTCGCTGCCATGAACGAGGTTTGGGAAGAGTGGGTCGCGG‑AAGGTCAGGCTCCTGCGCGAGCCTGTGTTGAGGCGGCTATGTTCAAGCCAAGCCTGCTCGTTGAGATGTCCGTAGTAGCAGCAATCGCCTGACGCTTAAAAGCGCACCCGGCAGGCGCCTATCCAGCCTGGTAACAGGTCAGAGGCGCTCGCCGGGCACCGCTAAA  >  NC_002947/5032615‑5032909

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