BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	258,977	(C)_5→6	coding (348/1056 nt)	phoE ←	outer membrane pore protein E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	258,972	1	.	C	100.0%	492.2 / NA	158	A118G (GCG→GGG)	phoE	outer membrane pore protein E
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (81/77); total (81/77)

Only 100 of 181 total aligned reads displayed.

TGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑CCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGTAT  >  NC_000913/258880‑259059

tGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACg                                                                                           <  2:1081189/93‑1 (MQ=255)

 ggCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑c                                                                                         >  2:1527451/1‑93 (MQ=255)

 ggCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑c                                                                                         >  2:3178242/1‑93 (MQ=255)

  gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cc                                                                                        <  2:2703134/93‑1 (MQ=255)

  gcgcGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cc                                                                                        >  1:353936/1‑93 (MQ=255)

     cgTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑c                                                                                    <  2:4430849/93‑1 (MQ=255)

      gTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑cc                                                                                   >  2:3976081/1‑92 (MQ=255)

      gTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑cc                                                                                   >  1:1739191/1‑92 (MQ=255)

       tttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑cca                                                                                  <  1:2003872/93‑3 (MQ=255)

        ttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑ccag                                                                                 <  2:1483305/93‑4 (MQ=255)

        ttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑ccag                                                                                 <  1:2046402/93‑4 (MQ=255)

        ttGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑ccag                                                                                 <  1:1662377/93‑4 (MQ=255)

         tGGTCATAAAGTTGTCGGTGTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑cccc‑ccagg                                                                                <  2:933379/93‑5 (MQ=255)

           gTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGtt                                                                               >  1:1107022/1‑93 (MQ=255)

           gTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGtt                                                                               >  1:571187/1‑93 (MQ=255)

             taTAAAGTTGTCGGTCTGCGCCGAGGAAGCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTAc                                                                             <  2:2604187/92‑1 (MQ=255)

             cATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTAc                                                                             >  1:2694113/1‑93 (MQ=255)

               tAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGa                                                                           >  1:492411/1‑93 (MQ=255)

                aaaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGAc                                                                          <  1:4174272/93‑1 (MQ=255)

                aaaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGAc                                                                          <  2:465683/93‑1 (MQ=255)

                aaaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGGTACGAc                                                                          >  2:310650/1‑93 (MQ=255)

                 aaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGAcc                                                                         >  1:4016853/1‑93 (MQ=255)

                 aaGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGAcc                                                                         >  1:2773244/1‑93 (MQ=255)

                  aGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCa                                                                        <  2:3279568/93‑1 (MQ=255)

                   gTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCAt                                                                       >  1:4468300/1‑93 (MQ=255)

                    ttGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATa                                                                      >  2:3570163/1‑93 (MQ=255)

                      gTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAAt                                                                    <  2:121421/93‑1 (MQ=255)

                      gTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAAt                                                                    <  1:4609040/93‑1 (MQ=255)

                      gTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAAt                                                                    <  1:2968444/93‑1 (MQ=255)

                       tCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATc                                                                   <  1:719805/93‑1 (MQ=255)

                        cGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCg                                                                  <  2:3618016/93‑1 (MQ=255)

                         ggTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGa                                                                 <  2:2612760/93‑1 (MQ=255)

                          gTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGaa                                                                >  2:1724091/1‑93 (MQ=255)

                           tCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGaaa                                                               >  1:4536621/1‑93 (MQ=255)

                           tCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGaaa                                                               <  1:1192924/93‑1 (MQ=255)

                           tCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGaaa                                                               >  1:4320549/1‑93 (MQ=255)

                             tGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGa                                                             <  1:2383731/93‑1 (MQ=255)

                              gcgcCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGaa                                                            <  2:2793905/93‑1 (MQ=255)

                               cgcCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCATGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAAc                                                           <  1:3301919/93‑1 (MQ=255)

                               cgcCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAAc                                                           <  2:3115648/93‑1 (MQ=255)

                                  cGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCa                                                        <  2:2428717/93‑1 (MQ=255)

                                   gAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCaa                                                       >  2:3798850/1‑93 (MQ=255)

                                    aGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCaaa                                                      <  2:3138928/93‑1 (MQ=255)

                                    aGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCaaa                                                      <  2:553188/93‑1 (MQ=255)

                                    aGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCaaa                                                      <  2:367652/93‑1 (MQ=255)

                                     ggAATCGCCGCCAAATTCCGGTACCATATCGGTCCAGTCTTCCACGTCTTACAACC‑CCCCGCAGGTTACGACCATAATCGAAAGAACCCaaag                                                     <  2:1176144/93‑2 (MQ=255)

                                        aTCGCCACCAAATTGCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCtt                                                  <  1:2714291/93‑1 (MQ=255)

                                        aTCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTGTGACCATAATCGAAAGAACCCAAATCtt                                                  <  2:3219528/93‑1 (MQ=37)

                                        aTCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCtt                                                  >  2:2145614/1‑93 (MQ=255)

                                           tccaccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTtat                                               <  2:719762/92‑1 (MQ=255)

                                           gCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTtat                                               <  2:4186674/93‑1 (MQ=255)

                                             caccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTatat                                             >  2:1806798/1‑93 (MQ=255)

                                             caccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTatat                                             >  1:1869803/1‑93 (MQ=255)

                                              cccaAATTCCGCGAACATATCGGTCCAGGCTTCCACGTCATCCAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATAtt                                            <  2:1278433/92‑1 (MQ=255)

                                              accaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATAtt                                            >  2:2242962/1‑93 (MQ=255)

                                              accaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATAtt                                            <  1:1058762/93‑1 (MQ=255)

                                              accaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATAtt                                            >  1:3949686/1‑93 (MQ=255)

                                              accaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATAtt                                            >  1:1271576/1‑93 (MQ=255)

                                               ccaAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATAttt                                           <  2:3668804/93‑1 (MQ=255)

                                                caAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTc                                          >  1:1009934/1‑93 (MQ=255)

                                                caAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTc                                          <  1:720835/93‑1 (MQ=255)

                                                caAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTc                                          >  2:4159821/1‑93 (MQ=255)

                                                caAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATAAAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTc                                          >  2:1128783/1‑93 (MQ=255)

                                                  aaTTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCaa                                        >  2:1124595/1‑93 (MQ=255)

                                                   aTTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAAc                                       <  1:901784/93‑1 (MQ=255)

                                                    ttCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAAcc                                      <  1:2188542/93‑1 (MQ=255)

                                                    ttCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAAcc                                      <  1:2923252/93‑1 (MQ=255)

                                                      ccGGGAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCg                                    >  2:230433/1‑93 (MQ=255)

                                                        gggAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGc                                  >  2:3689752/1‑93 (MQ=255)

                                                        gggAACATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGc                                  >  1:8345/1‑93 (MQ=255)

                                                            aCATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCaaaa                              <  1:3751881/93‑1 (MQ=255)

                                                            aCATATCGGTCCAGGCTGCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCaaaa                              <  1:4388638/93‑1 (MQ=255)

                                                             cATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAg                             <  1:246650/93‑1 (MQ=255)

                                                             cATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAg                             >  2:2964972/1‑93 (MQ=255)

                                                             cATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAg                             <  1:3831217/93‑1 (MQ=255)

                                                             cATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAg                             <  1:3173894/93‑1 (MQ=255)

                                                             cATATCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAg                             <  1:1975310/93‑1 (MQ=255)

                                                              atatCGGTCCAGGCTTCCACGTCATCCAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAGATCTTTATATTTCAACCCCGCAAAAGc                            >  2:1396191/1‑93 (MQ=255)

                                                              atatCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGc                            <  2:691942/93‑1 (MQ=255)

                                                               tatCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCg                           <  2:792360/93‑1 (MQ=255)

                                                               tatCGGTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCg                           >  1:1487817/1‑93 (MQ=255)

                                                                   ggTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGAc                       >  1:3939340/1‑93 (MQ=255)

                                                                   ggTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGAc                       <  1:4268499/93‑1 (MQ=255)

                                                                    gTCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGAcg                      >  1:2555081/1‑93 (MQ=255)

                                                                     tCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGAcgc                     >  2:3345629/1‑93 (MQ=255)

                                                                     tCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGAcgc                     >  1:1205789/1‑93 (MQ=255)

                                                                     tCCAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGAcgc                     >  2:2564961/1‑93 (MQ=255)

                                                                       cAGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGt                   >  1:1243692/1‑93 (MQ=255)

                                                                        aGGCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGtt                  >  1:1656797/1‑93 (MQ=255)

                                                                          gCTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCAT‑ATCGGAAAAACCCAAATCTTTATATTTCAAACCGGCAAAAGCGAGACGAGgtttt               >  2:4067407/1‑93 (MQ=255)

                                                                           cTTCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGttttt               <  2:2881361/93‑1 (MQ=255)

                                                                             tCCACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTg             <  1:3415083/93‑1 (MQ=255)

                                                                              ccACGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGc            <  2:3416122/93‑1 (MQ=255)

                                                                                aCGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCtg          >  1:101257/1‑93 (MQ=255)

                                                                                aCGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCtg          >  1:2975569/1‑93 (MQ=255)

                                                                                aCGTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCtg          >  1:2006161/1‑93 (MQ=255)

                                                                                  gTCATACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCtgtg        >  1:3544478/1‑93 (MQ=255)

                                                                                     aTACAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAg     >  2:553498/1‑93 (MQ=255)

                                                                                      aaCAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGt    <  2:2055623/92‑1 (MQ=255)

                                                                                        cAACGCCCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGTAt  >  1:1580678/1‑93 (MQ=255)

TGGCGCGTTTGGTCATAAAGTTGTCGGTCTGCGCCGAGGAATCGCCACCAAATTCCGGGAACATATCGGTCCAGGCTTCCACGTCATACAACG‑CCCC‑CAGGTTACGACCATAATCGAAAGAACCCAAATCTTTATATTTCAACCCGGCAAAAGCGAGACGCGTTTTTTGCTGTGCAGTAT  >  NC_000913/258880‑259059

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