Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,013,361	T→C	R103R (CGT→CGC)	metE →	cobalamin‑independent homocysteine transmethylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,013,361	0	T	C	100.0%	250.8 / NA	70	R103R (CGT→CGC)	metE	cobalamin‑independent homocysteine transmethylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (41/29);  total (41/29)

CTGCCGGTGGGCGATTTTGCCTGGTACGATCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGTGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGA  >  NC_000913/4013227‑4013495                                                                                                                                       |                                                                                                                                       ctgcCGGTGGGCGATTTTGCCTGGTACGATCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGcc                                                                                                                                   >  2:79067/1‑140 (MQ=255) ctgcCGGTGGGCGATTTTGCCTGGTACGATCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGcc                                                                                                                                   >  1:273576/1‑140 (MQ=255)       gTGGGCGATTTTGCCTGGTACGATCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTgg                                                                                                                             >  2:65164/1‑140 (MQ=255)           gCGATTTTGCCTGGTACGATCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGaa                                                                                                                         >  2:554162/1‑140 (MQ=255)                         gaCGATCATGTACTGACCACCAGGCTGCTGCGGGGTAAAGCTCCGGCGCGTCATCAGAACAAAGGTTGTTCGGTAGATATCGACACCCGGTTCCGTATTGGTCGTGGGCGCTCGCCGACTGGCGAACCTGCGGCGTCAGc                                                                                                           >  2:651836/2‑140 (MQ=255)                           cGATCATGTACTGACCTCCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCgg                                                                                                         <  1:258694/140‑1 (MQ=255)                              tCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGaaa                                                                                                      <  2:579499/140‑1 (MQ=255)                                  gTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGAc                                                                                                  >  2:236343/1‑140 (MQ=255)                                      tGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              <  1:79067/140‑1 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              >  2:978152/1‑138 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              <  1:978152/138‑1 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  2:513396/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  2:642852/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  2:138462/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  2:683678/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  2:137644/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  1:778189/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  1:573472/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  1:570332/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  1:459418/1‑140 (MQ=255)                                        accaccAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATg                                                                                            >  1:309907/1‑140 (MQ=255)                                             caGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              >  2:338890/1‑133 (MQ=255)                                             caGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              >  2:658368/1‑133 (MQ=255)                                             caGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              <  1:338890/133‑1 (MQ=255)                                             caGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCaaa                                                                                              <  1:658368/133‑1 (MQ=255)                                             caGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTa                                                                                       >  2:210707/1‑140 (MQ=255)                                              aGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTaa                                                                                      <  2:459418/140‑1 (MQ=255)                                                 ctgctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAacac                                                                                   <  1:137644/140‑1 (MQ=255)                                                  tgctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGa                                                                                                   >  1:641924/1‑123 (MQ=255)                                                  tgctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGa                                                                                                   <  2:641924/123‑1 (MQ=255)                                                  tgctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACAcc                                                                                  <  1:65164/140‑1 (MQ=255)                                                  tgctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACAcc                                                                                  <  2:309907/140‑1 (MQ=255)                                                    ctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCaa                                                                                <  2:573472/140‑1 (MQ=255)                                                    ctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCaa                                                                                <  1:236343/140‑1 (MQ=255)                                                    ctgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCaa                                                                                >  2:221839/1‑140 (MQ=255)                                                     tgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAaca                                                                                    <  1:176076/135‑1 (MQ=255)                                                     tgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAaca                                                                                    >  2:176076/1‑135 (MQ=255)                                                     tgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAAc                                                                               <  2:273576/140‑1 (MQ=255)                                                     tgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAAc                                                                               <  1:221839/140‑1 (MQ=255)                                                     tgctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAAc                                                                               <  1:554162/140‑1 (MQ=255)                                                       ctgGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTa                                                                             <  1:210707/140‑1 (MQ=255)                                                         gggTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATc                                                                           >  1:74756/1‑140 (MQ=255)                                                          ggTAACGTTCCGGCNCGTCATCAGAACNAANATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCa                                                                          >  1:443036/1‑140 (MQ=255)                                                          ggTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTa                                                                                       >  2:88866/1‑127 (MQ=255)                                                          ggTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTa                                                                                       <  1:88866/127‑1 (MQ=255)                                                          ggTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACAcc                                                                                  <  2:291642/132‑1 (MQ=255)                                                          ggTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACAcc                                                                                  >  1:291642/1‑132 (MQ=255)                                                              aCGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTa                                                                             >  1:901271/1‑133 (MQ=255)                                                              aCGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTa                                                                             <  2:901271/133‑1 (MQ=255)                                                               cGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACa                                                                     <  2:443036/140‑1 (MQ=255)                                                               cGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACa                                                                     <  1:513396/140‑1 (MQ=255)                                                                    cGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTg                                                                >  1:519660/1‑140 (MQ=255)                                                                    cGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTg                                                                >  2:533176/1‑140 (MQ=255)                                                                           tcatcaGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGt                                                         <  2:74756/140‑1 (MQ=255)                                                                             atcaGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTc                                                       >  2:655131/1‑140 (MQ=255)                                                                                      aaGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGc                                              <  1:655131/140‑1 (MQ=255)                                                                                         aTGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCaa                                           <  2:570332/140‑1 (MQ=255)                                                                                                       atCGACACCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGAc                             <  1:683678/140‑1 (MQ=255)                                                                                                          gACTCCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTg                          <  1:148844/140‑1 (MQ=255)                                                                                                            cacCCTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGa                        >  1:588667/1‑140 (MQ=255)                                                                                                               ccTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGc                     >  2:498785/1‑140 (MQ=255)                                                                                                               ccTGTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGc                     >  2:748618/1‑140 (MQ=255)                                                                                                                  gTTCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAgc                  >  2:890105/1‑140 (MQ=255)                                                                                                                    tCCGTATTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAgctg                <  1:890105/140‑1 (MQ=255)                                                                                                                         aTTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGa           >  1:983910/1‑140 (MQ=255)                                                                                                                         aTTGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGa           >  1:499328/1‑140 (MQ=255)                                                                                                                           tGGTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACg         >  2:56210/1‑140 (MQ=255)                                                                                                                             gTCGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGGCGGAGTTCGTTTAAGGCCAACAGTTCAAACTGACCTGGACGTAGCTGCTGGACGaa       >  2:708844/1‑140 (MQ=255)                                                                                                                               cGTGGACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGt     >  2:38408/1‑140 (MQ=255)                                                                                                                                  ggACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGa  >  1:408670/1‑140 (MQ=255)                                                                                                                                  ggACGCGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGa  >  1:149056/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                       CTGCCGGTGGGCGATTTTGCCTGGTACGATCATGTACTGACCACCAGTCTGCTGCTGGGTAACGTTCCGGCGCGTCATCAGAACAAAGATGGTTCGGTAGATATCGACACCCTGTTCCGTATTGGTCGTGGACGTGCGCCGACTGGCGAACCTGCGGCGGCAGCGGAAATGACCAAATGGTTTAACACCAACTATCACTACATGGTGCCGGAGTTCGTTAAAGGCCAACAGTTCAAACTGACCTGGACGCAGCTGCTGGACGAAGTGGA  >  NC_000913/4013227‑4013495

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