BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,587,936	A→G	R30R (CGT→CGC)	b1505 ←	putative outer membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,587,936	0	A	G	100.0%	807.6 / NA	245	R30R (CGT→CGC)	b1505	putative outer membrane protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (107/138); total (107/138)

Only 100 of 251 total aligned reads displayed.

TGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAT  >  NC_000913/1587846‑1588025

tGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACcgcg                                                                                         >  2:2702756/1‑93 (MQ=255)

 gACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCgcgc                                                                                        <  2:3837255/93‑1 (MQ=255)

  aCGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCt                                                                                       <  1:300216/93‑1 (MQ=255)

    gtgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTa                                                                                     >  2:2759831/1‑93 (MQ=255)

     tgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTAc                                                                                    >  2:786495/1‑93 (MQ=255)

     tgtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTAc                                                                                    <  1:6108131/93‑1 (MQ=255)

      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   <  1:4594621/93‑1 (MQ=255)

      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   <  1:3620443/93‑1 (MQ=255)

      gtTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACt                                                                                   <  2:5190213/93‑1 (MQ=255)

        tCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGt                                                                                 <  2:1776713/93‑1 (MQ=255)

         ccAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTa                                                                                <  1:6702313/93‑1 (MQ=255)

           aGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAGTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGa                                                                              <  2:1938625/93‑1 (MQ=255)

            gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGac                                                                             <  2:2198120/93‑2 (MQ=255)

            gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                             >  2:5098761/1‑93 (MQ=255)

            gTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGaa                                                                             >  2:4203359/1‑93 (MQ=255)

              accTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACa                                                                           <  2:3440968/91‑1 (MQ=255)

               aCTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAg                                                                          >  1:4789971/1‑93 (MQ=255)

                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         >  1:882717/1‑93 (MQ=255)

                cTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGa                                                                         <  1:5003550/93‑1 (MQ=255)

                  ttGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGAtt                                                                       <  2:4989842/93‑1 (MQ=255)

                   tGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTa                                                                      >  2:1736529/1‑93 (MQ=255)

                   tGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTa                                                                      >  1:2098116/1‑93 (MQ=255)

                    ggCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaa                                                                     >  2:1541104/1‑93 (MQ=255)

                    ggCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaa                                                                     >  2:482702/1‑93 (MQ=255)

                    ggCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaa                                                                     <  2:2416884/93‑1 (MQ=255)

                     gCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaaa                                                                    >  1:658983/1‑93 (MQ=255)

                     gCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaaa                                                                    >  1:4932060/1‑93 (MQ=255)

                     gCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTaaa                                                                    <  1:2612062/93‑1 (MQ=255)

                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  <  2:4516939/93‑1 (MQ=255)

                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  >  1:5961876/1‑93 (MQ=255)

                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  >  2:6005904/1‑93 (MQ=255)

                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  <  2:2593674/93‑1 (MQ=255)

                       gACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaat                                                                  <  1:2503179/93‑1 (MQ=255)

                        aCTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaata                                                                 >  1:3234342/1‑93 (MQ=255)

                         cTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataa                                                                <  1:2737166/93‑1 (MQ=255)

                          tGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataat                                                               <  2:3698125/93‑1 (MQ=255)

                           ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAaataata                                                              >  2:4577730/1‑92 (MQ=255)

                           ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATc                                                              >  1:6174175/1‑93 (MQ=255)

                           ggCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATc                                                              <  1:1499450/93‑1 (MQ=255)

                             cACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCaa                                                            <  1:4070766/93‑1 (MQ=255)

                              aCTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAAt                                                           <  1:3475266/93‑1 (MQ=255)

                              aCTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAAt                                                           <  2:3636447/93‑1 (MQ=255)

                              aCTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAAt                                                           >  1:896351/1‑93 (MQ=255)

                                 gAAAAATGTCGTACCGTAATTTCCAAGCTGCTGTGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATaaa                                                        >  2:432662/1‑93 (MQ=255)

                                  aaaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAt                                                       >  1:5886310/1‑93 (MQ=255)

                                  aaaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAt                                                       >  2:6027177/1‑93 (MQ=255)

                                  aaaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAt                                                       <  2:3976341/93‑1 (MQ=255)

                                  aaaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAt                                                       <  2:6614852/93‑1 (MQ=255)

                                   aaaaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAAtt                                                      <  2:2959604/93‑1 (MQ=255)

                                     aaTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgt                                                    >  1:3598958/1‑93 (MQ=255)

                                      aTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgtg                                                   <  1:4192072/93‑1 (MQ=255)

                                      aTGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATtgtg                                                   >  1:753294/1‑93 (MQ=255)

                                       tGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTgtgt                                                  >  2:5595490/1‑93 (MQ=255)

                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 <  2:354085/93‑1 (MQ=255)

                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 <  1:6253131/93‑1 (MQ=255)

                                        gTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTc                                                 <  2:4278661/93‑1 (MQ=255)

                                          cGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTg                                               <  1:4054591/93‑1 (MQ=255)

                                           gTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGc                                              <  2:5774696/93‑1 (MQ=255)

                                           gTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGc                                              <  2:6548216/93‑1 (MQ=255)

                                            tACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCt                                             <  2:588561/93‑1 (MQ=255)

                                            tACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCt                                             >  1:2197063/1‑93 (MQ=255)

                                              ccGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCa                                           <  1:3077417/93‑1 (MQ=255)

                                              ccGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCa                                           >  1:1698633/1‑93 (MQ=255)

                                               cGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAt                                          >  2:5860506/1‑93 (MQ=255)

                                               cGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAt                                          >  1:3841495/1‑93 (MQ=255)

                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         >  1:834834/1‑93 (MQ=255)

                                                gTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCAtt                                         >  1:1722411/1‑93 (MQ=255)

                                                  aaTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTgg                                       >  1:3185464/1‑93 (MQ=255)

                                                  aaTTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTgg                                       <  1:4445138/93‑1 (MQ=255)

                                                    tttCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGgtg                                     >  1:4290850/1‑93 (MQ=255)

                                                        cAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtctc                                 <  2:3066266/93‑1 (MQ=255)

                                                        cAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGtctc                                 <  2:1151531/93‑1 (MQ=255)

                                                          aGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCg                               >  2:2268159/1‑93 (MQ=255)

                                                          aGCTGCTGAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCg                               >  1:2776073/1‑93 (MQ=255)

                                                             tgctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTa                            <  1:6658412/93‑1 (MQ=255)

                                                              gctgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTag                           <  2:1470532/93‑1 (MQ=255)

                                                                tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         >  1:1341360/1‑93 (MQ=255)

                                                                tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         <  1:1878935/93‑1 (MQ=255)

                                                                tgAGAGATGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTagga                         >  2:4224026/1‑93 (MQ=255)

                                                                 gagagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTaggag                        <  2:128585/93‑1 (MQ=255)

                                                                   gagaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTATCTGCTCATTGGTGTCTCCGGTAGGAGGc                      >  1:2031879/1‑93 (MQ=255)

                                                                    agaTTCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCt                     <  1:6376912/93‑1 (MQ=255)

                                                                    agaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCt                     <  1:3273876/93‑1 (MQ=255)

                                                                     gaTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCtt                    >  1:4346893/1‑93 (MQ=255)

                                                                      aTGCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTg                   >  1:3444578/1‑93 (MQ=255)

                                                                        gCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGAc                 <  2:261528/93‑1 (MQ=255)

                                                                        gCTTATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGAc                 <  1:5100842/93‑1 (MQ=255)

                                                                          ttATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACgg               >  1:5218834/1‑93 (MQ=255)

                                                                           tATTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgt              >  2:5531555/1‑93 (MQ=255)

                                                                            aTTTGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGgtg             <  2:5728606/93‑1 (MQ=255)

                                                                             tttGTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGCGTCTCCGGTAGGAGGCTTGACGgcgt            >  2:3232362/1‑93 (MQ=255)

                                                                                gTTCCTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAAc         <  2:3167018/93‑1 (MQ=255)

                                                                                   ccTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCAc      <  1:1061209/93‑1 (MQ=255)

                                                                                   ccTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCAc      <  1:5668963/93‑1 (MQ=255)

                                                                                    cTGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACt     >  1:5528081/1‑93 (MQ=255)

                                                                                     tGACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTc    >  1:2369337/1‑93 (MQ=255)

                                                                                      gACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCa   >  2:6313361/1‑93 (MQ=255)

                                                                                      gACCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCa   <  1:157443/93‑1 (MQ=255)

                                                                                       aCCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAt  >  1:3711461/1‑93 (MQ=255)

                                                                                       aCCGCGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAt  <  2:6411328/93‑1 (MQ=255)

TGACGTGTTCCAGTAACTTTGGCGACTGGCACTGAAAAATGTCGTACCGTAATTTCCAAGCTGCTGAGAGATGCTTATTTGTTCCTGACCACGCTTACTGTAGAACAGATTAAAATAATCAATAAATTGTGTCTGCTCATTGGTGTCTCCGGTAGGAGGCTTGACGGTGTAACCACTCAT  >  NC_000913/1587846‑1588025

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