BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	486,646	G→A	S37S (TCG→TCA)	mscK →	potassium dependent, small conductance mechanosensitive channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	486,646	0	G	A	100.0%	1704.2 / NA	395	S37S (TCG→TCA)	mscK	potassium dependent, small conductance mechanosensitive channel
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (221/174); total (221/174)

Only 100 of 401 total aligned reads displayed.

CAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGG  >  NC_000913/486548‑486743

cAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCaaa                                                                                                 <  2:949491/101‑1 (MQ=255)

   tATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATgg                                                                                              >  2:950118/1‑101 (MQ=255)

   tATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATgg                                                                                              >  1:950115/1‑101 (MQ=255)

    aTTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGt                                                                                             <  1:949340/101‑1 (MQ=255)

    aTTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGt                                                                                             <  1:949039/101‑1 (MQ=255)

     ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATgg                                                                                              >  1:950125/1‑99 (MQ=255)

     ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTg                                                                                            <  2:949258/101‑1 (MQ=255)

     ttACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTg                                                                                            <  1:949738/101‑1 (MQ=255)

      tACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGa                                                                                           <  2:949381/101‑1 (MQ=255)

      tACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGa                                                                                           >  2:950127/1‑101 (MQ=255)

        cAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATc                                                                                         >  2:950131/1‑101 (MQ=255)

        cAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATc                                                                                         <  1:949725/101‑1 (MQ=255)

        cAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATc                                                                                         <  2:949744/101‑1 (MQ=255)

         aaaCGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCt                                                                                        <  1:949646/101‑1 (MQ=255)

         aaaCGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCt                                                                                        <  1:949786/101‑1 (MQ=255)

          aaCGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTg                                                                                       >  1:950137/1‑101 (MQ=255)

            cGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGcc                                                                                     <  2:949549/101‑1 (MQ=255)

            cGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGcc                                                                                     >  2:950141/1‑101 (MQ=255)

             gATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCg                                                                                    <  2:949500/101‑1 (MQ=255)

              aTCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGa                                                                                   >  1:950144/1‑101 (MQ=255)

              aTCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGa                                                                                   >  1:950149/1‑101 (MQ=255)

               tCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGAc                                                                                  >  2:950151/1‑101 (MQ=255)

                cACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACa                                                                                 <  1:949684/101‑1 (MQ=255)

                cACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACa                                                                                 >  1:950155/1‑101 (MQ=255)

                 aCGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACaa                                                                                >  2:950158/1‑101 (MQ=255)

                   ggCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACaaaa                                                                              >  1:950162/1‑101 (MQ=255)

                    gCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAg                                                                             >  1:950166/1‑101 (MQ=255)

                    gCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAg                                                                             <  1:949569/101‑1 (MQ=255)

                     cATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGcc                                                                                     >  2:950172/1‑92 (MQ=255)

                     cATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGc                                                                            <  1:949523/101‑1 (MQ=255)

                     cATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGc                                                                            <  2:949721/101‑1 (MQ=255)

                      aTTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCg                                                                           >  1:950176/1‑101 (MQ=255)

                       ttttGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCgg                                                                          >  1:950181/1‑101 (MQ=255)

                       ttttGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGAgg                                                                          <  1:949465/101‑1 (MQ=255)

                         ttGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGAc                                                                        >  2:950186/1‑101 (MQ=255)

                           gTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCt                                                                      >  1:950191/1‑101 (MQ=255)

                            ttttttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTg                                                                     >  2:950195/1‑101 (MQ=255)

                            ttttttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTg                                                                     <  2:949679/101‑1 (MQ=255)

                             tttttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGc                                                                    <  1:949804/101‑1 (MQ=255)

                              ttttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCa                                                                   >  2:950202/1‑101 (MQ=255)

                              ttttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCa                                                                   >  1:950198/1‑101 (MQ=255)

                               tttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAg                                                                  <  1:949675/101‑1 (MQ=255)

                                ttCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAgg                                                                 >  1:950207/1‑101 (MQ=255)

                                     cATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCaa                                                            >  1:950213/1‑101 (MQ=255)

                                     cATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCaa                                                            >  2:950214/1‑101 (MQ=255)

                                      aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAAc                                                           >  2:950222/1‑101 (MQ=255)

                                      aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAAc                                                           <  1:949666/101‑1 (MQ=255)

                                      aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAAc                                                           >  1:950221/1‑101 (MQ=255)

                                      aTTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAAc                                                           <  1:949606/101‑1 (MQ=255)

                                        ttATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGgcgc                                                              >  2:950232/1‑96 (MQ=255)

                                        ttATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACtt                                                         >  1:950229/1‑101 (MQ=255)

                                        ttATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACtt                                                         <  2:949750/101‑1 (MQ=255)

                                        ttATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACtt                                                         <  1:949699/101‑1 (MQ=255)

                                        ttATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACt                                                          <  2:949711/100‑1 (MQ=255)

                                         tATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTg                                                        >  1:950234/1‑101 (MQ=255)

                                             gCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                    >  2:950241/1‑101 (MQ=255)

                                             gCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTc                                                    <  2:949608/101‑1 (MQ=255)

                                              cTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCa                                                   >  1:950243/1‑101 (MQ=255)

                                               ttttgtttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCAc                                                  >  2:950246/1‑101 (MQ=255)

                                                    tttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaat                                             <  2:949584/101‑1 (MQ=255)

                                                    tttttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaat                                             >  2:950252/1‑101 (MQ=255)

                                                      tttgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaataa                                           >  2:950256/1‑101 (MQ=255)

                                                       ttgtCTTGTTATTCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACtaaataaa                                          <  2:949574/101‑1 (MQ=255)

                                                        tgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAAc                                         >  1:950262/1‑101 (MQ=255)

                                                        tgtCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAAc                                         >  2:950261/1‑101 (MQ=255)

                                                          tCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaa                                       >  2:950268/1‑101 (MQ=255)

                                                           cTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaa                                      >  1:950274/1‑101 (MQ=255)

                                                           cTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaa                                      >  2:950273/1‑101 (MQ=255)

                                                            ttGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaaa                                     <  2:949647/101‑1 (MQ=255)

                                                            ttGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaaa                                     <  1:949968/101‑1 (MQ=255)

                                                            ttGTTATGACAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaaa                                     >  2:950280/1‑101 (MQ=255)

                                                             tGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACaaaaa                                    >  1:950282/1‑101 (MQ=255)

                                                              gTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAg                                   <  2:949777/101‑1 (MQ=255)

                                                               ttATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGa                                  >  2:950286/1‑101 (MQ=255)

                                                                tATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGAt                                 <  2:949310/101‑1 (MQ=255)

                                                                tATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGAt                                 >  1:950291/1‑101 (MQ=255)

                                                                tATGCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGAt                                 <  2:949689/101‑1 (MQ=255)

                                                                   gCCAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCtt                              <  2:949845/101‑1 (MQ=255)

                                                                    ccAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCttt                             >  1:950297/1‑101 (MQ=255)

                                                                     cAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACACAAAGATCTTTc                            >  1:950301/1‑101 (MQ=255)

                                                                     cAGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTc                            <  2:949637/101‑1 (MQ=255)

                                                                      aGAACACGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCt                           >  2:950302/1‑101 (MQ=255)

                                                                         acacGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCt                        >  1:950308/1‑101 (MQ=255)

                                                                           acGGCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCa                      >  2:950311/1‑101 (MQ=255)

                                                                             ggCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAgg                    <  2:949733/101‑1 (MQ=255)

                                                                             ggCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAgg                    >  1:950315/1‑101 (MQ=255)

                                                                              gCGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGa                   >  2:950321/1‑101 (MQ=255)

                                                                               cGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGAc                  <  2:949433/101‑1 (MQ=255)

                                                                               cGTTTGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGAc                  >  1:950322/1‑101 (MQ=255)

                                                                                 tttGCGCGGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACaa                <  2:949702/101‑1 (MQ=255)

                                                                                      gcgGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTgg           <  1:949550/101‑1 (MQ=255)

                                                                                      gcgGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTgg           >  2:950328/1‑101 (MQ=255)

                                                                                       cgGGCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGt          <  2:949918/101‑1 (MQ=255)

                                                                                         ggCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTgc        >  1:950331/1‑101 (MQ=255)

                                                                                          gCGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTgca       <  2:950070/101‑1 (MQ=255)

                                                                                           cGTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTgcag      >  1:950337/1‑101 (MQ=255)

                                                                                            gTCATCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTg            <  1:949444/94‑1 (MQ=255)

                                                                                             tcatcaAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTgcagca    <  1:949847/101‑1 (MQ=255)

                                                                                               aTCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAgg  <  1:949768/101‑1 (MQ=255)

                                                                                               aTCAAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAgg  <  1:949866/101‑1 (MQ=255)

CAGTATTACAAACGATCACGGCATTTTGTTTTTTCAGCATTTATTGCTTTTGTTTTTGTCTTGTTATGCCAGAACACGGCGTTTGCGCGGGCGTCATCGAATGGTGATCTGCCGACAAAAGCGGACCTGCAGGCGCAACTTGACTCACTAAATAAACAAAAAGATCTTTCTGCTCAGGACAAACTGGTGCAGCAGG  >  NC_000913/486548‑486743

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: `‑`