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%	944.7 / NA	206	intergenic (+43/‑25)	b4316/fimD	Pil/FimD
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (108/98); total (108/98)

Only 100 of 222 total aligned reads displayed.

AAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542549‑4542728

  aTGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCgggg                                                                                       >  2:1543624/1‑93 (MQ=255)

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

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

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

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

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

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

           gcgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAg                                                                             >  2:4475378/1‑93 (MQ=255)

           gcgcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAg                                                                             >  1:812183/1‑93 (MQ=255)

             gcACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGcc                                                                           <  2:1391506/93‑1 (MQ=255)

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

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

               aCTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgg                                                                         <  1:2397339/93‑1 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggg                                                                        >  2:1245076/1‑92 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggg                                                                        >  1:3635080/1‑92 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  2:395821/93‑1 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  2:2442316/93‑1 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  1:1088210/93‑1 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  1:2258471/93‑1 (MQ=255)

                cTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCgga                                                                        <  1:1143764/93‑1 (MQ=255)

                 ttACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggag                                                                       <  1:4220091/93‑1 (MQ=255)

                 ttACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggag                                                                       <  2:1577506/93‑1 (MQ=255)

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

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

                   aCCCCCAAAATGACGGGCGTAATTGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCggagga                                                                     >  2:3795715/1‑93 (MQ=255)

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

                      cccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaata                                                                  >  1:1725343/1‑93 (MQ=255)

                      cccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaata                                                                  <  1:4247249/93‑1 (MQ=255)

                      cccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaata                                                                  <  1:1455473/93‑1 (MQ=255)

                       ccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataa                                                                 >  2:3682266/1‑93 (MQ=255)

                       ccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGaataa                                                                 <  2:3568541/93‑1 (MQ=255)

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

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

                         aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATg                                                               <  2:4296543/93‑1 (MQ=255)

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

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

                          aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGt                                                              >  2:1095848/1‑93 (MQ=255)

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

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

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

                           aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTc                                                             <  1:3243542/93‑1 (MQ=255)

                            aTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCa                                                            >  1:3311962/1‑93 (MQ=255)

                             tGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCat                                                           >  2:3226920/1‑93 (MQ=255)

                              gACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCata                                                          <  2:4008229/93‑1 (MQ=255)

                                cGGGCGTAATGGAATAACGCAGGGGGAAGTTGTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATc                                                        <  1:2567211/93‑1 (MQ=255)

                                    cGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGaa                                                    <  1:316167/93‑1 (MQ=255)

                                    cGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGaa                                                    <  2:1042386/93‑1 (MQ=255)

                                    cGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGaa                                                    <  2:2591473/93‑1 (MQ=255)

                                         tGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTaa                                               >  1:2213782/1‑93 (MQ=255)

                                          ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAg                                              <  2:2754963/93‑1 (MQ=255)

                                          ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAg                                              <  2:530407/93‑1 (MQ=255)

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

                                             aTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACt                                           <  1:3745537/93‑1 (MQ=255)

                                              tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACtt                                          >  1:2311347/1‑93 (MQ=255)

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

                                                aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTa                                        >  2:3431230/1‑93 (MQ=255)

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

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

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

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

                                                      ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCg                                  >  2:386006/1‑93 (MQ=255)

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

                                                       gggAATTTTTCGCCTCAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGa                                 <  2:4500405/93‑1 (MQ=255)

                                                        ggAATTTTTCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGaa                                <  2:2795273/93‑1 (MQ=255)

                                                        ggAATATTTCGCATGAATAAAAAGAAGTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAGTGTCATATCTGAATTTAAGACTTTACCAGCGaa                                <  2:1600602/93‑1 (MQ=255)

                                                          aaTTTTTCGCTTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAac                              <  2:1777039/93‑1 (MQ=255)

                                                                tCGCCTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACaa                        <  1:513279/93‑1 (MQ=255)

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

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

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

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

                                                                    cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGATACACACAATGCt                    >  1:1279261/1‑93 (MQ=255)

                                                                    cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    >  1:4144971/1‑93 (MQ=255)

                                                                    cTGAATAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCt                    >  2:361190/1‑93 (MQ=255)

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

                                                                       aaTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGc                 <  1:1484655/93‑1 (MQ=255)

                                                                       aaTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGc                 >  1:1828129/1‑93 (MQ=255)

                                                                       aaTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGc                 >  2:3319857/1‑93 (MQ=255)

                                                                        aTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCa                <  1:873087/93‑1 (MQ=255)

                                                                        aTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCa                >  2:414642/1‑93 (MQ=255)

                                                                        aTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCa                <  2:4402157/93‑1 (MQ=255)

                                                                        aTAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACAAGCGAAACACACAATGCTTGCa                <  1:2728051/93‑1 (MQ=255)

                                                                         tAAAAAGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCat               >  2:3192161/1‑93 (MQ=255)

                                                                          aaaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCata              >  2:1010158/1‑93 (MQ=255)

                                                                           aaaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCatat             <  2:1950281/93‑1 (MQ=255)

                                                                            aaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATAtt            >  2:1938613/1‑93 (MQ=255)

                                                                            aaaGAATTGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCGTGCATAtt            <  2:4568910/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                                   tGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGc     <  2:3073062/93‑1 (MQ=255)

                                                                                   tGACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGc     <  2:2529027/93‑1 (MQ=255)

                                                                                    gACTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCa    <  2:3895497/93‑1 (MQ=255)

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

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

                                                                                     aCTGCC‑GGGTGATTTTAAGCCGGACGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCAt   <  1:2616831/93‑1 (MQ=255)

                                                                                      cTGCC‑GGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATc  <  2:1852524/93‑1 (MQ=255)

AAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATC  >  NC_000913/4542549‑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: `‑`