Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,610,954	G→A	5.8%	G910G (GGC→GGT)	mcoA ←	Mn(II) copper oxidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,610,954	0	G	A	5.8%	207.9 / 2.2	69	G910G (GGC→GGT)	mcoA	Mn(II) copper oxidase A
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base A (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGA  >  NC_002947/3610810‑3611096                                                                                                                                                 |                                                                                                                                               cAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACTCACGccat                                                                                                                                            >  3:304177/1‑149 (MQ=255) cAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccat                                                                                                                                            >  4:353277/1‑149 (MQ=255)  aGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatc                                                                                                                                           >  1:162683/1‑149 (MQ=255)    cACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatccc                                                                                                                                         <  3:355225/149‑1 (MQ=255)       gTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatcccatc                                                                                                                                      >  7:409689/1‑149 (MQ=255)       gTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatc                                                                                                                                           >  4:46499/1‑144 (MQ=255)       gTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGccatc                                                                                                                                           <  3:46499/144‑1 (MQ=255)          cAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCt                                                                                                                                    <  2:162683/148‑1 (MQ=255)          cAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTg                                                                                                                                   >  2:5422/1‑149 (MQ=255)           aGGGTGATCGGCTTGCATGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTg                                                                                                                                   <  8:3838/148‑1 (MQ=255)               tGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGc                                                                                                                              <  7:237351/149‑1 (MQ=255)               tGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGc                                                                                                                              <  2:165593/149‑1 (MQ=255)                 aTCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCgg                                                                                                                            >  6:267360/1‑149 (MQ=255)                   cGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATAGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCgggg                                                                                                                          >  3:313052/1‑149 (MQ=255)                   cGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATAGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATACCATACTGGCAGCgggg                                                                                                                          >  2:219650/1‑149 (MQ=255)                     gCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTg                                                                                                                        <  3:381142/149‑1 (MQ=255)                       ttGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTggg                                                                                                                      >  2:247100/1‑149 (MQ=255)                           cTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCAt                                                                                                                   >  2:286553/1‑148 (MQ=255)                             gCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCAGCCTGGCCGCGGGGTGGGCATTAc                                                                                                                >  5:18957/1‑149 (MQ=255)                              cAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTAcc                                                                                                               <  8:409689/149‑1 (MQ=255)                               aGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCCGGCCTACCGAGGCCATGTAACGAACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCt                                                                                                              >  5:214134/1‑149 (MQ=255)                                   gTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCgggg                                                                                                                          >  7:350568/1‑133 (MQ=255)                                   gTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCgggg                                                                                                                          <  8:350568/133‑1 (MQ=255)                                   gTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGa                                                                                                           <  4:135028/148‑1 (MQ=255)                                    ttGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACACCATCCCATCCTGGCAGCGGGGTGGGCATTAccgtgacc                                                                                                         >  4:361005/1‑143 (MQ=255)                                    ttGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACACCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAAc                                                                                                         >  4:354943/1‑149 (MQ=255)                                       cGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGc                                                                                                      <  4:313052/149‑1 (MQ=255)                                               tCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATTCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCggg                                                                                                <  1:304217/147‑1 (MQ=255)                                                cGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGgtgt                                                                                             >  7:123579/1‑149 (MQ=255)                                                      ccACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAgg                                                                                       >  6:373283/1‑149 (MQ=255)                                                      ccACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAgg                                                                                       >  4:145684/1‑149 (MQ=255)                                                       cACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGt                                                                                      <  3:158801/149‑1 (MQ=255)                                                        aCGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTc                                                                                     <  5:373283/149‑1 (MQ=255)                                                         cGACAGGTGGTTTGTTGGCCTGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCGATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTcc                                                                                    <  1:219650/149‑1 (MQ=255)                                                         cGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTcc                                                                                    >  5:29816/1‑149 (MQ=255)                                                          gACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACGTGGAACTGCCCCGGGTGTTCAAGGTccc                                                                                   <  1:247100/149‑1 (MQ=255)                                                               gtggtTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAgcg                                                                               >  1:56079/1‑148 (MQ=255)                                                                   tttGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCgcagc                                                                           <  7:377893/148‑1 (MQ=255)                                                                      gttgGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCgcagcagca                                                                       >  5:179462/1‑149 (MQ=255)                                                                       ttgGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTg                                                                                                       >  2:269145/1‑116 (MQ=255)                                                                       ttgGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTg                                                                                                       <  1:269145/116‑1 (MQ=255)                                                                         ggCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATc                                                                     >  1:9919/1‑148 (MQ=255)                                                                           ccGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGaa                                                                  <  5:267360/149‑1 (MQ=255)                                                                           ccGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGaa                                                                  <  1:164866/149‑1 (MQ=255)                                                                             gggTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAaga                                                                <  6:29816/149‑1 (MQ=255)                                                                             gggTCGGCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGGGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAaga                                                                <  1:299394/149‑1 (MQ=255)                                                                                  gtcatcgtcgtcatcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCt                                                           <  1:286553/149‑1 (MQ=255)                                                                                          cgtcatcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACACCATCCCATGCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTgg                                                   <  3:361005/149‑1 (MQ=255)                                                                                           gtcatcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTgg                                                   <  1:29623/148‑1 (MQ=255)                                                                                           gtcatcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACACCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTgg                                                   <  3:354943/148‑1 (MQ=255)                                                                                               tcgGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGtgtt                                              >  3:242306/1‑149 (MQ=255)                                                                                                       ggTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGt                                      >  7:186393/1‑149 (MQ=255)                                                                                                        gTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGCGTCCTCGTGCTGGGTGTTGTGGCAggg                                     >  8:67920/1‑147 (MQ=255)                                                                                                           cGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTc                                  >  8:254261/1‑149 (MQ=255)                                                                                                            ggAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTcc                                 <  6:365796/149‑1 (MQ=255)                                                                                                            ggAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTcc                                 <  6:179462/149‑1 (MQ=255)                                                                                                             gAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCa                                >  4:165237/1‑149 (MQ=255)                                                                                                             gAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCa                                >  7:35756/1‑149 (MQ=255)                                                                                                                  ttGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAg                           <  8:123579/149‑1 (MQ=255)                                                                                                                      cAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGc                       <  4:242306/149‑1 (MQ=255)                                                                                                                        ggCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCt                     <  2:9919/149‑1 (MQ=255)                                                                                                                        ggCCTACAGAGGCCATGTAACGCACGCCATCACATCCTGGCAGCGGGGTTGGCATTACCTGGAACTGCCCGGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGTGTTTTGTGGCAGTGCTCCATGTAGGTGCCt                     <  4:304177/149‑1 (MQ=255)                                                                                                                            tACCGAGGCCATGTAACGCACGCCATCCCATCTTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGc                   <  5:276018/147‑1 (MQ=255)                                                                                                                            tACCGAGGCCATGTAACGCACGCCATCCCATCTTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGc                   >  6:276018/1‑147 (MQ=255)                                                                                                                              ccGAGGCCATGTAACGCACGCCATCCCTTCCTGGCAGCGGGGTGGGCATTAGCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCTTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGaa                <  6:214134/148‑1 (MQ=255)                                                                                                                              ccGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAAc               >  1:436261/1‑149 (MQ=255)                                                                                                                                      aTGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGaa       <  3:353277/149‑1 (MQ=255)                                                                                                                                       tGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAg      <  8:186393/149‑1 (MQ=255)                                                                                                                                         tAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGCGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGc     <  3:145684/148‑1 (MQ=255)                                                                                                                                           acgcacgcCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGa  >  7:148633/1‑149 (MQ=255)                                                                                                                                           acgcacgcCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGa  >  3:134075/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                               CAGCACGTTCAGGGTGATCGGCTTGCCTGCAGTCGTTGCGGCACTGTCGTTGGCCACGACAGGTGGTTTGTTGGCCGGGTCGTCATCGTCGTCATCGGTTTTGGTGCGGAACGTTGGCAGGCCTACCGAGGCCATGTAACGCACGCCATCCCATCCTGGCAGCGGGGTGGGCATTACCTGGAACTGCCCCGGGTGTTCAAGGTCCCAGCGCAGCAGCATCGAAGAGTCCTCGTGCTGGGTGTTGTGGCAGTGCTCCATGTAGGTGCCTGCGAACTCGCGGAAGCGGA  >  NC_002947/3610810‑3611096

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