BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,998,290	+G	coding (379/486 nt)	b4566 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,998,289	1	.	G	100.0%	699.0 / NA	193	L127P (CTG→CCG)	b4566	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (109/84); total (109/84)

Only 100 of 216 total aligned reads displayed.

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑G‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                           <  1:5635226/93‑1 (MQ=255)

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                           <  1:1328206/93‑1 (MQ=255)

gACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCa                                                                                           <  1:1712167/93‑1 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑g                                                                                         >  2:5499598/1‑93 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑g                                                                                         >  2:2644010/1‑93 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑g                                                                                         <  2:5258777/93‑1 (MQ=255)

 aCTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCAA‑g                                                                                         >  2:1973097/1‑93 (MQ=255)

  cTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑gg                                                                                        >  1:3971482/1‑92 (MQ=255)

   cgCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑gga                                                                                       <  2:4990243/92‑3 (MQ=255)

    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑ggac                                                                                      >  1:696906/1‑90 (MQ=255)

    gCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑ggac                                                                                      <  2:1972155/93‑4 (MQ=255)

     cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑ggacg                                                                                     >  2:4511728/1‑89 (MQ=255)

     cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑ggacg                                                                                     <  2:4440440/93‑5 (MQ=255)

     cACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑ggacg                                                                                     <  2:5796335/93‑5 (MQ=255)

       ccGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGa                                                                                   >  1:1571054/1‑93 (MQ=255)

        cGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGAt                                                                                  >  1:591351/1‑93 (MQ=255)

         gCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATa                                                                                 <  2:3562809/93‑1 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAg                                                                              >  2:4568982/1‑93 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAg                                                                              >  1:1789415/1‑93 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAg                                                                              >  2:5528021/1‑93 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAg                                                                              <  1:860584/93‑1 (MQ=255)

            aaCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAg                                                                              >  2:5825273/1‑93 (MQ=255)

             aCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGa                                                                             <  2:3285954/93‑1 (MQ=255)

             aCTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCCACATAAAACACGCCGCCAGG‑ACGGATAAAGa                                                                             >  2:5072575/1‑93 (MQ=255)

              cTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGaa                                                                            >  2:4155906/1‑93 (MQ=255)

               tGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGaaa                                                                           <  1:1154236/93‑1 (MQ=255)

                gAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAAt                                                                          >  2:1245867/1‑93 (MQ=255)

                  aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTc                                                                        >  2:654569/1‑93 (MQ=255)

                   tctcTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTcc                                                                       <  1:293892/93‑1 (MQ=255)

                    ctctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTccc                                                                      <  1:5090373/93‑1 (MQ=255)

                     tctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCt                                                                     >  2:2952643/1‑93 (MQ=255)

                      ctTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCtt                                                                    <  1:5832206/93‑1 (MQ=255)

                       tttAGTGGTTTTGAAGGCATCCAGACGAACCATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCGA‑GCACGGATAACGAAGTTCCCttt                                                                   >  1:3739607/1‑92 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTc                                                                   >  2:1360929/1‑93 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTc                                                                   >  2:3078234/1‑93 (MQ=255)

                       tttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTc                                                                   <  1:3180790/93‑1 (MQ=255)

                        ttACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCa                                                                  <  1:5765551/93‑1 (MQ=255)

                          aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCAcc                                                                >  2:3023099/1‑93 (MQ=255)

                          aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCAcc                                                                <  1:3249322/93‑1 (MQ=255)

                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCt                                                               <  2:3739674/93‑1 (MQ=255)

                           cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCt                                                               >  1:430062/1‑93 (MQ=255)

                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTg                                                              >  2:3986301/1‑93 (MQ=255)

                            tGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTg                                                              >  2:5594334/1‑93 (MQ=255)

                              gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGaa                                                            <  2:2416226/93‑1 (MQ=255)

                               ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAg                                                           >  1:1835307/1‑93 (MQ=255)

                                   gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCgag                                                       >  2:3735396/1‑93 (MQ=255)

                                      ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAtt                                                    >  1:1004224/1‑93 (MQ=255)

                                       gCATCCAGCCGAGCGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCCCCTGAAGCGAGAttt                                                   <  1:3297478/93‑1 (MQ=255)

                                        cATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTa                                                  <  1:942630/93‑1 (MQ=255)

                                          tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAc                                                >  1:2753039/1‑93 (MQ=255)

                                           ccAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAAcc                                               >  2:5020598/1‑92 (MQ=255)

                                            cAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGc                                              <  2:184354/93‑1 (MQ=255)

                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                             >  2:1667708/1‑93 (MQ=255)

                                             aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                                             >  2:3675634/1‑93 (MQ=255)

                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  1:4671524/1‑93 (MQ=255)

                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  1:6131783/1‑93 (MQ=255)

                                              gACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                                            >  1:1256320/1‑93 (MQ=255)

                                               acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAc                                           <  1:5793530/93‑1 (MQ=255)

                                                 gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                                         >  2:5698776/1‑93 (MQ=255)

                                                  aacgaTGCCGCCTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        <  2:4332744/93‑1 (MQ=255)

                                                  aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                                        >  2:2463816/1‑93 (MQ=255)

                                                    cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                      >  1:1594897/1‑93 (MQ=255)

                                                    cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                      <  1:2872189/93‑1 (MQ=255)

                                                    cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                      <  2:1552288/93‑1 (MQ=255)

                                                    cgaTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTct                                      <  2:2019583/93‑1 (MQ=255)

                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    >  1:2842281/1‑93 (MQ=255)

                                                      aTGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgc                                    >  2:2471379/1‑93 (MQ=255)

                                                          cGACTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc                                >  2:1388244/1‑93 (MQ=255)

                                                            aCTTTATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCg                              >  1:4749325/1‑93 (MQ=255)

                                                              tttATCGCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtt                            <  2:2134740/93‑1 (MQ=255)

                                                               ttATCTCCCGCGACATACAGCACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           <  2:5453192/93‑1 (MQ=255)

                                                               ttATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt                           >  1:3423079/1‑93 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:5370405/93‑1 (MQ=255)

                                                                tATCTCCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt                          <  2:256956/93‑1 (MQ=255)

                                                                 aTCTCCCGCAACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc                         <  1:2362459/93‑1 (MQ=255)

                                                                  tctcCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTca                        >  2:4673864/1‑92 (MQ=255)

                                                                  tctcCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTCTCg                        >  2:3174009/1‑93 (MQ=255)

                                                                   ctcCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGa                       <  2:4857955/93‑1 (MQ=255)

                                                                   ctcCCGCGACATAAAACACGCCGCCAGG‑ACGGAGAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGa                       <  2:5931943/93‑1 (MQ=255)

                                                                    tcCCGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt                      <  1:3252680/93‑1 (MQ=255)

                                                                      ccGCGCCATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    <  1:1357433/93‑1 (MQ=255)

                                                                      ccGCGACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt                    <  2:4737377/93‑1 (MQ=255)

                                                                        gcgACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTAc                  <  2:4682551/93‑1 (MQ=255)

                                                                         cgACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACa                 <  1:5505831/93‑1 (MQ=255)

                                                                          gACATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaa                >  1:4349542/1‑93 (MQ=255)

                                                                           aCATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               <  2:2898757/93‑1 (MQ=255)

                                                                           aCATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               >  1:4160831/1‑93 (MQ=255)

                                                                           aCATAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACaaa               <  1:3726579/93‑1 (MQ=255)

                                                                            cATAAAACAAGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAAc              >  2:4972595/1‑93 (MQ=255)

                                                                              tAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCg            >  1:2538520/1‑93 (MQ=255)

                                                                              tAAAACACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCCGACCGTTTTCGATGTACAAACCg            >  1:3502798/1‑93 (MQ=255)

                                                                               aaaaCACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGa           <  1:40147/93‑1 (MQ=255)

                                                                                 aaCACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgc         >  2:47850/1‑93 (MQ=255)

                                                                                 aaCACGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAACTTCTGCTGACCGTGTTCGAAGTACAAACCGAgc         >  2:6101965/1‑93 (MQ=255)

                                                                                  acacGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGAACAAACCGAgcg        >  2:2505681/1‑93 (MQ=255)

                                                                                  acacGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGATTTCGATGTACAAACCGAgcg        <  1:6038418/93‑1 (MQ=255)

                                                                                    acGCCGCCAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAgcggc      <  1:1645394/93‑1 (MQ=255)

                                                                                       ccgccAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCa   >  2:1147602/1‑93 (MQ=255)

                                                                                        cgccAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTTCAAACCGAGCGGCGCAt  <  1:3638168/93‑1 (MQ=255)

                                                                                        cgccAGG‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAt  >  1:242206/1‑93 (MQ=255)

GACTGCACCGCAAACTGAATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGCCGCCA‑G‑ACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGTACAAACCGAGCGGCGCAT  >  NC_000913/3998197‑3998376

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