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%	729.4 / NA	229	T309A (ACT→GCT)	dcuA	anaerobic dicarboxylate transport
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (126/103); total (126/103)

Only 100 of 235 total aligned reads displayed.

AGAAGAAGATGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGTATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGA  >  NC_000913/4363337‑4363516

agaagaagaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCAt                                                                                         >  2:2522879/1‑93 (MQ=255)

agaagaagaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCAt                                                                                         >  2:4096159/1‑93 (MQ=255)

 gaagaagaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATc                                                                                        <  2:3747125/93‑1 (MQ=255)

    gaagaTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAAAGAAAGCATCgtc                                                                                     >  2:1872087/1‑91 (MQ=255)

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

        aTGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGc                                                                                 <  1:3230127/93‑1 (MQ=255)

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

           aCGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCAc                                                                              >  2:4321753/1‑93 (MQ=255)

             ggCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGc                                                                            <  1:3810608/93‑1 (MQ=255)

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

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

                cagcagCCACGGATGACCCTGAATCTCTTCACCAGCGGTATCTTTGATCCAGTCGAGGTTGTTTGAAACGAAAGCATCGCCCAGCCACGCAac                                                                         <  2:929836/93‑1 (MQ=255)

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

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

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

                  gcagcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAacac                                                                       >  2:2341112/1‑93 (MQ=255)

                  gcagcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAacac                                                                       >  2:2774922/1‑93 (MQ=255)

                  gcagcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAacac                                                                       <  1:4184200/93‑1 (MQ=255)

                   cagcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACAcc                                                                      <  2:624677/93‑1 (MQ=255)

                    agcCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACAccc                                                                     >  1:2853609/1‑93 (MQ=255)

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

                    agcCACGGATGACCCTGAAACACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACAccc                                                                     <  1:2899092/93‑1 (MQ=255)

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

                      ccACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAg                                                                   <  1:1178025/93‑1 (MQ=255)

                      ccACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGACGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAg                                                                   <  1:3304077/93‑1 (MQ=255)

                        aCGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGa                                                                 >  2:4345740/1‑93 (MQ=255)

                        aCGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGa                                                                 >  1:2375562/1‑93 (MQ=255)

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

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

                           gATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATAc                                                              >  2:4623886/1‑93 (MQ=255)

                           gATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATAc                                                              <  1:3344103/93‑1 (MQ=255)

                           gATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATAc                                                              >  2:3394292/1‑93 (MQ=255)

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

                              gACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACaaa                                                           <  1:764506/93‑1 (MQ=255)

                              gACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACaaa                                                           <  1:1920002/93‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                      aTCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGgcg                                                   >  1:853563/1‑93 (MQ=255)

                                           ttCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCAt                                              >  1:4295351/1‑93 (MQ=255)

                                           ttCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCAt                                              >  1:3909463/1‑93 (MQ=255)

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

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

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

                                                 aGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGc                                        <  1:3992823/93‑1 (MQ=255)

                                                  gCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCt                                       <  1:2483093/93‑1 (MQ=255)

                                                   cGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCtt                                      <  1:3678684/93‑1 (MQ=255)

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

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

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

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

                                                      tATCTTTGATCAAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAAGCGCTCATACCTGCTTTGa                                   >  2:2242287/1‑93 (MQ=255)

                                                       aTCTTTGTTCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGaa                                  >  1:711618/1‑93 (MQ=255)

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

                                                         cTTTGATCCAGTGGATGTTGTTGAAAACGAAAGCATCGCCCAGCCACGCAACCCCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAgg                                <  2:3702557/93‑1 (MQ=255)

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

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

                                                          tttGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGt                               <  2:1260987/93‑1 (MQ=255)

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

                                                          tttGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACACAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGt                               >  2:298538/1‑93 (MQ=255)

                                                            tGATCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGc                             >  1:1868773/1‑93 (MQ=255)

                                                               tCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTgg                          >  1:1263356/1‑93 (MQ=255)

                                                               tCCAGTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTgg                          >  2:3217572/1‑93 (MQ=255)

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

                                                                  cgTCGATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGt                       <  2:2921108/92‑1 (MQ=255)

                                                                   gTNGATGTTGTTGGAAACGAAANCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGtt                      <  1:3438939/93‑1 (MQ=255)

                                                                      gATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAg                   >  2:1721325/1‑93 (MQ=255)

                                                                      gATGTTGTTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAg                   >  1:3757236/1‑93 (MQ=255)

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

                                                                        tgttgtTGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAggg                 >  2:3533768/1‑92 (MQ=255)

                                                                         gttgttGGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAgggt                <  1:1049680/93‑3 (MQ=255)

                                                                           tgttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGgtgt              <  2:3017598/93‑1 (MQ=255)

                                                                             ttgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGtgttg            >  1:2553036/1‑93 (MQ=255)

                                                                              tgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGtgttgt           >  1:3223746/1‑93 (MQ=255)

                                                                              tgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGtgttgt           <  1:249665/93‑1 (MQ=255)

                                                                              tgGAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGtgttgt           <  1:2647149/93‑1 (MQ=255)

                                                                               ggAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGtgttgtg          >  1:3515597/1‑92 (MQ=255)

                                                                               ggAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTc          >  2:117402/1‑93 (MQ=255)

                                                                               ggAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTc          <  1:1105514/93‑1 (MQ=255)

                                                                                gAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCg         >  1:3151833/1‑93 (MQ=255)

                                                                                gAAACGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCg         >  2:3476980/1‑93 (MQ=255)

                                                                                 aaaCGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCgg        <  2:4287385/93‑1 (MQ=255)

                                                                                 aaaCGAAAGCAGCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCgg        <  2:2310901/93‑1 (MQ=255)

                                                                                  aaCGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGt       >  2:793150/1‑93 (MQ=255)

                                                                                  aaCGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGt       >  2:1622383/1‑93 (MQ=255)

                                                                                    cGAAAGCATCGCCCTGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTAt     >  2:2864705/1‑93 (MQ=255)

                                                                                    cGAAAGCATCGCCCAGCGACGCAGCACCCAGGATACAAATACAGGCGCTCATGCCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTAt     <  1:140439/93‑1 (MQ=255)

                                                                                    cGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTAt     >  1:2427029/1‑93 (MQ=255)

                                                                                    cGAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGCAGTTGAGGGTGTTGTCGGTAt     <  2:2982829/93‑1 (MQ=255)

                                                                                     gAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTATc    <  2:4028969/93‑1 (MQ=255)

                                                                                     gAAAGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTATc    >  2:1576991/1‑93 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCTGTATCg   >  2:1043716/1‑93 (MQ=255)

                                                                                      aaaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTATCg   >  2:2900453/1‑93 (MQ=255)

                                                                                       aaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTATCGa  <  1:2111305/93‑1 (MQ=255)

                                                                                       aaGCATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGGTGTTGTCGGTATCGa  <  1:3576785/93‑1 (MQ=255)

AGAAGAAGATGACGGCCAGCAGCCACGGATGACCCTGAATCACTTCACCAGCGGTATCTTTGATCCAGTCGATGTTGTTGGAAACGAAAGTATCGCCCAGCCACGCAACACCCAGGATACAAATACAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGA  >  NC_000913/4363337‑4363516

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