Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,277,934	G→C	9.7%	V34V (GTG→GTC)	glcG →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,277,934	0	G	C	9.7%	161.1 / 9.7	62	V34V (GTG→GTC)	glcG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/31);  new base C (4/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.29e-01

CCACTGGGTCGAAATCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAGAGATGGTCAATGG  >  NC_002947/4277786‑4278074                                                                                                                                                       |                                                                                                                                                ccACTGGGTCGAAATCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑g                                                                                                                                                 <  2:143621/149‑1 (MQ=255)  cACTGGGTCGAAATCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑gg                                                                                                                                                >  4:344975/1‑149 (MQ=255)   aCTGGGTCGAAATCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTAGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGt                                                                                                                                               <  4:315964/149‑1 (MQ=255)       ggTCGAAATCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑cg                                                                                                                                           <  2:328150/147‑1 (MQ=255)            aaaTCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑CGT‑‑CGACCAcgg                                                                                                                                    >  4:366307/1‑149 (MQ=255)             aaTCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGCTCCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGAcgg                                                                                                                                    <  1:293513/149‑1 (MQ=255)              aTCGTCGAAGCAGCCTTGAACCCGGTGCCCGCCCATGCTTAGCAAGTTCGTTATGGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGA‑CCGCCAGTAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGAcggcg                                                                                                                                  <  2:185277/149‑1 (MQ=255)              aTCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGAcggcg                                                                                                                                  >  5:84864/1‑149 (MQ=255)              aTCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGAcggcg                                                                                                                                  <  3:268181/149‑1 (MQ=255)                 gtcgAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑CGT‑‑CGACCAcggcgat                                                                                                                                <  3:366307/148‑3 (MQ=255)                  tcgAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGC‑CGT‑‑CGTCGACcacggcga                                                                                                                              <  5:355765/149‑3 (MQ=255)                  tcgAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGC‑CGT‑‑CGTCGACcacggcga                                                                                                                              >  6:355765/1‑147 (MQ=255)                           ccTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGc                                                                                                                     >  8:144801/1‑149 (MQ=255)                           ccTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGc                                                                                                                     <  3:148051/149‑1 (MQ=255)                                     ggAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGc                                                                                                           <  1:25598/149‑1 (MQ=255)                                      gAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGcc                                                                                                          <  1:46868/149‑1 (MQ=255)                                           cgcccATGCTTAGCAAGTTCGTTATCGGCCAGTCCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTCGCCCTGGAACGCCTGGAc                                                                                                     >  2:169243/1‑149 (MQ=255)                                              ccATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGc                                                                                                  >  2:375466/1‑149 (MQ=255)                                                    ttAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGTCGGTAACGACAACGGCCGCCACCCACTGGCCCTGGAACGCCTGGACGGCTgc                                                                                               >  1:316168/1‑149 (MQ=255)                                                            ttCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAg                                                                                   >  4:217963/1‑149 (MQ=255)                                                               gTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGc                                                                                <  3:332209/149‑1 (MQ=255)                                                                ttATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGcc                                                                               <  1:237774/149‑1 (MQ=255)                                                                  aTCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCt                                                                              <  2:42628/148‑1 (MQ=255)                                                                   tCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTCGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTAc                                                                            <  1:169243/149‑1 (MQ=255)                                                                      gCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGGCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATc                                                                         >  1:268108/1‑149 (MQ=255)                                                                       ccAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGCCGGC‑GCGCAGCCGTCCAg                                                                                   <  5:383376/138‑1 (MQ=255)                                                                       ccAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGCCGGC‑GCGCAGCCGTCCAg                                                                                   >  6:383376/1‑138 (MQ=255)                                                                       ccAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATcg                                                                        >  5:338427/1‑149 (MQ=255)                                                                         aGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATcgcg                                                                      >  1:362443/1‑149 (MQ=255)                                                                                aaTCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATcgcg                                                                      >  5:357242/1‑142 (MQ=255)                                                                                aaTCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATcgcg                                                                      <  6:357242/142‑1 (MQ=255)                                                                                    gCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGc                                                           >  8:91026/1‑149 (MQ=255)                                                                                     cccGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑gg                                                                                                                                                >  5:208853/1‑66 (MQ=255)                                                                                     cccGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑gg                                                                                                                                                <  6:208853/66‑1 (MQ=255)                                                                                        gcgTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGc                                                       >  5:337030/1‑149 (MQ=255)                                                                                                aCCGCCGC‑CCGCCAGGAAGCCCAGGCCCGAGAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGCCGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCAccgccg                                                >  1:35362/1‑148 (MQ=255)                                                                                                    ccgc‑ccgcCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTg                                           >  7:155599/1‑149 (MQ=255)                                                                                                    ccgc‑ccgcCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTg                                           >  3:262171/1‑149 (MQ=255)                                                                                                     cgc‑ccgccAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTgg                                          >  1:83441/1‑149 (MQ=255)                                                                                                         ccgccAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGgc                                       <  4:262171/149‑1 (MQ=255)                                                                                                         ccgccAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGgc                                       <  2:268108/149‑1 (MQ=255)                                                                                                          cgccAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGgcg                                      <  4:303125/149‑1 (MQ=255)                                                                                                          cgccAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGgcg                                      >  3:303125/1‑149 (MQ=255)                                                                                                           gccAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGgcgc                                     <  3:217963/149‑1 (MQ=255)                                                                                                              aGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCaaa                                  <  2:362443/149‑1 (MQ=255)                                                                                                                  aGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAgaga                              <  6:84864/149‑1 (MQ=255)                                                                                                                   gCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAc                             <  3:344975/149‑1 (MQ=255)                                                                                                                   gCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAc                             <  2:316168/149‑1 (MQ=255)                                                                                                                     ccAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGGCGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAcgc                           >  6:954/1‑149 (MQ=255)                                                                                                                     ccAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAcgc                           <  1:375466/149‑1 (MQ=255)                                                                                                                     ccAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAcgc                           <  7:182239/149‑1 (MQ=255)                                                                                                                       aGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCgg                         <  7:91026/149‑1 (MQ=255)                                                                                                                         gCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGGCGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGa                       <  5:954/149‑1 (MQ=255)                                                                                                                         gCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGa                       <  6:338427/149‑1 (MQ=255)                                                                                                                             gACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGGAAAGAGCCGCGGGATTAt                   >  3:255436/1‑149 (MQ=255)                                                                                                                                aaTGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGaa                <  6:337030/149‑1 (MQ=255)                                                                                                                                aaTGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGaa                <  7:144801/149‑1 (MQ=255)                                                                                                                                         tGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAgaga            >  5:207779/1‑144 (MQ=255)                                                                                                                                         tGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAgaga            <  6:207779/144‑1 (MQ=255)                                                                                                                                           aCCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACACCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAGAGATGGTCaa     >  7:35901/1‑149 (MQ=255)                                                                                                                                            tggccGCCGT‑CGT‑‑CGACCACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAGAGATGGTCAAt    >  8:109376/5‑149 (MQ=255)                                                                                                                                             cATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAGAGATGGTCAAtg   >  7:19449/1‑149 (MQ=255)                                                                                                                                              gccGCCGT‑CGT‑‑CGACCACGGCCCCCACCCACAGGCCCTGGAACGCCTGCACGGCAGCGC‑GCCCGCCAGAGCCTACATCGCGCAGGAAAAGGCCCGCACCGCCGCACAGGCGCGCAAACAGACGCGGGAATATGAAGAGAAGGTCAAtgg  >  8:328439/3‑149 (MQ=255)                                                                                                                                              aTCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCTCTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAGAGATGGTCAAtgg  >  7:290739/1‑149 (MQ=255)                                                                                                                                                  ccGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAcgc                           >  1:65710/1‑120 (MQ=255)                                                                                                                                                  ccGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGAcgc                           <  2:65710/120‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                CCACTGGGTCGAAATCGTCGAAGCAGCCTTGAACCCGGAGCCCGCCCATGCTTAGCAAGTTCGTTATCGGCCAGACCGAAATCGCCCGCGTACTGACCGCCGC‑CCGCCAGGAAGCCCAGGCCCGACAATGGAGCGTGACCATCGCCGT‑GGT‑‑CGACGACGGCGGCCACCCACTGGCCCTGGAACGCCTGGACGGCTGCGC‑GCCGGCCAGTGCCTACATCGCGCTGGAAAAGGCCCGCACCGCCGCACTGGGGCGCAAAGAGACGCGGGATTATGAAGAGATGGTCAATGG  >  NC_002947/4277786‑4278074

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