Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,440,014	G→A	8.1%	A271V (GCC→GTC)	amn ←	AMP nucleosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,440,014	0	G	A	8.1%	175.9 / 3.2	61	A271V (GCC→GTC)	amn	AMP nucleosidase
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (3/2);  total (33/28)
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.97e-01

GCCGATGGTCTGCGACTGCCGCAGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGGTCGTCACGC  >  NC_002947/5439869‑5440157                                                                                                                                                  |                                                                                                                                                 gCCGATGGTCTGCGACTGCCGCAGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTcggcgg                                                                                                                                               >  2:114041/1‑149 (MQ=255)     aTGGTCTGCGACTGCCGCCGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGGCGGTGATGTTCTTGGCGCTGGATGCGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGCCGACGGCGAt                                                                                                                                           <  4:206148/149‑1 (MQ=255)      tGGTCTGCGACTGCCGCAGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATc                                                                                                                                          <  6:77828/149‑1 (MQ=255)      tGGTCTGCGACTGCCGCAGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGCCGACGGCGAtg                                                                                                                                          <  7:8086/149‑2 (MQ=255)         tCTGCGACTGCCGCAGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGGCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAgg                                                                                                                                       <  2:216123/149‑1 (MQ=255)                        ggCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCAACGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATc                                                                                                                        <  4:193296/149‑1 (MQ=255)                         gccgccGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCt                                                                                                                       <  4:44598/149‑1 (MQ=255)                                     aCCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGcc                                                                                                           <  6:382/149‑1 (MQ=255)                                       cGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGc                                                                                                            >  7:122959/1‑146 (MQ=255)                                          tcatcaGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGAcc                                                                                                      <  2:19090/149‑1 (MQ=255)                                            atcaGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCAc                                                                                                    >  3:264625/1‑149 (MQ=255)                                            atcaGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCAc                                                                                                    >  6:222147/1‑149 (MQ=255)                                              cagccagcAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGc                                                                                                  >  1:224630/1‑149 (MQ=255)                                               agccagcAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCt                                                                                                 >  6:19216/1‑149 (MQ=255)                                                     ccaTGCGGGGGCACCCAGGCCAGGTGGGTGGGGTATTTCTTGGCCGTGGATTGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGTTAGCGGTGCCAGACCACGCTGG‑ccac                                                                                          <  3:249414/147‑1 (MQ=255)                                                     cAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑ccac                                                                                          >  1:312229/1‑149 (MQ=255)                                                          cGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑ccacca                                                                                        >  1:197512/1‑146 (MQ=255)                                                          cGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑ccacca                                                                                        <  2:197512/146‑1 (MQ=255)                                                           gggCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGc                                                                                    <  3:56021/149‑1 (MQ=255)                                                             gcgcAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACGACGCTGG‑CCACCAGGGCCt                                                                                  >  8:226091/1‑149 (MQ=255)                                                              cgcAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTg                                                                                 <  8:278820/149‑1 (MQ=255)                                                               gcagcaCGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTgg                                                                                <  8:122959/149‑1 (MQ=255)                                                               gcagcaCGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGAGCAGGTGGTAGGCCGGCAGCTGGTAGCGGTGCCAGACCCCGCTGG‑CCACCAGGGCCTgg                                                                                >  4:177269/1‑149 (MQ=255)                                                                 agcaCGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGt                                                                              >  7:13947/1‑149 (MQ=255)                                                                 agcaCGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGt                                                                              >  1:38284/1‑149 (MQ=255)                                                                        ccAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCa                                                                       <  5:19216/149‑1 (MQ=255)                                                                          aGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATg                                                                     <  4:264625/149‑1 (MQ=255)                                                                              ggtCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTcc                                                                 <  8:13947/149‑1 (MQ=255)                                                                                            ttGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGc                                                                                    >  1:155191/1‑116 (MQ=255)                                                                                            ttGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGc                                                                                    <  2:155191/116‑1 (MQ=255)                                                                                              ggCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATc                                                  <  3:92171/148‑1 (MQ=255)                                                                                              ggCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATc                                                  <  1:114041/148‑1 (MQ=255)                                                                                                cGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCacg                                               <  2:224630/149‑1 (MQ=255)                                                                                                    ggATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGt                                           >  4:99743/1‑149 (MQ=255)                                                                                                     gATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGtt                                          >  8:46164/1‑149 (MQ=255)                                                                                                           cccACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCgg                                    <  3:99743/149‑1 (MQ=255)                                                                                                           cccACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCgg                                    <  5:270964/149‑1 (MQ=255)                                                                                                           cccACACCAATGTTGACCAAGGTAACGCAATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCgg                                    >  8:266411/1‑149 (MQ=255)                                                                                                            ccACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCtttt                                                       >  4:87253/1‑129 (MQ=255)                                                                                                            ccACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCtttt                                                       <  3:87253/129‑1 (MQ=255)                                                                                                             cacaCCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCgg                                    >  4:210850/1‑147 (MQ=255)                                                                                                               cacCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGcc                                                                        >  5:175475/1‑109 (MQ=255)                                                                                                               cacCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGcc                                                                        <  6:175475/109‑1 (MQ=255)                                                                                                                acCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGgcagca                               >  7:14271/1‑149 (MQ=255)                                                                                                                  cAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCAc                              >  3:185694/1‑148 (MQ=255)                                                                                                                   aaTGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATAAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCATGGCCTGGGTTTCGACATGGTCCATGCCCGTTTCGATCACGACGTTGCCCGGCAGCACCa                            >  5:9069/1‑149 (MQ=255)                                                                                                                     tGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGTCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCAta                          <  5:222147/149‑1 (MQ=255)                                                                                                                         gACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCAtacgta                      >  1:170327/1‑149 (MQ=255)                                                                                                                         gACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCAtacgt                       <  2:312229/148‑1 (MQ=255)                                                                                                                            cAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCAtacgtacg                    >  1:39186/1‑148 (MQ=255)                                                                                                                                 tAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCtttt                                                       >  3:317419/1‑108 (MQ=255)                                                                                                                                 tAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCtttt                                                       <  4:317419/108‑1 (MQ=255)                                                                                                                                  aaCGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGcgcg             >  4:169404/1‑149 (MQ=255)                                                                                                                                  aaCGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGcgcg             >  3:170830/1‑149 (MQ=255)                                                                                                                                   aCGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCgg            >  4:110273/1‑149 (MQ=255)                                                                                                                                    cctgaTCGCCGCCGACGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCGCCATACGTACGAAGCGCggg           >  5:122051/5‑149 (MQ=255)                                                                                                                                     ctgaTCGCCGCCGACGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGgt          >  1:188191/4‑149 (MQ=255)                                                                                                                                      tgaTCGCCGCCGACGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGGACCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGgt          >  1:97379/3‑149 (MQ=255)                                                                                                                                       cATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGgtcg        >  5:237320/1‑149 (MQ=255)                                                                                                                                        aTCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGgtcgt       >  2:162158/1‑149 (MQ=255)                                                                                                                                        aTCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGgtcgt       >  4:36986/1‑149 (MQ=255)                                                                                                                                             cGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGGTCGTCACGc  <  7:266411/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 GCCGATGGTCTGCGACTGCCGCAGGCCGCCGCAGTGACCGATCATCAGCCAGCAATGCGGGCGCAGCACGGCCAGGTGGTCGGTGATGTTCTTGGCGTTGGATGGGCCCACACCAATGTTGACCAAGGTAACGCCATCGCCGTCGGCGGCGATCAGGTGGTAGGCCGGCATCTGGTAGCGGTGCCAGACCACGCTGG‑CCACCAGGGCCTGGGTTTCGCCATGGTCCATGCCCTTTTCGATCACGACGTTGCCCGGCAGCACCATACGTACGAAGCGCGGGTCGTCACGC  >  NC_002947/5439869‑5440157

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