Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	396,856	C→G	7.6%	L30L (CTG→CTC)	PP_0329 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	396,856	0	C	G	7.6%	277.5 / 8.2	88	L30L (CTG→CTC)	PP_0329	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/41);  new base G (2/5);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGAT  >  NC_002947/396709‑397003                                                                                                                                                     |                                                                                                                                                      gTATCC‑AAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑a                                                                                                                                                     <  4:293977/148‑1 (MQ=255)      cAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGAtt                                                                                                                                                >  4:143733/1‑149 (MQ=255)       aaaGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGAtt                                                                                                                                                >  8:41752/1‑148 (MQ=255)         aGGTGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTT‑CAATCCTG‑CGGATT‑ggaa                                                                                                                                           <  7:281596/149‑4 (MQ=255)         aGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑Gaa                                                                                                                                            <  7:229195/149‑1 (MQ=255)         aGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTT‑CAATCCTG‑CGGATT‑ggaa                                                                                                                                           <  4:203986/149‑4 (MQ=255)         aGGGGGTGTAGGTGTGATTTTTTTGC‑TTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTT‑CAATCCTG‑CGGATT‑ggaac                                                                                                                                          <  8:200898/149‑5 (MQ=255)          gggggTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAAc                                                                                                                                           <  6:162296/149‑1 (MQ=255)          gggggTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTT‑CAATCCTG‑CGGATT‑ggaac                                                                                                                                          <  7:2110/149‑5 (MQ=255)             ggTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCt                                                                                                                                        >  7:27002/1‑149 (MQ=255)               tgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTc                                                                                                                                      <  2:272811/149‑1 (MQ=255)               tgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTc                                                                                                                                      <  2:42729/149‑1 (MQ=255)               tgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTc                                                                                                                                      <  4:5148/149‑1 (MQ=255)               tgtAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTc                                                                                                                                      >  1:40047/1‑149 (MQ=255)                       tgATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCCGATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGCGAGCGTTCCAATCCGC‑AGCATT‑GAACGCTTCTGAACGt                                                                                                                               >  1:273056/1‑149 (MQ=255)                       tgATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGtt                                                                                                                              >  6:126596/1‑149 (MQ=255)                       tgATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGtt                                                                                                                              <  5:126596/149‑1 (MQ=255)                          tttttttGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTc                                                                                                                            <  2:341026/148‑1 (MQ=255)                              tttGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAAt                                                                                                                        <  8:27002/148‑1 (MQ=255)                                 gCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCTAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGAc                                                                                                                     <  3:207956/148‑1 (MQ=255)                                  cTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAg                                                                                                                   >  5:340759/1‑149 (MQ=255)                                   tttCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGt                                                                                                                  >  5:341505/1‑149 (MQ=255)                                        gTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCa                                                                                                             <  6:340759/149‑1 (MQ=255)                                            ggCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGttt                                                                                                         >  5:136306/1‑149 (MQ=255)                                                tGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctg                                                                                                     <  2:369830/149‑1 (MQ=255)                                                  cTCACGACAAGCACTGAGTCC‑ATCTAATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGGCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGTTCGTTTCCAATGACAGTGGTCAGTTTGctggt                                                                                                   <  3:143733/149‑1 (MQ=255)                                                   tCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGtt                                                                                                                              <  2:332244/121‑1 (MQ=255)                                                   tCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGtt                                                                                                                              >  1:332244/1‑121 (MQ=255)                                                   tCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                                       <  4:218559/144‑1 (MQ=255)                                                   tCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGc                                                                                                       >  3:218559/1‑144 (MQ=255)                                                     aCGCCAAGCACTGAGTCC‑ATCTCCTCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGG‑AGGCTT‑GAACGCGGCTGATAGTGTCGAATGACAGTGGTCCGTTTGctggtctg                                                                                                >  6:292217/1‑149 (MQ=255)                                                       gACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtctgg                                                                                               >  8:6399/1‑148 (MQ=255)                                                        aCAAGCACTCCGTCC‑AGCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCGGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTa                                                                                             <  8:58547/149‑1 (MQ=255)                                                          aatcACTGAGTCA‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGATTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCAGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                                  >  7:171635/4‑142 (MQ=255)                                                          aaGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                                  <  8:166752/142‑1 (MQ=255)                                                          aaGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctggtc                                                                                                  >  7:166752/1‑142 (MQ=255)                                                          aaGCAATGATTCC‑ATCTCATCTGGTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGAGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTACAATGACAGTGGTCAGTTTGctggtc                                                                                                  <  8:171635/142‑1 (MQ=255)                                                            gCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGc                                                                                                                                                       >  3:357550/1‑89 (MQ=255)                                                            gCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGc                                                                                                                                                       <  4:357550/89‑1 (MQ=255)                                                            gCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGaaa                                                                                         >  1:31413/1‑149 (MQ=255)                                                                   gTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  <  2:31413/149‑1 (MQ=255)                                                                   gTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  <  7:41752/149‑1 (MQ=255)                                                                   gTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  <  7:6399/149‑1 (MQ=255)                                                                   gTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  <  3:239870/149‑1 (MQ=255)                                                                   gTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  <  6:341505/149‑1 (MQ=255)                                                                           tcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGATCTGGTAGAAATATTTgg                                                                                  <  4:299705/142‑1 (MQ=255)                                                                           tcatctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGATCTGGTAGAAATATTTgg                                                                                  >  3:299705/1‑142 (MQ=255)                                                                             atctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTg                                                                                                                 <  4:286590/109‑1 (MQ=255)                                                                             atctGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTg                                                                                                                 >  3:286590/1‑109 (MQ=255)                                                                                  ttGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTg                                                                    >  8:129253/1‑149 (MQ=255)                                                                                      gAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGtttt                                                                >  8:262457/1‑149 (MQ=255)                                                                                        aaaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGTAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGttttt                                                               <  2:40047/148‑1 (MQ=255)                                                                                         aaGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCa                                                             <  8:306906/149‑1 (MQ=255)                                                                                          aGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCAc                                                            <  7:129253/149‑1 (MQ=255)                                                                                             tGGGAACTTCCTGCTTCTTTGTCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctgg                                                                                                    <  5:261316/106‑1 (MQ=255)                                                                                             tGGGAACTTCCTGCTTCTTTGTCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGctgg                                                                                                    >  6:261316/1‑106 (MQ=255)                                                                                                 aaCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTGCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  >  1:117863/1‑120 (MQ=255)                                                                                                 aaCTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgg                                                                                  <  2:117863/120‑1 (MQ=255)                                                                                                    ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                                  >  2:16698/1‑149 (MQ=255)                                                                                                    ttCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCAGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTAAGTTTGCTGGTCTGGTGGAAATATTTGGGGGTATAAGGACTTTTTTTCACCGACGTGCTc                                                  >  7:77266/1‑149 (MQ=255)                                                                                                      ccTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTccc                                                >  5:357378/1‑149 (MQ=255)                                                                                                        tGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATGTTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGt                                              >  2:103262/1‑149 (MQ=255)                                                                                                             cttTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCt                                         >  6:341747/1‑149 (MQ=255)                                                                                                               ttGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTg                                        >  7:38256/1‑148 (MQ=255)                                                                                                                   cAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGtt                                   >  6:32507/1‑149 (MQ=255)                                                                                                                    aGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGCAAGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGtt                                   >  1:48857/1‑149 (MQ=255)                                                                                                                           tCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCg                            >  6:164712/1‑148 (MQ=255)                                                                                                                             aGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                                >  6:106270/1‑94 (MQ=255)                                                                                                                             aGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTgggg                                                                                <  5:106270/94‑1 (MQ=255)                                                                                                                              gCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                        <  2:273056/149‑1 (MQ=255)                                                                                                                               cAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                                  >  4:195084/1‑122 (MQ=255)                                                                                                                               cAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTc                                                  <  3:195084/122‑1 (MQ=255)                                                                                                                                aaGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCTGTTTGCTGGTCTGGTAGAAATATTTGGAGGTTTAAGGAAGTATTTTCACGGACTTGGACCCTTTGAGTCCATTTGGCCCAGaaaaaa                      >  8:97166/1‑147 (MQ=255)                                                                                                                                aaGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTTGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                        >  1:118585/1‑147 (MQ=255)                                                                                                                                aaGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGaaa                        >  8:137754/1‑147 (MQ=255)                                                                                                                                 agagGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAg                     <  6:136306/149‑1 (MQ=255)                                                                                                                                   agGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGAt                   <  7:137754/149‑1 (MQ=255)                                                                                                                                    gggAGCGTTCCCATCCTG‑AGGTGT‑GAACGCTTCTGATCGTTTCGTATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCTGGACTGTTTTTCACGGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGAt                   <  8:288213/148‑1 (MQ=255)                                                                                                                                    gggAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGt                                                                              >  2:281509/1‑89 (MQ=255)                                                                                                                                      gAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGAGAGTGGTCAGTTTGCTGGTCTGGTTGATATATTTTGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTTGACCCGGAAATAGATgg                 >  7:332494/1‑148 (MQ=255)                                                                                                                                      aaGCGTT‑CAATCCTG‑CGGATTGGAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGc                >  2:60098/2‑149 (MQ=255)                                                                                                                                       aGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCa                                      <  8:292276/126‑1 (MQ=255)                                                                                                                                           ttCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTggg           <  5:341747/149‑1 (MQ=255)                                                                                                                                            tCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTggg           >  2:219949/1‑148 (MQ=255)                                                                                                                                              cAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCAc        >  2:21055/1‑149 (MQ=255)                                                                                                                                              cAATCCCC‑AGGATT‑GAACCCTTCTGATCCTTTCCAATGACAGTGGTCAGTATGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTAACCGACGAGCTCACGTTCTCTCCAGTTGGCCCCGGAAATAGATGGCCTgccccc        >  2:100785/1‑144 (MQ=255)                                                                                                                                               aaTCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACt       <  6:357378/149‑1 (MQ=255)                                                                                                                                                   cGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCTCGGAAATAGATGGCCTGGGCACTACGa   <  1:16698/149‑1 (MQ=255)                                                                                                                                                   cGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGa   <  7:262457/149‑1 (MQ=255)                                                                                                                                                    gC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGAt  <  2:118585/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTATCCAAAGGGGGTGTAGGTGTGATTTTTTTGCTTTCCGTCAGGCGTGCTCACGACAAGCACTGAGTCC‑ATCTCATCTGTTGCGAAAAGCTGGGAACTTCCTGCTTCTTTGCCAGCGGACTCAGCAAGAGGGAGCGTTCCAATCCGC‑AGGATT‑GAACGCTTCTGATCGTTTCCAATGACAGTGGTCAGTTTGCTGGTCTGGTAGAAATATTTGGGGGTATCAGGACTGTTTTTCACCGACGTGCTCCCGTTGTCTGCAGTTGGCCCCGGAAATAGATGGCCTGGGCACTACGAT  >  NC_002947/396709‑397003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑