Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,298,035	T→C	E35G (GAA→GGA)	yraQ ←	putative inner membrane permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,298,035	0	T	C	100.0%	255.2 / NA	71	E35G (GAA→GGA)	yraQ	putative inner membrane permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/33);  total (38/33)

CGCTTTCCATACCGCGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTTCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATA  >  NC_000913/3297899‑3298165                                                                                                                                         |                                                                                                                                   cgcTTTCCATACCGCGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTccc                                                                                                                                  <  1:178260/139‑1 (MQ=255)     ttCCATACCGCGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                              <  1:1125318/139‑1 (MQ=255)              gcgAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATAcc                                                                                                                     >  2:706365/1‑139 (MQ=255)                               cgcgTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACa                                                                                                    >  1:861994/1‑139 (MQ=255)                                    aaTCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACgag                                                                                               >  1:681578/1‑139 (MQ=255)                                      tCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGagaa                                                                                             >  1:493024/1‑139 (MQ=255)                                            gccgccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGaaaa                                                                                       >  2:208531/1‑139 (MQ=255)                                            gccgccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGaaaa                                                                                       >  2:706937/1‑139 (MQ=255)                                                ccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                <  2:3490/94‑1 (MQ=255)                                                ccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                >  1:3490/1‑94 (MQ=255)                                                ccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAgcgc                                                                                   >  2:143131/1‑139 (MQ=255)                                                ccTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAgcgc                                                                                   <  1:355433/139‑1 (MQ=255)                                                       tGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTc                                                                            <  1:208531/139‑1 (MQ=255)                                                        gggTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTcc                                                                           >  1:902276/1‑139 (MQ=255)                                                               cATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTg                                                                    >  2:1103715/1‑139 (MQ=255)                                                                   cGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATg                                                                >  2:911274/1‑139 (MQ=255)                                                                    gCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTccc                                                                                                                                  >  2:80907/1‑72 (MQ=255)                                                                    gCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTccc                                                                                                                                  <  1:80907/72‑1 (MQ=255)                                                                    gCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATgg                                                               >  1:158177/1‑139 (MQ=255)                                                                    gCCTGCGCAAGGATAGATTCACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACAAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATgg                                                               >  2:416669/1‑139 (MQ=255)                                                                         cgcAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGtt                                                          >  1:305709/1‑139 (MQ=255)                                                                              ggATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTgccgc                                                     <  2:622132/139‑1 (MQ=255)                                                                              ggATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTgccgc                                                     <  2:1027840/139‑1 (MQ=255)                                                                                   gATTTACTGATACTGTGGGTTTCGGCAGCAGTAAACGCGTTGCCGTAGTAAGGTCCCCATTTCACATACAAGAGACAGGCAATGACAACGAGAAAGAAAAGCGAGGGTGTGCACCACTGAATGGGCGTTGGCGCGTgag                                                >  2:784969/1‑139 (MQ=255)                                                                                     tttACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGAt                                              <  2:861994/139‑1 (MQ=255)                                                                                        aCCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                              <  2:717447/56‑1 (MQ=255)                                                                                        aCCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCAtt                                                                                                                              >  1:717447/1‑56 (MQ=255)                                                                                        aCCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACaa                                                                                                   <  1:724114/83‑1 (MQ=255)                                                                                        aCCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACaa                                                                                                   >  2:724114/1‑83 (MQ=255)                                                                                         ccGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAg                                          >  2:166132/1‑139 (MQ=255)                                                                                         ccGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAg                                          <  1:409209/139‑1 (MQ=255)                                                                                              aCTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGa                                     <  2:305709/139‑1 (MQ=255)                                                                                               cTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  2:158177/139‑1 (MQ=255)                                                                                               cTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  2:493024/139‑1 (MQ=255)                                                                                               cTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  1:143131/139‑1 (MQ=255)                                                                                               cTGTGCGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAAAGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  1:706365/139‑1 (MQ=255)                                                                                                tgtgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAcc                                   <  1:1103715/139‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTc                                                                            >  2:1029995/1‑96 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTc                                                                            <  1:1029995/96‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGa                                                                   <  1:250054/105‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGa                                                                   >  2:250054/1‑105 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATg                                                                >  2:1059105/1‑108 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATg                                                                <  1:1059105/108‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                         <  1:1027161/131‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                         >  2:935827/1‑131 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                         >  2:1027161/1‑131 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGa                                         <  1:935827/131‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  2:161791/136‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  1:14944/136‑1 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    >  1:161791/1‑136 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    >  2:14944/1‑136 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAg                                 >  1:423167/1‑139 (MQ=255)                                                                                                  tgGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAAAGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAg                                 >  2:663661/1‑139 (MQ=255)                                                                                                   gggTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTccac                                                                         >  2:241761/1‑98 (MQ=255)                                                                                                   gggTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTccac                                                                         <  1:241761/98‑1 (MQ=255)                                                                                                    ggTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                >  2:528261/1‑42 (MQ=255)                                                                                                    ggTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCa                                                                                                                                <  1:528261/42‑1 (MQ=255)                                                                                                       ttCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATa                            >  1:22249/1‑139 (MQ=255)                                                                                                          ggCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCa                         <  2:902276/139‑1 (MQ=255)                                                                                                          ggCAGCAGTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCa                         >  1:327939/1‑139 (MQ=255)                                                                                                              gcagTAAACGCTTTGCCGTAGTAAGGTCCCCATTTca                                                                                                                           >  2:325442/1‑37 (MQ=255)                                                                                                              gcagTAAACGCTTTGCCGTAGTAAGGTCCCCATTTca                                                                                                                           <  1:325442/37‑1 (MQ=255)                                                                                                               cagTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATAcc                                                                                                                     <  1:515611/42‑1 (MQ=255)                                                                                                               cagTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATAcc                                                                                                                     >  2:515611/1‑42 (MQ=255)                                                                                                               cagTAAACGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGCTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTccc                    >  2:548952/1‑139 (MQ=255)                                                                                                                    aaCGCTTTGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGa               <  1:663661/139‑1 (MQ=255)                                                                                                                         tttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTcccc                   <  2:204716/130‑1 (MQ=255)                                                                                                                         tttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTcccc                   >  1:204716/1‑130 (MQ=255)                                                                                                                          ttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    >  2:399141/1‑112 (MQ=255)                                                                                                                          ttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGAc                                    <  1:399141/112‑1 (MQ=255)                                                                                                                          ttGCCGTAGTAAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGt         <  2:681578/139‑1 (MQ=255)                                                                                                                           tGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                >  2:1069269/1‑35 (MQ=255)                                                                                                                           tGCCGTAGTAAGGTCCCCATTTCACATACCAgaga                                                                                                                <  1:1069269/35‑1 (MQ=255)                                                                                                                                 agtaAGGTCCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATa  >  1:1031248/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                   CGCTTTCCATACCGCGAGGAAATAGATCATCGCGTAATCCAACGCCGCCTGCCATGGGTTAGCATCCGCCTGCGCAAGGATAGATTTACCGATACTGTGGGTTTCGGCAGCAGTAAACGCTTTGCCGTAGTAAGGTTCCCATTTCACATACCAGAGACCGGCAATGACAACGAGAAAGAAAAGCGCGGGTTTCCACCACTGAATGGGCGTTGCCGCCTGAGATGAAGACTGACCAGTCATAGCATTCCCCAGGAATAAGTTATGATA  >  NC_000913/3297899‑3298165

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