Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,786,987	T→C	8.5%	S134P (TCG→CCG) ‡	trxC →	thioredoxin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,786,987	0	T	C	8.5%	165.3 / 5.1	59	S134P (TCG→CCG) ‡	trxC	thioredoxin
Reads supporting (aligned to +/- strand):  ref base T (25/29);  new base C (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 = 4.04e-01

GCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTT  >  NC_002947/5786839‑5787128                                                                                                                                                     |                                                                                                                                              gCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGt                                                                                                                                               >  2:234700/1‑149 (MQ=255)     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)                                                                                                                                                     |                                                                                                                                              GCTTGCTGGCCGCTGCCGCCTGGCCAAGCTCGACAGCGAAGCCAACCGTAACCTCGCAGGGCAACTGGGCATACGCTCGATCCCCAGCCTGCTGCTGTTCAAGAACGGCCGCGAAATCAGCCGCCAGGCCGGGGCGTTTCCACTCCAGTCGCTGCTGGAGTGGCTGCGTAGCCAAGGGGTTTAAGCGTTTTCTTCAAGCAGCGAATGCAGCTCGACAAACTGCTGAGTCAGCTTGTGCCGCGGCTCCAGGTGAATCAGCGGCAGGCTGGCGTGGTGCGATTCGCGCATTT  >  NC_002947/5786839‑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: ‑