Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,456,311	T→G	8.6%	E51A (GAG→GCG)	eco ←	ecotin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,456,311	0	T	G	8.6%	159.4 / 5.5	56	E51A (GAG→GCG)	eco	ecotin
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base G (3/2);  total (28/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 = 8.89e-01

GCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCG  >  NC_002947/3456167‑3456459                                                                                                                                                  |                                                                                                                                                     gccgggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACtcgt                                                                                                                                                   <  5:367472/148‑1 (MQ=255)     ggccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCAAGGGAGTGTGCTTCCAGCTTGCCGCCCAGGCGCTGACGATTGCAGTCGCCCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGcc                                                                                                                                              >  2:56750/1‑149 (MQ=255)      gccgCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCt                                                                                                                                             >  3:33832/1‑149 (MQ=255)          cTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGc                                                                                                                                         <  3:364644/149‑1 (MQ=255)              cGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTTTGTGCTTCCATGTTGCCGCCCAGGCGCTGACGATTGCTGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGTCCGACTCGTCGGCCTGTGcctg                                                                                                                                      <  6:8583/148‑3 (MQ=255)                 tGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAgg                                                                                                                                  >  2:369077/1‑149 (MQ=255)                   tCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGt                                                                                                                                 <  2:369530/148‑1 (MQ=255)                       ggCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATa                                                                                                                            <  1:369077/149‑1 (MQ=255)                        gCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATaa                                                                                                                           >  1:28532/1‑149 (MQ=255)                              gtaATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGc                                                                                                                      <  4:33832/148‑1 (MQ=255)                                 aTTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCggg                                                                                                                  >  4:252490/1‑149 (MQ=255)                                  ttGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGt                                                                                                                 <  2:28532/149‑1 (MQ=255)                                         cccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCt                                                                                                          >  7:64706/1‑149 (MQ=255)                                          ccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCtt                                                                                                         >  3:295282/1‑149 (MQ=255)                                          ccAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCtt                                                                                                         >  4:280950/1‑149 (MQ=255)                                           cAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGc                                                                                                                                     <  7:338873/120‑1 (MQ=255)                                           cAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGc                                                                                                                                     >  8:338873/1‑120 (MQ=255)                                           cAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCttt                                                                                                        <  3:325390/149‑1 (MQ=255)                                                     gggTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCg                                                                                              >  8:34820/1‑149 (MQ=255)                                                        tgtgtgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTATATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCgggg                                                                                           >  7:192602/1‑149 (MQ=255)                                                            tgCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAgg                                                                                       >  2:125513/1‑149 (MQ=255)                                                              cTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgc                                                                                     >  3:329713/1‑149 (MQ=255)                                                               ttCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGCGAGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgc                                                                                     >  1:72293/1‑149 (MQ=255)                                                                 ccAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCTGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgcgc                                                                                   >  4:144524/1‑148 (MQ=255)                                                                 ccAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGgcgc                                                                                   <  4:31176/148‑1 (MQ=255)                                                                  cAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCg                                                                                              >  6:224278/1‑136 (MQ=255)                                                                  cAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCg                                                                                              <  5:224278/136‑1 (MQ=255)                                                                  cAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGcgcg                                                                                  <  7:34820/148‑1 (MQ=255)                                                                   aGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGAt                                                                                <  1:89822/149‑1 (MQ=255)                                                                      ttGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTc                                                                             <  2:291664/149‑1 (MQ=255)                                                                        gccgccCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTcc                                                                            <  3:252490/148‑1 (MQ=255)                                                                         ccgccCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCtt                                                                          >  6:383940/1‑149 (MQ=255)                                                                          cgccCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTc                                                                         <  1:56750/149‑1 (MQ=255)                                                                               aGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCt                                                                    >  4:155964/1‑149 (MQ=255)                                                                                   gcTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTggcg                                                                <  1:153945/149‑1 (MQ=255)                                                                                     tGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCTGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGcgcg                                                                                  <  7:243710/129‑1 (MQ=255)                                                                                     tGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGcgcg                                                                                  >  8:243710/1‑129 (MQ=255)                                                                                               aGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGt                                                    >  3:54182/1‑149 (MQ=255)                                                                                                                     catagAGGATTTCCAGTTTATAGGCAGACTCGTCGGACTGGGAAGGCAGGTTGATAACCTGCAGGGTTAAGCCTTTTTCCGCATAGGGGTAGAGCGCTATGTCCTTCAGGCTGGAGGCCATCGAAGGTGCGGAGGCGGCCAGGCTCAgg                               >  8:224494/5‑149 (MQ=255)                                                                                                                     cAGCGAGGAATTCCAGTTTATAGGCCGACGCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAgg                               <  5:191973/149‑1 (MQ=255)                                                                                                                      aGCGAGGATTTCCAGTTTATAGGCTGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAggg                              <  3:144524/149‑1 (MQ=255)                                                                                                                       gCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCCGCCGGGTGAAGCGTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGc                             >  4:98322/1‑149 (MQ=255)                                                                                                                                 ccAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAgg                                     >  5:258766/1‑131 (MQ=255)                                                                                                                                 ccAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAgg                                     <  6:258766/131‑1 (MQ=255)                                                                                                                                  cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATcgcg                   <  1:125513/148‑1 (MQ=255)                                                                                                                                  cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCgg                  <  5:383940/149‑1 (MQ=255)                                                                                                                                  cAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCgg                  >  3:183089/1‑149 (MQ=255)                                                                                                                                     ggcccaGGCCGACGAGTCGGCCTAGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTgcggc                                               <  2:71039/112‑1 (MQ=255)                                                                                                                                       cccaGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCAtt             >  6:154328/4‑149 (MQ=255)                                                                                                                                        ccaGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCa               >  2:22616/3‑146 (MQ=255)                                                                                                                                         taGGCCGACTCGTCGGCCTGGGCCGGCAGGTTGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGTCGTCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTg            <  7:224494/148‑1 (MQ=255)                                                                                                                                         taGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTg            <  2:72293/148‑1 (MQ=255)                                                                                                                                         taGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTg            <  3:280950/148‑1 (MQ=255)                                                                                                                                          aGGCCGACTCGTCGGCCTGTGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTggg          >  7:207178/1‑149 (MQ=255)                                                                                                                                          aGGCCGACGAGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTggg          >  3:20826/1‑149 (MQ=255)                                                                                                                                           ggCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTgggg         <  4:54182/149‑1 (MQ=255)                                                                                                                                           ggCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTggg          <  8:192602/148‑1 (MQ=255)                                                                                                                                                  tcgtcgGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGgcg  >  8:246837/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCCGGGCCGCTGACGTTGTCCAGGCGGTAGTAATTGTAGCCCCAGCCTTCCAGGGTGTGTGCTTCCAGGTTGCCGCCCAGGCGCTGACGATTGCAGTCGACCTGC‑AGCGTCTTGCCAGCGAGGATTTCCAGTTTATAGGCCGACTCGTCGGCCTGGGCCGGCAGGTGGATAACCTGCCGGGTGAAGCCTTTTTCCGCCTCGGGGTAGGGCGCGATGTCCTTCAGGCTGGCGGCCATCGCAGGTGCGGCGGCGGCCAGGCTCAGGGCCAGAATCGCGGTCATTGGGGTAGGGCG  >  NC_002947/3456167‑3456459

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