BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,495,173	T→A	A40A (GCA→GCT)	yfdY ←	DUF2545 family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,495,173	0	T	A	100.0%	427.0 / NA	125	A40A (GCA→GCT)	yfdY	DUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (54/71); total (54/71)

Only 100 of 129 total aligned reads displayed.

GCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCG  >  NC_000913/2495077‑2495269

  gAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGagcagc                                                                                             <  1:752958/100‑1 (MQ=255)

     tcttcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAg                                                                                          <  2:904173/97‑1 (MQ=255)

      catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                         >  2:2406203/1‑100 (MQ=255)

      catcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgg                                                                                         >  2:3079256/1‑100 (MQ=255)

       atcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                        <  1:1327220/100‑1 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCCTTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       <  2:3223184/100‑1 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAgaga                                                                                       >  1:904161/1‑97 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       <  1:3116449/100‑1 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       >  2:1564156/1‑100 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       >  2:148505/1‑100 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       >  2:3120686/1‑100 (MQ=255)

        tcaAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGAt                                                                                       >  2:3477427/1‑100 (MQ=255)

        tcaAGCCACAAGACGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGa                                                                                        <  1:929977/99‑1 (MQ=255)

         caAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATc                                                                                      >  1:2358725/1‑100 (MQ=255)

          aaGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCa                                                                                     >  2:177348/1‑100 (MQ=255)

          aaGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCa                                                                                     >  1:2531048/1‑100 (MQ=255)

          aaGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCa                                                                                     >  1:2415444/1‑100 (MQ=255)

           aGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAc                                                                                    <  2:655781/100‑1 (MQ=255)

            gCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCAcc                                                                                   <  2:2336949/100‑1 (MQ=255)

             ccACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCa                                                                                  <  2:757892/100‑1 (MQ=255)

                caAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGc                                                                               <  1:652803/100‑1 (MQ=255)

                caAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGc                                                                               <  2:2890440/100‑1 (MQ=255)

                 aaGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGcc                                                                              <  2:149915/100‑1 (MQ=255)

                  atccGGGGCAAGAAACATTAAGATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATTGTGATGTAATAAATGAGAGCATCGAGGATCACCATGCCg                                                                             <  1:477725/98‑1 (MQ=255)

                  aGCCGGGGCAAGAAACATAAAGATGCCCGGCACTAGTGCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGCGAGCAGCGAGGATCACCATGcc                                                                              <  1:2021358/99‑1 (MQ=255)

                    ccGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaa                                                                           <  2:186614/100‑1 (MQ=255)

                     cGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaag                                                                          <  2:1888001/100‑1 (MQ=255)

                       gggCAAGAAACATCAATATCCGCGGCACTAGCTCATTACCGGCAGTTAACCACATGGGGAAGTAATAGACGAGAGCAGCGAGGACCACCATGc                                                                               <  1:2375510/93‑1 (MQ=255)

                       gggCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaagaa                                                                        <  2:1543254/100‑1 (MQ=255)

                       gggCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCgaagaa                                                                        <  1:360355/100‑1 (MQ=255)

                         gCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGaagaaa                                                                       <  2:3056027/99‑1 (MQ=255)

                          cAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGTAGCGAGGATCACCATGCCGaagaaaga                                                                     <  1:1854766/100‑1 (MQ=255)

                            aGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAaagaag                                                                   <  2:2156088/100‑1 (MQ=255)

                            aGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAaagaag                                                                   <  1:40849/100‑1 (MQ=255)

                              aaaCATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAagaagan                                                                 <  1:2644833/100‑2 (MQ=255)

                                 cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAcc                                                              >  1:2683781/1‑100 (MQ=255)

                                 cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAcc                                                              <  1:3224171/100‑1 (MQ=255)

                                 cATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGAAATAAATGAGAGCAGCGAGGAACACCATGCCGAAGAAAGAAGAGAcc                                                              >  1:1804586/1‑100 (MQ=255)

                                  aTAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGAccg                                                             >  2:942010/1‑100 (MQ=255)

                                       tatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGAt                                                        <  1:2007770/100‑1 (MQ=255)

                                       tatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGAt                                                        <  2:238618/100‑1 (MQ=255)

                                       tatCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGAt                                                        <  1:2007633/100‑1 (MQ=255)

                                        atCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATa                                                       >  2:261649/1‑100 (MQ=255)

                                         tCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATAc                                                      <  1:1625345/100‑1 (MQ=255)

                                           cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    <  2:3312226/100‑1 (MQ=255)

                                           cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    <  1:173449/100‑1 (MQ=255)

                                           cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    <  2:3315356/100‑1 (MQ=255)

                                           cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    <  1:2713594/100‑1 (MQ=255)

                                           cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    >  2:1739801/1‑100 (MQ=255)

                                           cccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCa                                                    >  2:173451/1‑100 (MQ=255)

                                            ccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCAc                                                   <  2:1497983/100‑1 (MQ=255)

                                            ccGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCAc                                                   <  2:631848/100‑1 (MQ=255)

                                              gTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGt                                                 <  2:3406657/100‑1 (MQ=255)

                                                cACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTc                                               <  2:2531082/100‑1 (MQ=255)

                                                cACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTc                                               <  1:148329/100‑1 (MQ=255)

                                                    aGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTa                                           >  2:2719357/1‑100 (MQ=255)

                                                     gTTCATTACCGTCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTAc                                          >  1:1012694/1‑100 (MQ=255)

                                                       tCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCt                                        <  1:1021044/100‑1 (MQ=255)

                                                        cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       >  2:3383028/1‑100 (MQ=255)

                                                        cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       >  1:2293596/1‑100 (MQ=255)

                                                        cATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTg                                       <  2:2443206/100‑1 (MQ=255)

                                                           tcttccgatctAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGAcc                                    <  2:909029/90‑1 (MQ=255)

                                                           tACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGAcc                                    >  1:1232716/1‑100 (MQ=255)

                                                           tACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGAcc                                    <  1:523954/100‑1 (MQ=255)

                                                            aCCGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCg                                   >  2:2802464/1‑100 (MQ=255)

                                                             ccGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                  >  2:1434648/1‑100 (MQ=255)

                                                              cGCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGa                                  <  1:112074/99‑1 (MQ=255)

                                                               gCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGTGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATg                                >  1:1642666/1‑100 (MQ=255)

                                                               gCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCATTACCTGACCGATg                                <  2:2358757/100‑1 (MQ=255)

                                                               gCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATg                                >  2:807033/1‑100 (MQ=255)

                                                               gCCAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATg                                >  1:600798/1‑100 (MQ=255)

                                                                 cAGTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTa                              <  1:3382976/100‑1 (MQ=255)

                                                                   gTTAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGc                            >  1:3245641/1‑100 (MQ=255)

                                                                    ttAACCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGcc                           >  2:2437797/1‑100 (MQ=255)

                                                                      aaCCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGc                         <  1:1564134/100‑1 (MQ=255)

                                                                      aaCCACATGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGc                         >  1:3097746/1‑100 (MQ=255)

                                                                         cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                      >  1:558519/1‑100 (MQ=255)

                                                                         cacaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTa                      >  2:3407624/1‑100 (MQ=255)

                                                                          acaTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTAc                     <  2:2311121/100‑1 (MQ=255)

                                                                            aTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGc                   <  2:283288/100‑1 (MQ=255)

                                                                            aTGGTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGc                   <  2:1946136/100‑1 (MQ=255)

                                                                              ggTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAg                 <  1:1105205/100‑1 (MQ=255)

                                                                              ggTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAg                 >  1:1445629/1‑100 (MQ=255)

                                                                              ggTGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAg                 >  2:531478/1‑100 (MQ=255)

                                                                                tGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGAc               <  1:3407572/100‑1 (MQ=255)

                                                                                tGAAGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGAc               <  2:2862107/100‑1 (MQ=255)

                                                                                 ctaGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACa              <  2:762031/98‑1 (MQ=255)

                                                                                  aaGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACaa             >  1:1891711/1‑100 (MQ=255)

                                                                                   aGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACaag            >  1:762020/1‑99 (MQ=255)

                                                                                   aGTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAAt            <  1:2510052/100‑1 (MQ=255)

                                                                                    gTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAATGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATa           <  1:2437764/100‑1 (MQ=255)

                                                                                    gTAATAAATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATa           <  2:3245693/100‑1 (MQ=255)

                                                                                     taataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATac          >  1:3379781/1‑100 (MQ=255)

                                                                                       ataaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATacac        <  1:894121/100‑1 (MQ=255)

                                                                                        taaATGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATAcaca       >  2:1546335/1‑100 (MQ=255)

                                                                                         aaaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATAcaca       >  1:926231/1‑99 (MQ=255)

                                                                                          aaTGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAgg     <  2:1071536/100‑1 (MQ=255)

                                                                                            tGAGAGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGc   >  1:820881/1‑100 (MQ=255)

                                                                                             gagaGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCg  >  2:1342342/1‑100 (MQ=255)

                                                                                             gagaGCAGCGAGGAGCACCATGCCGTAGAAAGAAGCGACCGCCTATACCACGTCCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGc   <  2:1253099/99‑1 (MQ=255)

GCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAATAAATGAGTGCAGCGAGGATCACCATGCCGAAGAAAGAAGAGACCGCCGATACCACGTTCAGTACCTGACCGATGTAGCCGCTTACGCAGACAATACACAGGGCG  >  NC_000913/2495077‑2495269

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