BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,555	(G)_5→6	coding (55/618 nt)	b3646 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,550	1	.	G	100.0%	716.5 / NA	226	T17R (ACG→AGG)	b3646	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base G (118/108); total (118/108)

Only 100 of 250 total aligned reads displayed.

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

 acaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGATCCATGCTGTTACACATTTTGTATTTGGTTGGCATTAGTGCCGAAGCCATGAC‑g                                                                                        >  2:821546/1‑93 (MQ=255)

 acaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑g                                                                                        <  1:313347/93‑1 (MQ=255)

 acaAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑g                                                                                        >  1:4675744/1‑93 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                      <  1:329823/93‑1 (MQ=255)

   aaaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggg                                                                                      <  2:698965/93‑1 (MQ=255)

    aaTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑gggg                                                                                     >  1:3168101/1‑93 (MQ=255)

     aTCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggg                                                                                    >  1:2081922/1‑93 (MQ=255)

        gCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggggcg                                                                                 >  1:4625747/1‑90 (MQ=255)

         cTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑ggggggcgc                                                                                >  1:505121/1‑89 (MQ=255)

                gCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTggcggc                                                                          >  2:349566/1‑93 (MQ=255)

                    gacagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGga                                                                      >  2:3375713/1‑93 (MQ=255)

                      cagaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacg                                                                    >  1:2421300/1‑93 (MQ=255)

                       agaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacga                                                                   <  1:1496407/93‑1 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGTCCGgacgac                                                                  >  1:741624/1‑93 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  <  2:1813393/93‑1 (MQ=255)

                        gaTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  <  2:1173948/93‑1 (MQ=255)

                        gaTTTTTAGGAAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGgacgac                                                                  <  2:2688872/93‑1 (MQ=255)

                         aTTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGacgacg                                                                 <  2:2828697/93‑1 (MQ=255)

                            tttAGGGAGAGAACCATGCTGTTACACATTTTGCATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGTCGGACGACGgcg                                                              >  2:116445/1‑93 (MQ=255)

                             ttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgcgc                                                             <  2:2352711/93‑1 (MQ=255)

                             ttAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGgcgc                                                             >  1:690823/1‑93 (MQ=255)

                              tAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCa                                                            <  1:4169454/93‑1 (MQ=255)

                               aGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCAt                                                           >  2:2999213/1‑93 (MQ=255)

                                gggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATg                                                          <  2:1289147/93‑1 (MQ=255)

                                gggAGAGAACCATGCTGTTACACATTTTGCATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCCTg                                                          >  2:1106680/1‑93 (MQ=255)

                                 ggAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATgg                                                         >  2:570905/1‑93 (MQ=255)

                                  gagagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGa                                                        >  2:4487793/1‑93 (MQ=255)

                                   agagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGAt                                                       >  1:4491009/1‑93 (MQ=255)

                                    gagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGGATGGaaa                                                      >  2:4019458/1‑91 (MQ=255)

                                    gagaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATa                                                      <  2:1316632/93‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGTATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    <  2:3203707/93‑1 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    >  2:2625233/1‑93 (MQ=255)

                                      gaACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACa                                                    <  1:4353663/93‑1 (MQ=255)

                                        aCCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGAAGACGGCGAATGGATACAtt                                                  >  2:1880160/1‑93 (MQ=255)

                                         ccATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACAttt                                                 <  2:3111755/93‑1 (MQ=255)

                                          cATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTg                                                >  2:378122/1‑93 (MQ=255)

                                          cATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTg                                                <  2:1091270/93‑1 (MQ=255)

                                          cATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTg                                                >  2:750218/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               <  1:970075/93‑1 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               <  2:2437447/93‑1 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               >  2:4439479/1‑93 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgg                                               <  2:4631510/93‑1 (MQ=255)

                                           aTGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTgc                                               >  1:4441446/1‑92 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  2:3241000/1‑93 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  1:1429462/1‑93 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  1:3628365/1‑93 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           <  1:1737971/93‑1 (MQ=255)

                                               tGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTa                                           >  2:1436792/1‑93 (MQ=255)

                                                gTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTaa                                          >  2:3186707/1‑93 (MQ=255)

                                                gTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTaa                                          >  1:2279892/1‑93 (MQ=255)

                                                gTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTaa                                          >  1:2502133/1‑93 (MQ=255)

                                                  tACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatt                                        >  1:4049025/1‑93 (MQ=255)

                                                  tACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAatt                                        <  1:2157196/93‑1 (MQ=255)

                                                     acaTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatt                                     >  2:1611102/1‑93 (MQ=255)

                                                      caTTTTTTATTTGGTGGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  2:4435567/93‑1 (MQ=255)

                                                      caTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    >  2:469500/1‑93 (MQ=255)

                                                      aaTTTTGTATTTGGTTGGCATGACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattatta                                    <  2:228068/92‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   <  1:4570606/93‑1 (MQ=255)

                                                       aTTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattat                                   >  1:1328023/1‑93 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  2:1794866/93‑1 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  1:1979852/93‑1 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  >  2:4113140/1‑93 (MQ=255)

                                                        ttttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAattattatt                                  <  2:3736602/93‑1 (MQ=255)

                                                         tttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTg                                 <  1:2974259/93‑1 (MQ=255)

                                                          ttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGc                                <  2:879740/93‑1 (MQ=255)

                                                          ttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGc                                >  2:2101636/1‑93 (MQ=255)

                                                          ttGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGc                                >  1:2773521/1‑93 (MQ=255)

                                                           tGTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCg                               >  1:968232/1‑93 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  2:3029378/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  1:3225733/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  1:3528170/93‑1 (MQ=255)

                                                             tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             >  2:3582202/1‑93 (MQ=255)

                                                             tATTTGGTTGGCATGACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGAc                             <  2:4155609/93‑1 (MQ=255)

                                                               tttGGTTGGCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCg                           >  1:2872573/1‑93 (MQ=255)

                                                                  ggttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCa                        >  1:1810955/1‑93 (MQ=255)

                                                                   gttggCATTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGTCAc                       >  1:390970/1‑93 (MQ=255)

                                                                   gttggCATTACTGCCGAAGCCAGGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCAc                       <  2:3961678/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:1141452/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:2586331/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:368454/93‑1 (MQ=255)

                                                                         aTTACTGCCGAAGCAATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAt                 <  1:1014969/93‑1 (MQ=255)

                                                                          ttACTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAAtt                <  2:4330826/93‑1 (MQ=255)

                                                                            aCTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTgg              <  2:4256472/93‑1 (MQ=255)

                                                                             cTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGc             <  1:638443/93‑1 (MQ=255)

                                                                             cTGCCGAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGc             >  2:3391365/1‑93 (MQ=255)

                                                                                  gAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        <  2:505121/93‑1 (MQ=255)

                                                                                  gAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        >  1:1925394/1‑93 (MQ=255)

                                                                                  gAAGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgg        >  1:4004644/1‑93 (MQ=255)

                                                                                   aaGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAggg       <  1:456139/93‑1 (MQ=255)

                                                                                   aaGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAggg       >  1:4021276/1‑93 (MQ=255)

                                                                                   aaGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAggg       >  2:2061991/1‑93 (MQ=255)

                                                                                    aGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAgggg      >  1:1245815/1‑92 (MQ=255)

                                                                                    aGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      <  2:2793717/93‑1 (MQ=255)

                                                                                    aGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      >  2:453421/1‑93 (MQ=255)

                                                                                    aGCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGt      >  1:3969614/1‑93 (MQ=255)

                                                                                     gCCATGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTc     <  2:2486331/93‑1 (MQ=255)

                                                                                        aTGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  1:3349703/93‑1 (MQ=255)

                                                                                        aTGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  1:4224635/93‑1 (MQ=255)

                                                                                        aTGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  2:1902758/93‑1 (MQ=255)

                                                                                        aTGACGGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGt  <  1:2633901/93‑1 (MQ=255)

CACAAATCGCTACCCTGCCAGACAGATTTTTAGGGAGAGAACCATGCTGTTACACATTTTGTATTTGGTTGGCATTACTGCCGAAGCCATGAC‑GGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGT  >  NC_000913/3816458‑3816637

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