Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,011,115	T→A	6.9%	L59Q (CTG→CAG)	bcsG →	endoglucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,011,115	0	T	A	6.9%	175.2 / 2.9	58	L59Q (CTG→CAG)	bcsG	endoglucanase
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base A (2/2);  total (30/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 = 6.06e-01

GCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTATCTGCTCGAGCTGCTGGGCCG  >  NC_002947/3010973‑3011259                                                                                                                                                |                                                                                                                                                 gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTa                                                                                                                                             <  2:8065/148‑1 (MQ=255) gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTAt                                                                                                                                            >  4:74125/1‑149 (MQ=255) gCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTAt                                                                                                                                            >  1:251394/1‑149 (MQ=255)           cTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCaa                                                                                                                                  <  3:74125/149‑1 (MQ=255)             ggACTGGGAGGGTGGAACCTCTCTTTTTTTGCAAAGGTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCCGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAAct                                                                                                                                <  8:253998/149‑1 (MQ=255)                   ggAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGc                                                                                                                          >  6:198772/1‑149 (MQ=255)                   ggAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGc                                                                                                                          >  2:187823/1‑149 (MQ=255)                    gAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTAGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCg                                                                                                                         >  8:57541/1‑149 (MQ=255)                      gggTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGt                                                                                                                       <  7:73333/149‑1 (MQ=255)                      gggTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGt                                                                                                                       >  7:232817/1‑149 (MQ=255)                        gTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGc                                                                                                                     >  8:242179/1‑149 (MQ=255)                         tGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGcc                                                                                                                    <  4:47352/149‑1 (MQ=255)                          ggAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTTTTCGATGCGTTTCTGCTTGTCCCTTTGCCGTGCAAGG‑TGCTGCGTATTGACCGGCAACTGATGGAGGTGca                                                                                                                    >  8:111932/1‑147 (MQ=255)                              ccTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCg                                                                                                               <  2:251394/149‑1 (MQ=255)                                 ctATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCaa                                                                                                                                  >  8:196410/1‑127 (MQ=255)                                 ctATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCaa                                                                                                                                  <  7:196410/127‑1 (MQ=255)                                  tATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTAt                                                                                                           >  5:102283/1‑149 (MQ=255)                                    tttttttGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGTCGGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTa                                                                                                            <  1:187823/147‑1 (MQ=255)                                            cAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCtttt                                                                                                  <  5:198772/148‑1 (MQ=255)                                            cAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCtttt                                                                                                  <  6:16934/148‑1 (MQ=255)                                                 tttCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAATA‑CGCAGCCACTTGCCCGGCaaagg                                                                                                                               >  3:260500/1‑111 (MQ=37)                                                 tttCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAATA‑CGCAGCCACTTGCCCGGCaaagg                                                                                                                               <  4:260500/114‑4 (MQ=37)                                                   tCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAg                                                                                           <  5:273623/148‑1 (MQ=255)                                                        tgGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACAcc                                                                                     >  5:141120/1‑149 (MQ=255)                                                          gTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTg                                                                                   >  4:120452/1‑149 (MQ=255)                                                               gcTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACctgcctgc                                                                              >  8:7067/1‑149 (MQ=255)                                                                    gCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCAccc                                                                         <  8:232817/149‑1 (MQ=255)                                                                         cTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAg                                                                    <  7:57541/149‑1 (MQ=255)                                                                          tCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCg                                                                                                               >  7:169620/1‑105 (MQ=255)                                                                          tCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCg                                                                                                               <  8:169620/105‑1 (MQ=255)                                                                               ttCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAATA‑CGCAGCCACTTGCCCGGCaaa                                                                                                                                 >  1:123346/1‑81 (MQ=21)                                                                               ttCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAATA‑CGCAGCCACTTGCCCGGCaaa                                                                                                                                 <  2:123346/82‑2 (MQ=21)                                                                               ttCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCAc                                                                           >  1:33834/1‑136 (MQ=255)                                                                               ttCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCAc                                                                           <  2:33834/136‑1 (MQ=255)                                                                               ttCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCct                                                              >  5:29369/1‑149 (MQ=255)                                                                                 ccAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCct                                                              <  5:114535/147‑1 (MQ=255)                                                                                  cAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGc                                                                              >  1:33749/1‑130 (MQ=255)                                                                                  cAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGc                                                                              <  2:33749/130‑1 (MQ=255)                                                                                  cAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCctgc                                                            <  7:7067/148‑1 (MQ=255)                                                                                             gAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCgg                                                >  1:171388/1‑149 (MQ=255)                                                                                             gAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCgg                                                >  7:130100/1‑149 (MQ=255)                                                                                                ccTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCggggg                                             >  4:30786/1‑149 (MQ=255)                                                                                                 cTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCt                                                                                >  7:104119/1‑113 (MQ=255)                                                                                                 cTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCt                                                                                <  8:104119/113‑1 (MQ=255)                                                                                                         cGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGAc                                     <  3:120452/148‑1 (MQ=255)                                                                                                         cGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGAc                                     <  6:102283/148‑1 (MQ=255)                                                                                                         cGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACt                                    >  1:138662/1‑149 (MQ=255)                                                                                                            tGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCt                                                                                <  4:179203/102‑1 (MQ=255)                                                                                                            tGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTtct                                 >  6:214025/1‑149 (MQ=255)                                                                                                             gCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCTTATTGCCCGGCAACTGCTGGCGGGGCCAATCGTTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTtctc                                <  7:111932/149‑1 (MQ=255)                                                                                                               gTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTTTCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGt                              >  6:149495/1‑149 (MQ=255)                                                                                                                     gCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAAACGGTATCGCGCTTTTCTACCAGCACACCTCGCTCCCACCCTACAGCCGCCTGCTGGCGCAGCCGGGGCTGCTGGACTTCTCGTGGGAc                         >  6:85182/1‑148 (MQ=255)                                                                                                                         gTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTAt                      <  5:214025/147‑1 (MQ=255)                                                                                                                         gTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTAt                      <  6:141120/147‑1 (MQ=255)                                                                                                                         gTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTATCt                    >  5:12466/1‑149 (MQ=255)                                                                                                                             cTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACt                        >  1:108200/1‑141 (MQ=255)                                                                                                                                     ggCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTATCTGCTCGAgctgct        <  8:130100/149‑1 (MQ=255)                                                                                                                                           t‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTATCTGCTCGAGCTGCTGGGCCg  >  2:51277/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                 GCCCGCTGGCCTGGACTGGGAGGGTGGAACCTCTATTTTTTTGCAAAGTTTCTGCTGGTCGCGCTGGGCATGCTCGACTTCCAGGCGCTGCCGAACCTGTTGTTCGCTGCGTTTCTGCTTGTCCCTTTGCCGGGCAAGT‑GGCTGCGTATTGCCCGGCAACTGCTGGCGGTGCCAATCGGTATCGCGCTTTTCTACCAGGACACCTGGCTGCCACCCTTCAGCCGCCTGCTGGCGCAGCCGGGGGTGCTGGACTTCTCGTGGGACTATCTGCTCGAGCTGCTGGGCCG  >  NC_002947/3010973‑3011259

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑