BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,063,840	C→T	G312D (GGC→GAC)	b3877 ←	putative permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,063,840	0	C	T	100.0%	833.9 / NA	229	G312D (GGC→GAC)	b3877	putative permease
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (106/123); total (106/123)

Only 100 of 237 total aligned reads displayed.

AACGCCAGGCAGGAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGCCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGAAAAATCCC  >  NC_000913/4063750‑4063929

aaCGCCAGGCAGGAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGccgtcg                                                                                         >  1:5813584/1‑93 (MQ=255)

   gCCAGGCAGGAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGata                                                                                      <  1:1016082/93‑1 (MQ=255)

   gCCAGGCAGGAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGata                                                                                      >  1:2225589/1‑93 (MQ=255)

   gCCAGGCAGGAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGata                                                                                      >  2:6030643/1‑93 (MQ=255)

     caggcaggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGAtata                                                                                    >  1:6155475/1‑93 (MQ=255)

       ggcaggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATaaa                                                                                  >  1:1271768/1‑93 (MQ=255)

       ggcaggAGAACGCCACGAAGCTCACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATaaa                                                                                  <  1:4828252/93‑1 (MQ=255)

        gcaggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAAc                                                                                 <  1:2672314/93‑1 (MQ=255)

         caggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACt                                                                                >  2:5183299/1‑93 (MQ=255)

         caggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACt                                                                                <  1:5598694/93‑1 (MQ=255)

          aggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACtt                                                                               <  2:6253024/93‑1 (MQ=255)

           ggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGCGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATCAACttt                                                                              >  2:3521735/1‑93 (MQ=255)

           ggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGCTATAAACttt                                                                              >  2:3923002/1‑93 (MQ=255)

           ggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACttt                                                                              <  2:5322279/93‑1 (MQ=255)

           ggAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACttt                                                                              >  2:3734319/1‑93 (MQ=255)

            gagaACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTttc                                                                             <  1:4973556/93‑1 (MQ=255)

            gagaACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTttc                                                                             >  2:916231/1‑93 (MQ=255)

             agaACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTttct                                                                            <  2:3090160/93‑1 (MQ=255)

                aCGCCACGAAGCTTACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCTCCCAAATCAGCAGGCCGTCGATATAAACTTtcttct                                                                         <  2:4368768/93‑1 (MQ=255)

                 cGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTcttctt                                                                        <  1:5335095/93‑1 (MQ=255)

                  gCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTg                                                                       <  1:5008193/93‑1 (MQ=255)

                  gCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTg                                                                       <  1:3849817/93‑1 (MQ=255)

                   ccACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGc                                                                      >  2:878985/1‑93 (MQ=255)

                       gAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCaaa                                                                  >  1:1191322/1‑93 (MQ=255)

                        aaGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCaaaa                                                                 <  2:3096647/93‑1 (MQ=255)

                        aaGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCaaaa                                                                 >  1:2943271/1‑93 (MQ=255)

                         aGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAAc                                                                >  1:553253/1‑93 (MQ=255)

                          gCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACg                                                               >  2:1267551/1‑93 (MQ=255)

                            tGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTc                                                             >  2:3349161/1‑93 (MQ=255)

                            tGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTc                                                             >  1:667137/1‑93 (MQ=255)

                            tGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTc                                                             >  2:5615794/1‑93 (MQ=255)

                             gACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTct                                                            <  1:1728492/93‑1 (MQ=255)

                             gACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTct                                                            >  2:4865419/1‑93 (MQ=255)

                              aCCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctg                                                           <  2:1597367/93‑1 (MQ=255)

                              aCCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctg                                                           <  2:1491643/93‑1 (MQ=255)

                              aCCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctg                                                           <  2:1211270/93‑1 (MQ=255)

                                  aaCCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctgactg                                                       >  1:3621736/1‑93 (MQ=255)

                                  aaCCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctgactg                                                       <  1:6233048/93‑1 (MQ=255)

                                   aCCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctgactga                                                      <  2:5659752/93‑1 (MQ=255)

                                   aCCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTctgactga                                                      <  2:3119558/93‑1 (MQ=255)

                                     cgccgccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGagg                                                    <  1:5248542/93‑1 (MQ=255)

                                     cgccgccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGagg                                                    <  1:947234/93‑1 (MQ=255)

                                     cgccgccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGagg                                                    <  1:5896726/93‑1 (MQ=255)

                                     cgccgccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGagg                                                    <  2:650261/93‑1 (MQ=255)

                                      gccgccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGaggc                                                   <  1:354212/93‑1 (MQ=255)

                                        cgccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGaggcag                                                 <  1:2765919/93‑1 (MQ=255)

                                         gccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGaggcagg                                                <  2:4180613/93‑1 (MQ=255)

                                         gccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGaggcagg                                                >  1:4343659/1‑93 (MQ=255)

                                         gccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGaggcagg                                                <  1:1943675/93‑1 (MQ=255)

                                          ccgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGaggcaggc                                               >  1:2581481/1‑93 (MQ=255)

                                           cgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAggcaggca                                              >  2:445779/1‑93 (MQ=255)

                                           cgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAggcaggca                                              >  2:2312082/1‑93 (MQ=255)

                                           cgAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAggcaggca                                              >  2:1195754/1‑93 (MQ=255)

                                             aagaaATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATc                                            >  1:6004954/1‑93 (MQ=255)

                                             aagaaATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATc                                            >  1:3280639/1‑93 (MQ=255)

                                              agaaATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCAGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCa                                           >  1:4441476/1‑93 (MQ=255)

                                                aaaTAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAgg                                         <  2:2499238/93‑1 (MQ=255)

                                                 aaTAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGa                                        <  2:3439330/93‑1 (MQ=255)

                                                 aaTAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGa                                        <  1:20890/93‑1 (MQ=255)

                                                  aTAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGaa                                       >  2:2763131/1‑93 (MQ=255)

                                                   tAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAAt                                      >  2:4651391/1‑93 (MQ=255)

                                                    aGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATa                                     <  2:4576342/93‑1 (MQ=255)

                                                     gTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATAc                                    <  2:3999740/93‑1 (MQ=255)

                                                     gTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATAc                                    <  2:3562446/93‑1 (MQ=255)

                                                        gAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGcc                                 <  2:557054/93‑1 (MQ=255)

                                                        gAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGcc                                 <  2:5117039/93‑1 (MQ=255)

                                                         aGCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCg                                <  1:2349599/93‑1 (MQ=255)

                                                          gCAGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGa                               <  1:5962940/93‑1 (MQ=255)

                                                            aGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAATCTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATg                             <  2:1998724/93‑1 (MQ=255)

                                                            aGATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATg                             >  2:4118867/1‑93 (MQ=255)

                                                             gATCGCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGa                            <  1:870764/93‑1 (MQ=255)

                                                                 gCCCAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAAt                        <  1:2752711/93‑1 (MQ=255)

                                                                 ccccAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAAt                        <  2:2433456/92‑1 (MQ=255)

                                                                  cccAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATa                       >  1:2813826/1‑93 (MQ=255)

                                                                   ccAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATAc                      >  1:3066094/1‑93 (MQ=255)

                                                                   ccAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATAc                      >  2:5734201/1‑93 (MQ=255)

                                                                   ccAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATAc                      <  2:5869505/93‑1 (MQ=255)

                                                                   ccAGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATAc                      <  1:4077268/93‑1 (MQ=255)

                                                                    cAGCACCCAAATCAGCAGGCCATCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACa                     >  1:2987763/1‑93 (MQ=255)

                                                                     aGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAg                    >  1:187842/1‑93 (MQ=255)

                                                                     aGCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAg                    <  1:4145459/93‑1 (MQ=255)

                                                                      gCACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGc                   <  1:5693249/93‑1 (MQ=255)

                                                                       cACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGcc                  >  2:1851604/1‑93 (MQ=255)

                                                                       cACCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGcc                  >  2:2277683/1‑93 (MQ=255)

                                                                        aCCCAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGccc                 >  1:486193/1‑93 (MQ=255)

                                                                          ccAAATCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCAt               <  2:4220785/93‑1 (MQ=255)

                                                                            aaaTCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGc             <  1:972436/93‑1 (MQ=255)

                                                                             aaTCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCt            <  1:2620765/93‑1 (MQ=255)

                                                                             aaTCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCt            >  2:6159110/1‑93 (MQ=255)

                                                                             aaTCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCt            >  1:3735269/1‑93 (MQ=255)

                                                                              aTCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTg           <  2:5367040/93‑1 (MQ=255)

                                                                               tCAGCAGGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGa          >  2:2371011/1‑93 (MQ=255)

                                                                               tCAGCAAGCCGTCGATGTGAACTTACTACAAGCCAAAACAACTGACTGAGGCCGGCATCAGGAATACGCCGATGAAAATACAGCCGATGCTGa          <  2:6309647/93‑1 (MQ=255)

                                                                                cagcagGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGaa         <  1:5052189/93‑1 (MQ=255)

                                                                                 agcagGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGaaa        >  2:3506700/1‑93 (MQ=255)

                                                                                 agcagGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGaaa        <  2:3205482/93‑1 (MQ=255)

                                                                                  gcagGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGaaaa       <  1:5189530/93‑1 (MQ=255)

                                                                                  gcagGCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGaaaa       >  1:4310231/1‑93 (MQ=255)

                                                                                     ggCCGTCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGAAAAATc    >  2:1813054/1‑93 (MQ=255)

                                                                                       ccgtcgATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGAAAAATccc  >  2:5996967/1‑93 (MQ=255)

AACGCCAGGCAGGAGAACGCCACGAAGCTGACCGAACCGCCGCCGAAGAAATAGTTGAGCAGATCGCCCAGCACCCAAATCAGCAGGCCGCCGATATAAACTTTCTTCTTGCCAAAACGTCTGACTGAGGCAGGCATCAGGAATACGCCGATGAAAATACAGCCCATGCTGAAAAATCCC  >  NC_000913/4063750‑4063929

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