BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,363,427	T→C	T309A (ACT→GCT)	dcuA ←	anaerobic dicarboxylate transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,363,427	0	T	C	100.0%	1316.6 / NA	404	T309A (ACT→GCT)	dcuA	anaerobic dicarboxylate transport
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (231/173); total (231/173)

Only 100 of 418 total aligned reads displayed.

AGAAGAAGATGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGTATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGA  >  NC_000913/4363337‑4363516

  aagaagaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCg                                                                                       >  2:5575389/1‑93 (MQ=255)

  aagaagaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCg                                                                                       <  1:7152731/93‑1 (MQ=255)

      agaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCa                                                                                   >  1:1246368/1‑93 (MQ=255)

       gaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAg                                                                                  >  2:8777564/1‑93 (MQ=255)

        aTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGc                                                                                 >  2:1826430/1‑93 (MQ=255)

          gACGGCCAGCAGCCGCGGATGACCCTGAATCACTTCGCCAGCGGTATCTTTGATCCAGTGGATGTTGTTGGAAACGAAAGCATCGCCCAGCCa                                                                               <  1:2024004/93‑1 (MQ=255)

          gACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCa                                                                               <  2:4405095/93‑1 (MQ=255)

            cGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACg                                                                             >  2:6336685/1‑93 (MQ=255)

             ggCCAGCAGCCACGGATGACCCTGAATCACCTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAACATCGCCCACCCCCGc                                                                            >  1:3268015/1‑93 (MQ=255)

              gCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCa                                                                           >  1:5797907/1‑93 (MQ=255)

              gCCAGCAGCCACGGATGACCCTGAAGCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTCGTTAGAAACGAAAGCATCGCCCAGCCACGCa                                                                           <  2:5182211/93‑1 (MQ=255)

                cagcagCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAac                                                                         <  2:1365689/93‑1 (MQ=255)

                cagcagCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAac                                                                         >  1:8610923/1‑93 (MQ=255)

                cagcagCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAac                                                                         >  1:758633/1‑93 (MQ=255)

                cagcagCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAac                                                                         <  1:1688849/93‑1 (MQ=255)

                 agcagcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAaca                                                                        <  2:7816982/93‑1 (MQ=255)

                  gcagcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGGCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAacac                                                                       >  2:2696368/1‑93 (MQ=255)

                    agcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACAccc                                                                     <  1:7875639/93‑1 (MQ=255)

                     gcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCa                                                                    <  2:4704828/93‑1 (MQ=255)

                     gcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCa                                                                    >  1:1365592/1‑93 (MQ=255)

                        aCGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGa                                                                 <  1:9771430/93‑1 (MQ=255)

                         cGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGAt                                                                >  2:2056651/1‑93 (MQ=255)

                         cGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGAt                                                                <  2:2419031/93‑1 (MQ=255)

                          ggATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATa                                                               >  2:6443407/1‑93 (MQ=255)

                          ggATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATa                                                               <  2:4913318/93‑1 (MQ=255)

                            aTGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACa                                                             >  2:797666/1‑93 (MQ=255)

                                cccTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATa                                                         <  1:6573610/93‑1 (MQ=255)

                                cccTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATa                                                         <  2:3424885/93‑1 (MQ=255)

                                 ccTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATAc                                                        <  2:6905657/93‑1 (MQ=255)

                                  cTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACa                                                       <  1:3085285/93‑1 (MQ=255)

                                  cTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACa                                                       >  2:2798441/1‑93 (MQ=255)

                                  cTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACa                                                       <  1:3569722/93‑1 (MQ=255)

                                   ttaaTCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCGCGCAACTCCCAGGATACAAATACAg                                                      <  1:8740018/91‑1 (MQ=255)

                                   tGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAg                                                      >  2:5710791/1‑93 (MQ=255)

                                   tGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAg                                                      >  2:5931983/1‑93 (MQ=255)

                                    gAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAgg                                                     >  2:9273408/1‑93 (MQ=255)

                                     aaTCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGgc                                                    >  1:7765148/1‑93 (MQ=255)

                                       tCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGgcgc                                                  >  1:4463591/1‑93 (MQ=255)

                                       tCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGgcgc                                                  >  2:1508035/1‑93 (MQ=255)

                                          cTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCa                                               >  2:5042216/1‑93 (MQ=255)

                                          cTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCa                                               <  1:5557996/93‑1 (MQ=255)

                                          cTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCa                                               <  2:7425564/93‑1 (MQ=255)

                                           ttCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCAt                                              <  2:8442721/93‑1 (MQ=255)

                                             cACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATAc                                            <  2:7329837/93‑1 (MQ=255)

                                              aCCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATAcc                                           <  1:5905496/93‑1 (MQ=255)

                                               ccAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCt                                          >  2:9702167/1‑93 (MQ=255)

                                                cAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTg                                         <  1:241223/93‑1 (MQ=255)

                                                  gCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCt                                       >  2:2290107/1‑93 (MQ=255)

                                                   cGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCtt                                      >  2:1672588/1‑93 (MQ=255)

                                                    ggTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGCTACAAATACAGGCGCTCATACCTGCttt                                     >  2:6530742/1‑93 (MQ=255)

                                                    ggTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCttt                                     <  2:4515527/93‑1 (MQ=255)

                                                    ggTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCttt                                     >  2:4015642/1‑93 (MQ=255)

                                                    ggTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCttt                                     <  2:3726652/93‑1 (MQ=255)

                                                     gTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTg                                    >  2:1714846/1‑93 (MQ=255)

                                                     gTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTg                                    >  2:1871742/1‑93 (MQ=255)

                                                       aTCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCACCACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGaa                                  <  2:6047370/93‑1 (MQ=255)

                                                       aTCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGaa                                  >  1:1062845/1‑93 (MQ=255)

                                                       aTCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGaa                                  >  2:4775449/1‑93 (MQ=255)

                                                       aTCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGaa                                  >  2:2209726/1‑93 (MQ=255)

                                                        tCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAg                                 <  1:894127/93‑1 (MQ=255)

                                                         cTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCATCACCCAGGATACGAGTACAGGCGCTCATCCCAGCTTTGAAgg                                >  2:5399025/1‑93 (MQ=255)

                                                         cTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAgg                                >  1:3452598/1‑93 (MQ=255)

                                                         cTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAgg                                <  2:529093/93‑1 (MQ=255)

                                                          tttGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAggg                               >  1:5210299/1‑92 (MQ=255)

                                                          tttGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGt                               >  2:4645680/1‑93 (MQ=255)

                                                           ttGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTg                              >  2:9138042/1‑93 (MQ=255)

                                                           ttGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTg                              >  2:1618344/1‑93 (MQ=255)

                                                             gATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCt                            >  1:2896697/1‑93 (MQ=255)

                                                              aTCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCCGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTg                           <  1:9909209/93‑1 (MQ=255)

                                                              aTCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTg                           >  2:8087072/1‑93 (MQ=255)

                                                                 cAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAg                        >  1:5698517/1‑93 (MQ=255)

                                                                  aGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGt                       <  1:6911108/93‑1 (MQ=255)

                                                                  aGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGt                       >  1:2290176/1‑93 (MQ=255)

                                                                  aGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGt                       <  1:1817045/93‑1 (MQ=255)

                                                                    tCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTg                     >  2:9831750/1‑93 (MQ=255)

                                                                     tgATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGa                    <  1:8270708/92‑1 (MQ=255)

                                                                      gATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAg                   <  2:7195363/93‑1 (MQ=255)

                                                                       aTGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAgg                  >  2:8307972/1‑93 (MQ=255)

                                                                       aTGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAgg                  >  1:2627145/1‑93 (MQ=255)

                                                                        tgttgtTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGa                 <  1:9085647/93‑1 (MQ=255)

                                                                         gttgttGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAt                <  1:9533955/93‑1 (MQ=255)

                                                                          ttgttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAtg               >  1:1136293/1‑93 (MQ=255)

                                                                           tgttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAtgt              <  1:8873632/93‑1 (MQ=255)

                                                                           tgttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAtgt              <  1:4110978/93‑1 (MQ=255)

                                                                            gttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAtgtt             >  2:389334/1‑93 (MQ=255)

                                                                            gttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAtgtt             <  1:481921/93‑1 (MQ=255)

                                                                              tgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGAtgttgt           <  2:7541693/93‑1 (MQ=255)

                                                                              gggAAGCGAAAGCATCGCCCGGCCACGCAGCACCCAGGATCAAAAAACAGGCGCTCATCCCTGCTTTGAAGGTGCT‑GAGTTGAGGATGTTGTc          <  1:1479878/92‑1 (MQ=255)

                                                                               ggAAACGAAAGCATCGCCCAGCCACGCAACACCGAGGATACAAATCCAGACGCTCATCCCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTc          <  1:4241091/93‑1 (MQ=255)

                                                                                gAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCg         >  1:9220163/1‑93 (MQ=255)

                                                                                 aaaCGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCgg        <  2:1014300/93‑1 (MQ=255)

                                                                                     gAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATc    <  1:6336686/93‑1 (MQ=255)

                                                                                     gAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATc    <  1:543069/93‑1 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCg   >  1:2539065/1‑93 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCg   >  1:6260698/1‑93 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCg   <  1:459815/93‑1 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCg   >  1:904565/1‑93 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCg   >  2:3070068/1‑93 (MQ=255)

                                                                                       aaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGa  >  1:7345571/1‑93 (MQ=255)

                                                                                       aaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGa  >  2:9511031/1‑93 (MQ=255)

AGAAGAAGATGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGTATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGA  >  NC_000913/4363337‑4363516

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