Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,786,988	C→A	8.5%	S134* (TCG→TAG) ‡	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,786,988	0	C	A	8.5%	179.3 / 7.0	59	S134* (TCG→TAG) ‡	trxC	thioredoxin
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base A (2/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

GCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTT  >  NC_002947/5786843‑5787128                                                                                                                                                  |                                                                                                                                             gCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCg                                                                                                                                             <  1:364713/147‑1 (MQ=255) gCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCg                                                                                                                                             <  2:99973/147‑1 (MQ=255)     gCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGctgctg                                                                                                                                       >  4:317016/1‑149 (MQ=255)     gCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGctgctg                                                                                                                                       <  6:27980/149‑1 (MQ=255)      ccGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGctgctg                                                                                                                                       >  8:122508/1‑148 (MQ=255)       cGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGa                                                                                                                                     <  8:256958/149‑1 (MQ=255)        gCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAg                                                                                                                                    <  5:340183/149‑1 (MQ=255)         cTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                   <  1:317455/149‑1 (MQ=255)         cTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGCAGCGACTGGAGt                                                                                                                                   <  5:258052/149‑1 (MQ=255)            ccgccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTg                                                                                                                                  <  8:155572/147‑1 (MQ=255)             cgccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCt                                                                                                                               <  8:275683/149‑1 (MQ=255)              gccTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTg                                                                                                                              >  5:162574/1‑149 (MQ=255)               ccTGGCCAAGCTAGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGATGTTCAAGAACTGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCAATCCAGTCGCTGCTGGAGTGGCTGc                                                                                                                             <  4:237532/149‑1 (MQ=255)                cTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCg                                                                                                                            >  4:97828/1‑149 (MQ=255)                             aCAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGtt                                                                                                               >  7:303818/1‑149 (MQ=255)                              cAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGttt                                                                                                              <  7:122508/149‑1 (MQ=255)                              cAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGttt                                                                                                              <  4:123911/149‑1 (MQ=255)                                  gAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGGCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGc                                                                                                          <  7:73782/149‑1 (MQ=255)                                  gAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGc                                                                                                          <  8:303818/149‑1 (MQ=255)                                     gCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGtt                                                                                                       <  5:382147/149‑1 (MQ=255)                                        aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttc                                                                                                    >  1:373065/1‑149 (MQ=255)                                        aaCCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttc                                                                                                    >  2:101609/1‑149 (MQ=255)                                           cGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttctt                                                                                                  >  6:95324/1‑148 (MQ=255)                                            gTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTttcttc                                                                                                 >  2:33678/1‑148 (MQ=255)                                                     cAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGaa                                                                                       <  7:283854/149‑1 (MQ=255)                                                     cAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGaa                                                                                       <  7:156425/149‑1 (MQ=255)                                                     cAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGaa                                                                                       <  6:162574/149‑1 (MQ=255)                                                      aGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAAt                                                                                      >  3:159547/1‑149 (MQ=255)                                                         gCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCgctgct                                                                                                                                        >  3:323241/1‑96 (MQ=255)                                                         gCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCgctgct                                                                                                                                        <  4:323241/96‑1 (MQ=255)                                                         gCAACTGGGCATACGCTCGATCCCCACCCTGCAGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGGGAATGCa                                                                                   <  6:108801/149‑1 (MQ=255)                                                                 gaaTACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGAct                                                                           <  5:262150/147‑2 (MQ=255)                                                                 gCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGAc                                                                            >  4:190227/1‑148 (MQ=255)                                                                      cGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActg                                                                      <  3:317016/149‑1 (MQ=255)                                                                        cTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAActgct                                                                    <  3:97828/149‑1 (MQ=255)                                                                           gATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAg                                                                 >  7:183016/1‑149 (MQ=255)                                                                                ccAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGc                                                            >  7:9479/1‑149 (MQ=255)                                                                                 cAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCt                                                           >  4:78821/1‑149 (MQ=255)                                                                                             gTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCt                                               >  5:13820/1‑149 (MQ=255)                                                                                              ttCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTc                                              >  7:93981/1‑149 (MQ=255)                                                                                               tCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTc                                              >  2:293671/1‑148 (MQ=255)                                                                                                 aaGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAg                                           >  7:321376/1‑149 (MQ=255)                                                                                                         ccGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAAt                                     >  8:100408/1‑147 (MQ=255)                                                                                                            cgAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCg                                <  4:159547/149‑1 (MQ=255)                                                                                                              aaaTCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCagcagc                              <  8:321376/149‑1 (MQ=255)                                                                                                                aTCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAg                            <  1:234700/149‑1 (MQ=255)                                                                                                                              cGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCCGGTGAATCAGCGGCAGGCTGGCGTGGTGCg              <  1:33678/149‑1 (MQ=255)                                                                                                                               ggggCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGa             <  8:9479/149‑1 (MQ=255)                                                                                                                                gggCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCACGGTGAATCAGCGGCAGGCTGGCGTGGTGCGAt            >  6:81165/1‑149 (MQ=255)                                                                                                                                 tacgcagccACTCCAGCAGCGACTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGAtt           >  4:362807/8‑149 (MQ=255)                                                                                                                                    gTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc        >  4:356254/1‑149 (MQ=255)                                                                                                                                    gTTTCAACTCCCGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTATGCGTTTTCTTCCATCACCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCGGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgc        <  2:373065/149‑1 (MQ=255)                                                                                                                                     tttCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgcg       >  6:93493/1‑149 (MQ=255)                                                                                                                                      ttCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTcgcg       >  2:21521/1‑148 (MQ=255)                                                                                                                                        ccACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCa     <  8:183016/148‑1 (MQ=255)                                                                                                                                        ccACTCCAGCAGCGACTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAt    >  3:263040/1‑149 (MQ=255)                                                                                                                                          aCTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCAttt  <  1:101609/149‑1 (MQ=255)                                                                                                                                              cAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCa                                                                                                >  6:208195/1‑51 (MQ=255)                                                                                                                                              cAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCa                                                                                                <  5:208195/51‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             GCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTT  >  NC_002947/5786843‑5787128

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