BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,557,200	T→C	V137A (GTG→GCG)	b2442 →	putative prophage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,557,200	0	T	C	100.0%	783.2 / NA	241	V137A (GTG→GCG)	b2442	putative prophage integrase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (120/121); total (120/121)

Only 100 of 249 total aligned reads displayed.

CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289

cccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCgg                                                                                         >  1:3718380/1‑93 (MQ=255)

cccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCgg                                                                                         <  2:3192397/93‑1 (MQ=255)

  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       <  2:416865/93‑1 (MQ=255)

   ttATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGcggcgg                                                                                      >  1:6684013/1‑93 (MQ=255)

    tATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGa                                                                                     >  1:4252916/1‑93 (MQ=255)

    tATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGa                                                                                     <  2:6051481/93‑1 (MQ=255)

     aTTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGaa                                                                                    >  1:262854/1‑93 (MQ=255)

        gCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGAc                                                                                 <  1:6003684/93‑1 (MQ=255)

        gCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGAc                                                                                 <  1:1574151/93‑1 (MQ=255)

         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACt                                                                                >  1:510272/1‑93 (MQ=255)

         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACt                                                                                <  1:5144934/93‑1 (MQ=255)

          cAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTg                                                                               <  1:838773/93‑1 (MQ=255)

           aaTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTgg                                                                              >  2:73207/1‑93 (MQ=255)

            aTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGt                                                                             >  2:6413968/1‑93 (MQ=255)

            aTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGt                                                                             >  2:2718512/1‑93 (MQ=255)

               gAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTg                                                                          <  2:4597415/93‑1 (MQ=255)

                aaTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGc                                                                         >  2:6255757/1‑93 (MQ=255)

                aaTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGc                                                                         <  2:4642700/93‑1 (MQ=255)

                aaTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGc                                                                         >  1:3595928/1‑93 (MQ=255)

                 aTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCt                                                                        <  1:3666381/93‑1 (MQ=255)

                   tAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAc                                                                      <  1:2149788/93‑1 (MQ=255)

                    aaCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAcc                                                                     <  2:3207116/93‑1 (MQ=255)

                    aaCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAcc                                                                     <  2:1561031/93‑1 (MQ=255)

                    aaCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTAcc                                                                     <  1:4075582/93‑1 (MQ=255)

                      cgcAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAg                                                                   <  2:2769365/91‑1 (MQ=255)

                      cccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAg                                                                   >  2:3695292/1‑93 (MQ=255)

                      cccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAg                                                                   <  2:5904302/93‑1 (MQ=255)

                      cccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAg                                                                   <  1:1943994/93‑1 (MQ=255)

                       ccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTTCTACCAgg                                                                  >  2:900678/1‑93 (MQ=255)

                         aGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGt                                                                <  2:6599329/93‑1 (MQ=255)

                          gTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAgggct                                                               >  1:3320104/1‑91 (MQ=255)

                          gTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGtt                                                               <  1:3220943/93‑1 (MQ=255)

                           tGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTg                                                              <  2:1250726/93‑1 (MQ=255)

                           tGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTg                                                              >  1:6241129/1‑93 (MQ=255)

                            ggAACAAGCCCGCACGGAAAAAGCCATTGATGTGCTCAACATGGCGGCCGGATTTAAGCAGGCGGCGGAAGACCGGTTTGCCACCAGGGATgg                                                             >  1:1358632/1‑93 (MQ=255)

                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             >  1:5718379/1‑93 (MQ=255)

                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             >  2:6706881/1‑93 (MQ=255)

                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             <  2:3647502/93‑1 (MQ=255)

                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGa                                                            >  2:2296908/1‑93 (MQ=255)

                              aaCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGGTGGAg                                                           >  1:5316710/1‑93 (MQ=255)

                                cAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAggg                                                         <  2:5601842/93‑1 (MQ=255)

                                 aaGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggt                                                        <  2:336025/93‑1 (MQ=255)

                                  aGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggtg                                                       >  1:973878/1‑93 (MQ=255)

                                   gCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggtgg                                                      <  1:1260384/93‑1 (MQ=255)

                                   gCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggtgg                                                      >  1:664468/1‑93 (MQ=255)

                                    cccGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGggtggt                                                     <  1:3544886/93‑1 (MQ=255)

                                     ccGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTc                                                    >  1:6134073/1‑93 (MQ=255)

                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  1:2497411/1‑93 (MQ=255)

                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  2:5063943/1‑93 (MQ=255)

                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGGGGTCGGAAt                                               >  1:1689122/1‑93 (MQ=255)

                                            gAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATcc                                             <  1:4460576/93‑1 (MQ=255)

                                            gAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATcc                                             <  2:5660647/93‑1 (MQ=255)

                                            gAAAAAGCCATTGATGCGCTCAACATGGCGCCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATcc                                             <  2:653888/93‑1 (MQ=255)

                                             aaaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCt                                            >  2:5545961/1‑93 (MQ=255)

                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           >  2:4729750/1‑93 (MQ=255)

                                               aaaGCCATTGATGCGCTCAGCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTAt                                          <  2:4489608/93‑1 (MQ=255)

                                               aaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTAt                                          >  1:5008541/1‑93 (MQ=255)

                                               aaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGCTGGAGGGTGGTCGGAATCCTAt                                          >  1:2725568/1‑93 (MQ=255)

                                                aaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGGGGTCGGAATCCTATg                                         >  1:5591782/1‑93 (MQ=255)

                                                 aGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGACTCCTATGc                                        >  2:3007964/1‑93 (MQ=255)

                                                 aGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATACTATGc                                        >  1:4661160/1‑93 (MQ=255)

                                                   ccATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaa                                      >  1:3607081/1‑93 (MQ=255)

                                                    cATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaa                                     >  1:3458726/1‑93 (MQ=255)

                                                    cATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaa                                     >  2:1787434/1‑93 (MQ=255)

                                                     aTTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaaa                                    >  2:694740/1‑93 (MQ=255)

                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAAc                                   >  2:3973228/1‑93 (MQ=255)

                                                        gATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAg                                 >  2:4743303/1‑93 (MQ=255)

                                                          tGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTCTGCAAAACAggc                               >  2:2548812/1‑92 (MQ=255)

                                                          tGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGt                               <  2:517473/93‑1 (MQ=255)

                                                           gcgcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTc                              >  1:489830/1‑93 (MQ=255)

                                                            cgcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCa                             <  2:2061782/93‑1 (MQ=255)

                                                             gcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGGGGGAATCCTATCCAAAACAGGTCAg                            >  2:1648642/1‑93 (MQ=255)

                                                             gcTCAAAATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAg                            >  1:5552870/1‑93 (MQ=255)

                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           >  1:6007988/1‑93 (MQ=255)

                                                               tCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGAt                          >  1:5473961/1‑93 (MQ=255)

                                                               tCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGAt                          >  2:5367459/1‑93 (MQ=255)

                                                                cAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATc                         <  1:1070644/93‑1 (MQ=255)

                                                                 aaCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGCGGTCGGAATCCTATGCAAAACAGGTCAGATCg                        >  1:6482029/1‑93 (MQ=255)

                                                                  aCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCgg                       >  2:2279660/1‑93 (MQ=255)

                                                                  aCATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCgg                       >  1:5895352/1‑93 (MQ=255)

                                                                   cATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGc                      >  2:316662/1‑93 (MQ=255)

                                                                   cATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGGGGTCGGAATCCTATGCAAAACAGGTCAGATCGGc                      >  2:3472815/1‑93 (MQ=255)

                                                                    aTGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCa                     <  2:595043/93‑1 (MQ=255)

                                                                     tGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCAc                    <  1:3272540/93‑1 (MQ=255)

                                                                       gcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTa                  <  1:751971/93‑1 (MQ=255)

                                                                         ggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGa                >  1:4046799/1‑93 (MQ=255)

                                                                           cAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaa              >  1:2471699/1‑93 (MQ=255)

                                                                            aGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaa             >  1:2382390/1‑93 (MQ=255)

                                                                            aGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaa             >  2:2017125/1‑93 (MQ=255)

                                                                            aGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGGTGGAGGGGGGTCGGAATCCTATGCAAAACAGGGCAGATCGGCACTAGaaaa             >  1:4378855/1‑93 (MQ=255)

                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            >  2:3862199/1‑93 (MQ=255)

                                                                              gATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAg           >  2:1201248/1‑93 (MQ=255)

                                                                               aTTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGa          <  2:5811087/93‑1 (MQ=255)

                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         >  2:5749909/1‑93 (MQ=255)

                                                                                  tAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGt       <  2:1446626/93‑1 (MQ=255)

                                                                                   aagaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTc      >  2:170306/1‑93 (MQ=255)

                                                                                    agaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     >  2:6338961/1‑93 (MQ=255)

                                                                                    agaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     >  2:2604566/1‑93 (MQ=255)

                                                                                    agaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     >  1:3519625/1‑93 (MQ=255)

                                                                                     gaagGCGGCGGAAGACTGGTGTGCTTCCAGGGTTGGAGGGTGGTCGGAATCATATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTa    <  1:182979/93‑1 (MQ=255)

CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289

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