BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,606,004	T→C	G318G (GGT→GGC)	lsrB →	Autoinducer‑2 ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,606,004	0	T	C	100.0%	1196.2 / NA	289	G318G (GGT→GGC)	lsrB	Autoinducer‑2 ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (143/146); total (143/146)

Only 100 of 297 total aligned reads displayed.

GGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTAT  >  NC_000913/1605906‑1606101

ggATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCaa                                                                                                 <  1:2698302/101‑1 (MQ=255)

 gATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAAt                                                                                                >  1:2698817/1‑101 (MQ=255)

   tCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATgg                                                                                              <  2:2698307/101‑1 (MQ=255)

    cAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGc                                                                                             >  2:2698822/1‑101 (MQ=255)

     aaTGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGatggcaatggca                                                                                            >  2:2698826/1‑101 (MQ=255)

      aTGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCAt                                                                                           >  2:2698830/1‑101 (MQ=255)

      aTGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCAt                                                                                           <  2:2698344/101‑1 (MQ=255)

       tGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATc                                                                                          >  1:2698833/1‑101 (MQ=255)

         aaaaCGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGt                                                                                        >  2:2698835/1‑101 (MQ=255)

           aaCGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTAc                                                                                      >  2:2698840/1‑101 (MQ=255)

            aCGGGCGACAAGCTGGACATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGGAGCGGATGGCAATGGCATCGTACt                                                                                     >  2:2698843/1‑101 (MQ=255)

             cGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTg                                                                                    >  1:2698844/1‑101 (MQ=255)

             cGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTg                                                                                    <  2:2698300/101‑1 (MQ=255)

                  gACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTAcc                                                                               >  1:2698849/1‑101 (MQ=255)

                   aCAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCg                                                                              >  2:2698850/1‑101 (MQ=255)

                    cAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCgg                                                                             <  1:2698201/101‑1 (MQ=255)

                    cAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCgg                                                                             <  1:2698060/101‑1 (MQ=255)

                     aaGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGa                                                                            >  1:2698854/1‑101 (MQ=255)

                     aaGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGa                                                                            >  2:2698855/1‑101 (MQ=255)

                        cTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAgcg                                                                         <  2:2698397/101‑1 (MQ=255)

                         tGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAgcgc                                                                        >  2:2698859/1‑101 (MQ=255)

                          ggATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGcgcg                                                                       >  2:2698864/1‑101 (MQ=255)

                          ggATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGcgcg                                                                       >  1:2698863/1‑101 (MQ=255)

                           gATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAg                                                                           >  2:2698869/1‑96 (MQ=255)

                            atatCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTg                                                                     <  2:2698265/101‑1 (MQ=255)

                            atatCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTg                                                                     >  1:2698872/1‑101 (MQ=255)

                            atatCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTg                                                                     >  1:2698870/1‑101 (MQ=255)

                             tatCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGa                                                                    >  2:2698877/1‑101 (MQ=255)

                             tatCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGa                                                                    >  1:2698876/1‑101 (MQ=255)

                              atCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGat                                                                   <  1:2698484/101‑1 (MQ=255)

                               tCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAgc                                                                          >  1:2698882/1‑93 (MQ=255)

                                cAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGatat                                                                 >  2:2698883/1‑101 (MQ=255)

                                cAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGatat                                                                 <  1:2698400/101‑1 (MQ=255)

                                 aaGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATAtt                                                                >  2:2698887/1‑101 (MQ=255)

                                 aaGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATAtt                                                                <  2:2698364/101‑1 (MQ=255)

                                  aGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTc                                                               >  1:2698889/1‑101 (MQ=255)

                                    ggCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTcaa                                                             <  2:2698527/101‑1 (MQ=255)

                                     gCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTcaac                                                            <  1:2698327/101‑1 (MQ=255)

                                      cGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTcaaca                                                           <  2:2698330/101‑1 (MQ=255)

                                       gTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTcaacaa                                                          >  2:2698894/1‑101 (MQ=255)

                                        tAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACaaa                                                         <  1:2698781/101‑1 (MQ=255)

                                        tAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACaaa                                                         >  2:2698897/1‑101 (MQ=255)

                                         aGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAag                                                        >  1:2698900/1‑101 (MQ=255)

                                         aGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAag                                                        <  1:2698121/101‑1 (MQ=255)

                                         aGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAag                                                        >  2:2698901/1‑101 (MQ=255)

                                             cAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGaa                                                    >  2:2698906/1‑101 (MQ=255)

                                              aGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCg                                                                              <  2:2698506/74‑1 (MQ=255)

                                              aGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAat                                                   >  1:2698908/1‑101 (MQ=255)

                                               ggTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAata                                                  <  1:2698458/101‑1 (MQ=255)

                                               ggTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAata                                                  <  2:2698239/101‑1 (MQ=255)

                                                gTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAatat                                                 <  1:2698283/101‑1 (MQ=255)

                                                 ttGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAAt                                                                                                >  1:2698912/1‑53 (MQ=255)

                                                 ttGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATc                                                <  1:2698425/101‑1 (MQ=255)

                                                 ttGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATc                                                <  2:2698408/101‑1 (MQ=255)

                                                  tGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCg                                               <  1:2698382/101‑1 (MQ=255)

                                                  tGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCg                                               >  2:2698913/1‑101 (MQ=255)

                                                   gAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCgg                                              <  1:2698112/101‑1 (MQ=255)

                                                    aaGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGc                                             >  1:2698916/1‑101 (MQ=255)

                                                     aGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCa                                            <  2:2698639/101‑1 (MQ=255)

                                                     aGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCa                                            <  1:2698570/101‑1 (MQ=255)

                                                     aGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCa                                            >  2:2698919/1‑101 (MQ=255)

                                                        ctcGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAAt                                         >  2:2698924/1‑101 (MQ=255)

                                                          cGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATAc                                       <  2:2698561/101‑1 (MQ=255)

                                                          cGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATAc                                       <  2:2698377/101‑1 (MQ=255)

                                                          cGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATAc                                       <  1:2698525/101‑1 (MQ=255)

                                                          cGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATAc                                       <  1:2698555/101‑1 (MQ=255)

                                                           gCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACg                                      <  2:2698433/101‑1 (MQ=255)

                                                            ccAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGa                                     >  1:2698927/1‑101 (MQ=255)

                                                              aaaCAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGAtt                                   <  1:2698465/101‑1 (MQ=255)

                                                              aaaCAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGAtt                                   >  2:2698930/1‑101 (MQ=255)

                                                               aaCAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGAttt                                  <  2:2698480/101‑1 (MQ=255)

                                                                 cAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCt                                <  2:2698446/101‑1 (MQ=255)

                                                                 cAGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCt                                >  1:2698933/1‑101 (MQ=255)

                                                                  aGCGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTg                               <  2:2698358/101‑1 (MQ=255)

                                                                    cGTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGAt                             >  2:2698934/1‑101 (MQ=255)

                                                                     gTTCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGAtg                            <  1:2698508/101‑1 (MQ=255)

                                                                       tCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGAtgtg                          <  2:2698136/101‑1 (MQ=255)

                                                                       tCAGGGCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGAtgtg                          >  1:2698938/1‑101 (MQ=255)

                                                                          gggCTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCAt                       >  1:2698941/1‑101 (MQ=255)

                                                                             cTATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTAc                    <  2:2698424/101‑1 (MQ=255)

                                                                              tATGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACt                   >  2:2698944/1‑101 (MQ=255)

                                                                               aTGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACtt                  <  1:2698240/101‑1 (MQ=255)

                                                                                tGACTACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTa                 <  2:2698604/101‑1 (MQ=255)

                                                                                    tACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCg             <  1:2698414/101‑1 (MQ=255)

                                                                                    tACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCg             >  1:2698948/1‑101 (MQ=255)

                                                                                    tACGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCg             <  1:2698659/101‑1 (MQ=255)

                                                                                     aCGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCgg            <  1:2698606/101‑1 (MQ=255)

                                                                                     aCGAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCgg            >  2:2698950/1‑101 (MQ=255)

                                                                                       gAAGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAg          >  2:2698954/1‑101 (MQ=255)

                                                                                        aaGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGt         <  1:2698540/101‑1 (MQ=255)

                                                                                        aaGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGt         >  1:2698955/1‑101 (MQ=255)

                                                                                        aaGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGt         <  1:2698365/101‑1 (MQ=255)

                                                                                         aGCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTa        <  1:2698435/101‑1 (MQ=255)

                                                                                          gCGGATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTaa       >  1:2698958/1‑101 (MQ=255)

                                                                                            ggATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGt     <  1:2698560/101‑1 (MQ=255)

                                                                                             gATGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGtt    >  2:2698962/1‑101 (MQ=255)

                                                                                              atggcaatggcaTCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTa   <  2:2698741/101‑1 (MQ=255)

                                                                                              atggcaatggcaTCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTa   <  1:2698453/101‑1 (MQ=255)

                                                                                               tGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTAt  >  1:2698966/1‑101 (MQ=255)

                                                                                               tGGCAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTAt  <  2:2698211/101‑1 (MQ=255)

GGATCAATGAAAACGGGCGACAAGCTGGATATCAAGGGCGTAGGTCAGGTTGAAGTCTCGCCAAACAGCGTTCAGGGCTATGACTACGAAGCGGATGGTAATGGCATCGTACTGTTACCGGAGCGCGTGATATTCAACAAAGAGAATATCGGCAAATACGATTTCTGATGTGCATTACTTAACCGGAGTAAGTTAT  >  NC_000913/1605906‑1606101

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