BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,542,643	Δ1 bp	intergenic (+46/‑22)	b4316 → / → fimD	Pil/FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,640	0	G	.	100.0%	941.9 / NA	203	intergenic (+43/‑25)	b4316/fimD	Pil/FimD
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (104/99); total (104/99)

Only 100 of 217 total aligned reads displayed.

TAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542548‑4542728

tAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCg                                                                                          >  2:3600637/1‑93 (MQ=255)

tAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCg                                                                                          >  2:2177067/1‑93 (MQ=255)

 aaaTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgg                                                                                         <  2:5449035/93‑1 (MQ=255)

  aaTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCATGAATAAAAAGAATTGACTGCCggg                                                                                        <  1:1504800/93‑1 (MQ=255)

    tGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggtg                                                                                      >  1:870427/1‑91 (MQ=255)

       ttATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtt                                                                                  >  1:6209047/1‑93 (MQ=255)

        tATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGTGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAttt                                                                                 <  2:1491814/93‑1 (MQ=255)

        tATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAttt                                                                                 >  1:5998339/1‑93 (MQ=255)

        tATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAttt                                                                                 <  2:5733979/93‑1 (MQ=255)

         aTGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtttt                                                                                <  1:4902028/93‑1 (MQ=255)

         aTGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGAtttt                                                                                >  1:1981406/1‑93 (MQ=255)

           ggCGCACTTACCCCAAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTaa                                                                              <  2:2816485/93‑1 (MQ=255)

             cgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGc                                                                            <  1:3158480/93‑1 (MQ=255)

              gcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                                           >  1:4467846/1‑93 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          >  2:3255132/1‑93 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          >  1:1378500/1‑93 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          >  1:2677117/1‑93 (MQ=255)

               cACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCg                                                                          <  2:5357661/93‑1 (MQ=255)

                aCTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                                         >  1:4109749/1‑93 (MQ=255)

                aCTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                                         >  2:5394006/1‑93 (MQ=255)

                 cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        >  2:1949048/1‑93 (MQ=255)

                 cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        >  2:3723539/1‑93 (MQ=255)

                   tACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagg                                                                      <  2:357860/93‑1 (MQ=255)

                     cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaa                                                                    >  1:516424/1‑93 (MQ=255)

                     cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaa                                                                    >  1:4716274/1‑93 (MQ=255)

                     cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaa                                                                    >  1:1206529/1‑93 (MQ=255)

                      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                                   <  1:209656/93‑1 (MQ=255)

                      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                                   >  2:2472566/1‑93 (MQ=255)

                      ccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaat                                                                   >  1:5506444/1‑93 (MQ=255)

                         cAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataat                                                                >  1:6290984/1‑93 (MQ=255)

                         aaaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataat                                                                <  2:1649739/92‑1 (MQ=255)

                          aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                                               >  1:4277314/1‑93 (MQ=255)

                          aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                                               >  1:5245624/1‑93 (MQ=255)

                          aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGACGGAGGAATAATg                                                               >  2:4585410/1‑93 (MQ=255)

                            aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTc                                                             >  1:4011891/1‑93 (MQ=255)

                             aTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCa                                                            <  2:4439571/93‑1 (MQ=255)

                                aCGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCatat                                                         <  1:2425555/93‑1 (MQ=255)

                                 cGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGGCTGCC‑GGGTGATTTTAAGCCGGAGGAATAACGTCATATc                                                        <  2:6161123/93‑1 (MQ=255)

                                      gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAt                                                   <  1:1792018/93‑1 (MQ=255)

                                       tAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAtt                                                  >  1:3266662/1‑93 (MQ=255)

                                        aaTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAAttt                                                 >  1:5895559/1‑93 (MQ=255)

                                         aTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTa                                                <  1:1236478/93‑1 (MQ=255)

                                            gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGa                                             <  1:1884806/93‑1 (MQ=255)

                                            gAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGa                                             >  2:3449331/1‑93 (MQ=255)

                                             aaTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGAc                                            <  1:2728545/93‑1 (MQ=255)

                                              aTAACGCAGGGGGAATTTTTCGCCTGAATAAAAATAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACt                                           <  2:2283030/93‑1 (MQ=255)

                                              aTAACACAGGGGGAATTTTTCGCTTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACt                                           <  2:4183258/93‑1 (MQ=255)

                                               tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACtt                                          <  1:3439660/93‑1 (MQ=255)

                                                aaCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                                         <  2:6209047/93‑1 (MQ=255)

                                                aaCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                                         <  1:1321191/93‑1 (MQ=255)

                                                aaCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                                         >  1:801190/1‑93 (MQ=255)

                                                 aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa                                        <  2:264944/93‑1 (MQ=255)

                                                  cGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAc                                       <  2:1525051/93‑1 (MQ=255)

                                                   gCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                                      <  2:1199006/93‑1 (MQ=255)

                                                   gCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                                      >  2:1775698/1‑93 (MQ=255)

                                                   gCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTCTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                                      >  2:1247954/1‑93 (MQ=255)

                                                    cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCa                                     <  2:4272630/93‑1 (MQ=255)

                                                      gggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGc                                   >  2:2698672/1‑93 (MQ=255)

                                                        gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGa                                 <  1:6021660/93‑1 (MQ=255)

                                                        gggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGa                                 >  2:5992597/1‑93 (MQ=255)

                                                          gAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaaa                               >  1:4133419/1‑93 (MQ=255)

                                                            aTTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca                             <  2:2923122/93‑1 (MQ=255)

                                                             tttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacac                            >  1:4770778/1‑93 (MQ=255)

                                                             tttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacac                            <  1:5038932/93‑1 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           >  2:557187/1‑93 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           <  2:4667711/93‑1 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           >  2:4603035/1‑93 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           <  1:2208489/93‑1 (MQ=255)

                                                              ttttCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca                           <  2:355397/93‑1 (MQ=255)

                                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACaa                        >  1:3907475/1‑93 (MQ=255)

                                                                 tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACaa                        >  2:2834404/1‑93 (MQ=255)

                                                                  cGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAAt                       <  2:1452661/93‑1 (MQ=255)

                                                                   gCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg                      >  2:5481975/1‑93 (MQ=255)

                                                                   gCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg                      >  2:3844578/1‑93 (MQ=255)

                                                                   gCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg                      >  1:4198794/1‑93 (MQ=255)

                                                                   gCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg                      >  2:6026059/1‑93 (MQ=255)

                                                                    ccTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGc                     >  2:5290020/1‑93 (MQ=255)

                                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    <  1:1035465/93‑1 (MQ=255)

                                                                     cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    <  1:3654415/93‑1 (MQ=255)

                                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt                   <  2:5638814/93‑1 (MQ=255)

                                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt                   >  2:23541/1‑93 (MQ=255)

                                                                      tGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCtt                   >  1:3025178/1‑93 (MQ=255)

                                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata              <  1:4851998/93‑1 (MQ=255)

                                                                           aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata              <  2:5846038/93‑1 (MQ=255)

                                                                            aaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCatat             >  2:3947965/1‑93 (MQ=255)

                                                                              aaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTc           <  1:4631153/93‑1 (MQ=255)

                                                                               aGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCg          >  1:1704723/1‑93 (MQ=255)

                                                                               aGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCg          >  1:6017275/1‑93 (MQ=255)

                                                                                gAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGt         >  2:1071937/1‑93 (MQ=255)

                                                                                gAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGt         >  2:5344140/1‑93 (MQ=255)

                                                                                gAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGt         >  1:1542632/1‑93 (MQ=255)

                                                                                 aaTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTa        >  1:4415644/1‑93 (MQ=255)

                                                                                  aTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTaa       <  1:4970959/93‑1 (MQ=255)

                                                                                  aTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTaa       >  1:1749733/1‑93 (MQ=255)

                                                                                  aTTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTaa       <  1:1626057/93‑1 (MQ=255)

                                                                                  aTTGACTGCC‑GGGCGATTTTAAGCCGGAGGAATAATGTCATGTCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTaa       >  2:3323974/1‑93 (MQ=255)

                                                                                   ttGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAg      <  2:503590/93‑1 (MQ=255)

                                                                                     gACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCa    >  2:787970/1‑93 (MQ=255)

                                                                                      aCTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCAt   >  2:5315202/1‑93 (MQ=255)

                                                                                      aCTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCAt   <  1:4268067/93‑1 (MQ=255)

TAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542548‑4542728

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