Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,033,217	G→A	G25R (GGG→AGG)	trkH →	potassium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,033,217	0	G	A	100.0%	233.0 / NA	60	G25R (GGG→AGG)	trkH	potassium transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/28);  total (32/28)

GGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTC  >  NC_000913/4033074‑4033355                                                                                                                                                |                                                                                                                                           ggTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGt                                                                                                                                     >  1:596306/1‑151 (MQ=255)  gTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGt                                                                                                                                     >  1:617896/1‑150 (MQ=255)  gTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTa                                                                                                                                    <  2:606999/151‑1 (MQ=255)   ttCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGt                                                                                                                                     >  1:726933/1‑149 (MQ=255)    tCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAg                                                                                                                                   <  1:534139/150‑1 (MQ=255)         gCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCa                                                                                                                             <  2:313650/151‑1 (MQ=255)                  aGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCggg                                                                                                                    >  1:23209/1‑151 (MQ=255)                    cGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGAt                                                                                                                  <  2:260003/151‑1 (MQ=255)                    agATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGAt                                                                                                                  >  1:349050/2‑151 (MQ=255)                                  tCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCgcgc                                                                                                    >  1:146237/1‑151 (MQ=255)                                  tCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCgcgc                                                                                                    >  2:824072/1‑151 (MQ=255)                                       aCTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTCCCGGCATGGAGCGGGCCGCTCctttc                                                                                               <  1:881389/151‑2 (MQ=255)                                       aCTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCtttt                                                                                                >  2:543766/1‑150 (MQ=255)                                         ttCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTAcc                                                                                             >  1:538131/1‑151 (MQ=255)                                              tttGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGAc                                                                                        >  2:530125/1‑151 (MQ=255)                                              tttGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGAc                                                                                        >  2:14227/1‑151 (MQ=255)                                                    gAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAg                                                                                          >  1:313291/1‑143 (MQ=255)                                                    gAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAg                                                                                          <  2:313338/143‑1 (MQ=255)                                                    gAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCttttt                                                                                  >  2:657792/1‑151 (MQ=255)                                                     aaCTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCtttttt                                                                                 <  1:517517/151‑1 (MQ=255)                                                         aaGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTc                                                                              <  2:152430/150‑1 (MQ=255)                                                                cGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGc                                                                      >  2:251991/1‑151 (MQ=255)                                                                  gCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCa                                                                    <  1:328125/151‑1 (MQ=255)                                                                    agagATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCAtt                                                                  <  1:880671/151‑1 (MQ=255)                                                                    agagATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCAtt                                                                  <  2:409861/151‑1 (MQ=255)                                                                      agaTGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTg                                                                 >  1:566380/1‑150 (MQ=255)                                                                         tGCATTTTCGCGCCATTACACGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGctc                                                             >  1:73986/1‑151 (MQ=255)                                                                           cATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTa                                                           <  2:503143/151‑1 (MQ=255)                                                                             ttttCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATg                                                         >  1:795605/1‑151 (MQ=255)                                                                               ttCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCt                                                       >  1:558607/1‑151 (MQ=255)                                                                                 cgcgCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCtgt                                                     <  2:521851/151‑1 (MQ=255)                                                                                 cgcgCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCtgt                                                     >  2:773993/1‑151 (MQ=255)                                                                                  gcgcCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCtgtg                                                    <  1:718593/151‑1 (MQ=255)                                                                                  gcgcCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCtgtg                                                    >  2:813156/1‑151 (MQ=255)                                                                                        ttACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGcc                                              >  1:135180/1‑151 (MQ=255)                                                                                            ccGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAAc                                          >  2:30403/1‑151 (MQ=255)                                                                                                aTCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGc                                       >  1:408363/1‑150 (MQ=255)                                                                                                  cGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCaaa                                    <  2:17747/151‑1 (MQ=255)                                                                                                  cGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCaa                                     <  1:145787/150‑1 (MQ=255)                                                                                                         ctactGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGaa                               >  2:850020/1‑149 (MQ=255)                                                                                                           actGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGc                           <  1:28900/151‑1 (MQ=255)                                                                                                           actGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGc                           <  2:369566/151‑1 (MQ=255)                                                                                                            ctGGTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCg                          >  1:709826/1‑151 (MQ=255)                                                                                                              ggTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCg                          <  2:491770/149‑1 (MQ=255)                                                                                                              ggTCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGaa                        <  1:672136/151‑1 (MQ=255)                                                                                                                tCATCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACt                      >  2:15010/1‑151 (MQ=255)                                                                                                                  aTCTTATTTTCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACtt                     >  1:296148/1‑150 (MQ=255)                                                                                                                      tagtttCAGTGGCCAGGAGTACCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGTGCTTTTACCCAGACCTTATTTGTCGCCCTCGCCATTGGCTCGATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAAATTAAAt                 <  1:151302/147‑1 (MQ=255)                                                                                                                          ttCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCg             <  2:596382/150‑1 (MQ=255)                                                                                                                          ttCAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCg             <  2:23208/150‑1 (MQ=255)                                                                                                                            cAGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTg           >  1:638512/1‑150 (MQ=255)                                                                                                                             aGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAg         <  2:521248/151‑1 (MQ=255)                                                                                                                             aGGGACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAg         <  1:808597/151‑1 (MQ=255)                                                                                                                              gggACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAgg        <  2:423924/151‑1 (MQ=255)                                                                                                                               agACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAggg       >  1:320243/2‑151 (MQ=255)                                                                                                                                gACCATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAgggg      <  2:727023/151‑1 (MQ=255)                                                                                                                                   cATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGTAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGttt   >  1:3176/1‑151 (MQ=255)                                                                                                                                   cATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGttt   <  1:530059/151‑1 (MQ=255)                                                                                                                                   cATGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGttt   >  2:5457/1‑151 (MQ=255)                                                                                                                                    aTGATTATCCCTAGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTc  >  1:212081/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                           GGTTCATTGCAATTGTTAGCGATTGAAGAATAATCCCCACTTCGTTTTGCAGAACTAAGGAAGCGGCAGAGATGCATTTTCGCGCCATTACCCGAATCGTTGGACTACTGGTCATCTTATTTTCAGGGACCATGATTATCCCTGGGCTGGTAGCACTCATCTACCGGGATGGAGCGGGCCGCGCTTTTACCCAGACCTTTTTTGTCGCCCTCGCCATTGGCTCTATGCTGTGGTGGCCGAACCGCAAAGAGAAAGGCGAACTTAAATCCCGTGAGGGGTTTC  >  NC_000913/4033074‑4033355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑