Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,549	A→C	intergenic (‑154/‑123)	metJ ← / → metB	transcriptional repressor, S‑adenosylmethionine‑binding/cystathionine gamma‑synthase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,549	0	A	C	100.0%	271.5 / NA	68	intergenic (‑154/‑123)	metJ/metB	transcriptional repressor, S‑adenosylmethionine‑binding/cystathionine gamma‑synthase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (30/38);  total (30/38)

CTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTATTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCG  >  NC_000913/4128402‑4128696                                                                                                                                                    |                                                                                                                                                    cTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCttt                                                                                                                                                  >  2:154538/1‑151 (MQ=255)    aaTCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTggg                                                                                                                                               >  1:506752/1‑151 (MQ=255)       ccTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGAtt                                                                                                                                            <  1:634461/151‑1 (MQ=255)       ccTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGAtt                                                                                                                                            <  1:576069/151‑1 (MQ=255)       ccTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGAt                                                                                                                                             <  1:1177731/150‑1 (MQ=255)        ctctTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGAttt                                                                                                                                           >  2:659705/1‑151 (MQ=255)        ctctTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGAttt                                                                                                                                           >  1:271934/1‑151 (MQ=255)        ctcctcGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGAttt                                                                                                                                           >  2:1040325/1‑151 (MQ=255)             cGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTCCCGTGACATCGTGTAATGCGCCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCa                                                                                                                                      <  2:616673/151‑1 (MQ=255)                cAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATc                                                                                                                                   >  2:378368/1‑151 (MQ=255)                 aaTAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCt                                                                                                                                  >  1:476550/1‑151 (MQ=255)                  aTAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCt                                                                                                                                  >  2:65875/1‑150 (MQ=255)                    aaaTTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCtat                                                                                                                                >  1:334840/1‑150 (MQ=255)                         gagaCCAGACCACACAGTTGATGTGGGTACTGACCGTAGACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCa                                                                                                                           <  1:1196619/150‑1 (MQ=255)                             ccagaccaCACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAaagaa                                                                                                                      <  2:805513/151‑1 (MQ=255)                              cagaccaCACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAaagaag                                                                                                                     <  2:119350/151‑1 (MQ=255)                              cagaccaCACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAaaga                                                                                                                       >  2:126085/1‑149 (MQ=255)                              cagaccaCACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAaaga                                                                                                                       >  2:297850/1‑149 (MQ=255)                                    acacaGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGttt                                                                                                                  >  1:989037/1‑148 (MQ=255)                                     cacaGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGAt                                                                                                              <  1:479435/151‑1 (MQ=255)                                       caGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGt                                                                                                            <  2:911795/151‑1 (MQ=255)                                        aGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGt                                                                                                            <  2:308312/150‑1 (MQ=255)                                           tGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCa                                                                                                         <  2:134779/150‑1 (MQ=255)                                              tgtgGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGAt                                                                                                      <  1:576767/150‑1 (MQ=255)                                                        gACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGt                                                                                            <  1:164927/150‑1 (MQ=255)                                                        gACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGt                                                                                            <  1:543039/150‑1 (MQ=255)                                                        gACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGt                                                                                            <  1:1206221/150‑1 (MQ=255)                                                         aCCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTcc                                                                                          <  1:839478/151‑1 (MQ=255)                                                              aaaCCCGCATAGTTTACCGTACAGGCGTTACTGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTaa                                                                                     <  1:300210/151‑1 (MQ=255)                                                                    gCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAAt                                                                               <  1:176628/151‑1 (MQ=255)                                                                     cATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATg                                                                              <  1:1094560/151‑1 (MQ=255)                                                                     cATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATg                                                                              >  1:324491/1‑151 (MQ=255)                                                                        aGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGttt                                                                           <  1:308842/151‑1 (MQ=255)                                                                        aGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGtt                                                                            <  2:322827/150‑1 (MQ=255)                                                                             aCCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTActct                                                                      >  2:345965/1‑151 (MQ=255)                                                                             aCCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTActct                                                                      <  2:820031/151‑1 (MQ=255)                                                                                     ggCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTg                                                              >  1:606537/1‑151 (MQ=255)                                                                                     ggCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTg                                                              <  2:211917/151‑1 (MQ=255)                                                                                      gCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTg                                                              >  1:207032/1‑150 (MQ=255)                                                                                      gCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTg                                                              >  2:1017691/1‑150 (MQ=255)                                                                                          tACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACAtt                                                         >  2:1097074/1‑151 (MQ=255)                                                                                           aCCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACAttt                                                        <  1:473133/151‑1 (MQ=255)                                                                                            ccGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACAttt                                                        >  2:906620/1‑150 (MQ=255)                                                                                             cGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTc                                                       >  2:534207/1‑150 (MQ=255)                                                                                                  cATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGAc                                                 >  1:612883/1‑151 (MQ=255)                                                                                                     cGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACaa                                               >  2:1062994/1‑150 (MQ=255)                                                                                                         tAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGccc                                          >  1:977740/1‑151 (MQ=255)                                                                                                           aTGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAg                                        <  2:430498/151‑1 (MQ=255)                                                                                                            tGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGt                                                                   >  2:821050/1‑123 (MQ=255)                                                                                                            tGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGt                                                                   <  1:821092/123‑1 (MQ=255)                                                                                                             gCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAgg                                       <  1:630233/150‑1 (MQ=255)                                                                                                             gCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAgg                                       <  1:786906/150‑1 (MQ=255)                                                                                                               aCCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGaa                                    <  1:1108910/151‑1 (MQ=255)                                                                                                               aCCTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGa                                     >  2:628728/1‑150 (MQ=255)                                                                                                                ccTGTCGGCGTGATAATGCATAGAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGACAAGGGAAc                                   <  1:479461/151‑1 (MQ=255)                                                                                                                 cTGTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACGATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACt                                  >  2:1160258/1‑151 (MQ=255)                                                                                                                   gTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACtt                                 <  2:982492/150‑1 (MQ=255)                                                                                                                   gTCGGCGTGATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACt                                  <  1:42542/149‑1 (MQ=255)                                                                                                                           gATAATGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACAt                         <  2:372376/150‑1 (MQ=255)                                                                                                                              aaTGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGAcg                     >  1:311447/1‑151 (MQ=255)                                                                                                                                tGCATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGAcgcg                   >  1:495665/1‑151 (MQ=255)                                                                                                                                  cATATAATTTTAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTa                 <  1:426018/151‑1 (MQ=255)                                                                                                                                            tAACGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGcca        <  1:1184094/150‑1 (MQ=255)                                                                                                                                             aaCGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGccacc      >  1:790102/1‑151 (MQ=255)                                                                                                                                             aaCGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGccacc      >  2:725663/1‑151 (MQ=255)                                                                                                                                               cGGCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCAt    <  2:495641/151‑1 (MQ=255)                                                                                                                                                ggCTCTTTGGGATTTGCGCAATCTATACGCAAAGAAGTGTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATc   <  2:747811/151‑1 (MQ=255)                                                                                                                                                 gCTCTTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCg  <  1:63964/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCATATAATTTTAACGGCTATTTGGGATTTGCTCAATCTATACGCAAAGAAGTTTAGATGTCCAGATGTATTGACGTCCATTAACACAATGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCG  >  NC_000913/4128402‑4128696

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