BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,992,862	A→G	Q164R (CAG→CGG)	dapF →	diaminopimelate epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,992,862	0	A	G	100.0%	650.1 / NA	209	Q164R (CAG→CGG)	dapF	diaminopimelate epimerase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (104/105); total (104/105)

Only 100 of 220 total aligned reads displayed.

AAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTT  >  NC_000913/3992774‑3992951

aaaGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTc                                                                                       >  1:4575049/1‑93 (MQ=255)

   gCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgat                                                                                    >  2:1496326/1‑93 (MQ=255)

   gCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgat                                                                                    >  1:1180412/1‑93 (MQ=255)

     ggAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatga                                                                                  >  1:2270269/1‑93 (MQ=255)

     ggAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatga                                                                                  >  2:1774658/1‑93 (MQ=255)

      gAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCgatgat                                                                                 >  1:444033/1‑93 (MQ=255)

       aaaaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGatgatg                                                                                <  1:4144589/93‑1 (MQ=255)

        aaaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGt                                                                               <  2:2733204/93‑1 (MQ=255)

         aaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCTTGATTCGGGTCGATGATGTc                                                                              <  2:4513759/93‑1 (MQ=255)

         aaGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTc                                                                              >  2:4443853/1‑93 (MQ=255)

          aGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGTATTGCGTGATTCGGGTCGATGATGTCg                                                                             <  1:2661215/93‑1 (MQ=255)

           gACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGa                                                                            >  2:3269754/1‑93 (MQ=255)

            gccTATATTATGCGCGTCGGCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGTAAATCCCCATTGCGGGATTCGGGTCGATGATGTCGAt                                                                           <  2:651888/92‑1 (MQ=255)

            aCCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGaa                                                                           >  1:2237331/1‑92 (MQ=255)

             ccTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCACTGCGTGATTCGGGTCGCTGATGTCGATa                                                                          >  2:318164/1‑93 (MQ=255)

              cTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATAc                                                                         >  2:838944/1‑93 (MQ=255)

               tataTTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATAcc                                                                        >  2:2029943/1‑93 (MQ=255)

                 tattatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgc                                                                      <  2:1120520/93‑1 (MQ=255)

                  attatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgcg                                                                     >  1:2827271/1‑93 (MQ=255)

                  attatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACcgcg                                                                     >  2:3592432/1‑93 (MQ=255)

                    tatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCgcggc                                                                   >  1:3255463/1‑93 (MQ=255)

                    tatGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCgcggc                                                                   >  1:2981370/1‑93 (MQ=255)

                       gcgcgcCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGt                                                                >  2:1956403/1‑93 (MQ=255)

                          cgccgcCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTata                                                             <  2:4086985/93‑3 (MQ=255)

                           gccgccGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGaa                                                            <  2:1907517/93‑1 (MQ=255)

                            ccgccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGaag                                                           <  1:3396985/93‑2 (MQ=255)

                             cgccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAgcagc                                                          >  1:318828/1‑89 (MQ=255)

                             cgccgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGaagc                                                          <  1:1348420/93‑3 (MQ=255)

                               ccgAGCAGACAGTCTTATGCGGCGTGGTGTCGATAGAAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATCCCGCGGCGGTAGAAgcgc                                                        <  2:3319067/93‑1 (MQ=255)

                                cgAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCt                                                       >  2:2566355/1‑93 (MQ=255)

                                 gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                      <  2:1311447/93‑1 (MQ=255)

                                 gAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCtt                                                      >  2:1519641/1‑93 (MQ=255)

                                    cAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGt                                                   <  1:296783/93‑1 (MQ=255)

                                    cAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGt                                                   <  2:2955804/93‑1 (MQ=255)

                                     aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTc                                                  <  2:2535665/93‑1 (MQ=255)

                                     aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTc                                                  <  1:1501045/93‑1 (MQ=255)

                                      gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                 >  1:539785/1‑93 (MQ=255)

                                      gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                 >  2:352609/1‑93 (MQ=255)

                                      gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTcc                                                 >  2:3794960/1‑93 (MQ=255)

                                       aCAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCt                                                <  1:2204575/93‑1 (MQ=255)

                                        cAATCTTATGCGGCGTGGTGTCGATGGTAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTg                                               <  1:655578/93‑1 (MQ=255)

                                         aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                              <  1:1782363/93‑1 (MQ=255)

                                         aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                              <  2:2192576/93‑1 (MQ=255)

                                         aaTCTTAGGCGGCGTTGTGTCGATGGGAAATCCGCCTTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGt                                              <  1:55036/93‑1 (MQ=255)

                                          aTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGtt                                             >  1:3335261/1‑93 (MQ=255)

                                           tCTTATGCGGCGTGGTGTCGATGGTAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                            >  1:169841/1‑93 (MQ=255)

                                           tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                            >  2:238925/1‑93 (MQ=255)

                                           tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                            <  2:2406717/93‑1 (MQ=255)

                                           tCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTc                                            >  1:1054836/1‑93 (MQ=255)

                                            cTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCt                                           <  2:2931941/93‑1 (MQ=255)

                                             ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                          <  1:491688/93‑1 (MQ=255)

                                             ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                          >  2:2785593/1‑93 (MQ=255)

                                             ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                          <  1:2870309/93‑1 (MQ=255)

                                             ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTg                                          <  1:4050459/93‑1 (MQ=255)

                                                tGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaa                                       <  2:1745152/93‑1 (MQ=255)

                                                tGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGCCGGGAGAAGCGCTTGGTCCTGTTCTGGaa                                       >  1:2080810/1‑93 (MQ=255)

                                                 gcggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGaaa                                      <  1:4282405/93‑1 (MQ=255)

                                                   ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                    >  2:2206543/1‑93 (MQ=255)

                                                   ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGGAGAAGCGCTTGGTCCTGTTCTGGAAAGc                                    >  1:37578/1‑93 (MQ=255)

                                                    gcgTGGTGTCGATGGGAAATCCGCACTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGcc                                   >  2:3959961/1‑93 (MQ=255)

                                                     cgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCa                                  >  1:954656/1‑93 (MQ=255)

                                                      gtggtgTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCAc                                 <  2:781964/93‑1 (MQ=255)

                                                      gtggtgTCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCAc                                 <  2:3409869/93‑1 (MQ=255)

                                                        tgtgtCGATGGGAAATCCGCATTGCGTGATTCGGGGCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGa                               <  1:67660/92‑1 (MQ=255)

                                                         gtgtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAg                              <  2:3831249/93‑1 (MQ=255)

                                                         gtgtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAg                              >  1:3499459/1‑93 (MQ=255)

                                                           gtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCg                            <  2:230916/93‑1 (MQ=255)

                                                           gtCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCAACGacct                            >  1:2584749/1‑90 (MQ=255)

                                                           gtCGATGCGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCg                            <  1:2314554/93‑1 (MQ=255)

                                                            tCGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                           <  1:1016340/93‑1 (MQ=255)

                                                            gcGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGt                           <  2:1455295/92‑1 (MQ=255)

                                                             cGATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGtt                          <  1:1667991/93‑1 (MQ=255)

                                                              gATGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                         <  1:2854782/93‑1 (MQ=255)

                                                              gATGGGAAATCCGCATTGCGTGATTCGGGGCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGttt                         <  2:643062/93‑1 (MQ=255)

                                                                tGGGAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGTAAAGCCACGAGCGTTTTc                       <  1:1459610/93‑1 (MQ=255)

                                                                 gggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                      >  1:1586164/1‑93 (MQ=255)

                                                                 gggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTcc                      <  1:1141635/93‑1 (MQ=255)

                                                                  ggAAATCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCg                     >  1:395359/1‑93 (MQ=255)

                                                                    aaaTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGa                   <  2:2531624/93‑1 (MQ=255)

                                                                      aTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                 <  1:4006534/93‑1 (MQ=255)

                                                                      aTCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgc                 <  2:1725919/93‑1 (MQ=255)

                                                                       tCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgtg                >  2:123181/1‑91 (MQ=255)

                                                                       tCCGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcg                >  2:2087631/1‑93 (MQ=255)

                                                                        ccGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc               <  1:1420926/93‑1 (MQ=255)

                                                                        ccGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc               <  1:4343613/93‑1 (MQ=255)

                                                                        ccGCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgc               <  2:715247/93‑1 (MQ=255)

                                                                          gCATTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAgcgcgc             >  2:494297/1‑93 (MQ=255)

                                                                            aTTGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCa           >  1:723582/1‑93 (MQ=255)

                                                                             atGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGGGTTTTCCGGAGCGCGCCaa          <  2:4583257/92‑1 (MQ=255)

                                                                              tGCGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAat         <  1:249611/93‑1 (MQ=255)

                                                                                cGTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAatat       <  1:2597833/93‑1 (MQ=255)

                                                                                 gTGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATc      <  1:3039460/93‑1 (MQ=255)

                                                                                  tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg     <  2:3776252/93‑1 (MQ=255)

                                                                                  tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg     >  2:537892/1‑93 (MQ=255)

                                                                                  tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg     <  1:364833/93‑1 (MQ=255)

                                                                                  tGATTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCg     <  1:709157/93‑1 (MQ=255)

                                                                                    aTTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt   <  1:1983960/93‑1 (MQ=255)

                                                                                    aTTCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGt   <  2:140940/93‑1 (MQ=255)

                                                                                     ttCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt  <  1:1860384/93‑1 (MQ=255)

                                                                                     ttCGGGTCGATGATGTCGATACCGCGGCGGTAGAAGCGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGtt  <  1:1100488/93‑1 (MQ=255)

AAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTT  >  NC_000913/3992774‑3992951

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