BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,544,170	T→C	I502I (ATT→ATC)	fimD →	FimD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,544,170	0	T	C	100.0%	423.8 / NA	150	I502I (ATT→ATC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (89/61); total (89/61)

Only 100 of 150 total aligned reads displayed.

GTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAC  >  NC_000913/4544097‑4544240

gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACCACATCGa                                                                      >  2:3764960/1‑76 (MQ=255)

gTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGa                                                                      >  1:5574026/1‑76 (MQ=255)

 ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaa                                                                     >  1:1216153/1‑76 (MQ=255)

 ttACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaa                                                                     <  2:2767065/76‑1 (MQ=255)

  tACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGaaa                                                                    >  2:2559632/1‑76 (MQ=255)

    ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAaca                                                                  <  1:2174171/76‑1 (MQ=255)

    ccGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACAGCGAAaca                                                                  >  2:4357170/1‑76 (MQ=255)

    ccGTTATGTGAACAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAaca                                                                  <  2:4465742/76‑1 (MQ=255)

     cGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAacac                                                                 >  2:4098671/1‑76 (MQ=255)

     cGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAacac                                                                 >  1:808737/1‑76 (MQ=255)

      gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAAcaca                                                                >  2:1222992/1‑76 (MQ=255)

        tATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAgg                                                              <  1:3020072/76‑1 (MQ=255)

         aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                             <  1:154488/76‑1 (MQ=255)

         aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                             <  2:2380615/76‑1 (MQ=255)

             gACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGa                                                         <  2:1727452/76‑1 (MQ=255)

              aCCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAg                                                        >  2:2172314/1‑76 (MQ=255)

               ccAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGt                                                       >  1:5093897/1‑76 (MQ=255)

                cAGCGGATATTTTAATTTCGCTGATACACCATACAGTCGAATGAATGGCTACAACATCGAAACCAAGGACGGAGtt                                                      >  1:5796681/1‑76 (MQ=255)

                cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                      <  2:4554473/76‑1 (MQ=255)

                cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                      <  2:5727570/76‑1 (MQ=255)

                cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                      >  1:5783794/1‑76 (MQ=255)

                cAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGtt                                                      >  1:1231687/1‑76 (MQ=255)

                 aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTa                                                     <  1:2582264/76‑1 (MQ=255)

                 aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAAACTCGAAACACCGGCCGGAGTTa                                                     >  2:271974/1‑76 (MQ=255)

                  gCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                    >  1:2907070/1‑76 (MQ=255)

                  gCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAt                                                    >  1:3030681/1‑76 (MQ=255)

                   cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                   <  2:5624059/76‑1 (MQ=255)

                   cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                   <  2:5135054/76‑1 (MQ=255)

                    ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                  >  1:1784831/1‑76 (MQ=255)

                    ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                  >  1:5393292/1‑76 (MQ=255)

                     gATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCa                                                 >  1:532492/1‑76 (MQ=255)

                      atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAg                                                <  2:1827402/76‑1 (MQ=255)

                      atatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAg                                                >  1:4131550/1‑76 (MQ=255)

                        atTTTAATTTCGCTGCTACAACATACAGTCGCATGAATCGCTACCCACTCGAACCACAGGACGGAGTTATTCAGGt                                              >  2:945165/1‑76 (MQ=255)

                        atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                              >  2:935891/1‑76 (MQ=255)

                        atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                              <  1:3275300/76‑1 (MQ=255)

                        atTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGt                                              >  1:916526/1‑76 (MQ=255)

                         ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                             <  2:374375/60‑1 (MQ=255)

                         ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                             >  1:374375/1‑60 (MQ=255)

                          tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                            >  2:4377106/1‑76 (MQ=255)

                          tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                            <  2:5726874/76‑1 (MQ=255)

                          tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTa                                            >  2:4922956/1‑76 (MQ=255)

                            tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAg                                          >  1:1205207/1‑76 (MQ=255)

                             aaTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGc                                         >  2:420022/1‑76 (MQ=255)

                             aaTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGc                                         <  1:5294533/76‑1 (MQ=255)

                                ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCTAAACACAGGACGGAGTTATTCAGGTTAAGCCGa                                      <  1:2119068/76‑1 (MQ=255)

                                ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGCCGGAGTTATTCCGGTTAACCCGa                                      >  2:561825/1‑76 (MQ=255)

                                 tCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaa                                     >  1:1523806/1‑76 (MQ=255)

                                 tCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGGTAAGCCGaa                                     >  1:962882/1‑76 (MQ=255)

                                  cGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGaaa                                    >  1:2598897/1‑76 (MQ=255)

                                    cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                  <  2:3497775/76‑1 (MQ=255)

                                    cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                  <  2:3527064/76‑1 (MQ=255)

                                    cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                  >  2:4657690/1‑76 (MQ=255)

                                    cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                  >  2:4908075/1‑76 (MQ=255)

                                    cTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAtt                                  >  1:3878245/1‑76 (MQ=255)

                                    cTGATACAACATACAGTCGAATGAATGGCTACAACAGCGAAACACAACACGGAGTTACTCAGGTTAAGCCGAAAtt                                  >  2:4270717/1‑76 (MQ=255)

                                     tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTc                                 >  1:2705717/1‑76 (MQ=255)

                                     tGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATGCAGGTTAAGCCGAAATTc                                 >  1:3517651/1‑76 (MQ=255)

                                      gATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCa                                >  2:1386746/1‑76 (MQ=255)

                                       aTACAACATGCAGTCGAATGAATGGCTACAACATCGAAAGACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAc                               <  2:5285113/76‑1 (MQ=255)

                                       aTACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAc                               >  1:1145456/1‑76 (MQ=255)

                                        tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                              <  2:4131549/76‑1 (MQ=255)

                                        tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                              >  1:4971198/1‑76 (MQ=255)

                                        tACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCAcc                              >  2:2916134/1‑76 (MQ=255)

                                         acaacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                             >  2:4271385/1‑76 (MQ=255)

                                          caacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGa                            >  2:3008752/1‑76 (MQ=255)

                                          caacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGa                            <  1:2287293/76‑1 (MQ=255)

                                           aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                           >  1:482577/1‑76 (MQ=255)

                                           aacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGAc                           >  2:4169669/1‑76 (MQ=255)

                                             caTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTa                         >  2:5497248/1‑76 (MQ=255)

                                               tACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAtt                       >  2:4835728/1‑76 (MQ=255)

                                                aCAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTa                      <  2:735893/76‑1 (MQ=255)

                                                 cAGTCGAATGAATGGCTACAACATCGAAACACGGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                     <  2:5116902/76‑1 (MQ=255)

                                                 cAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                     >  1:2423563/1‑76 (MQ=255)

                                                 cAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                     <  2:2646765/76‑1 (MQ=255)

                                                 cAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTAc                     <  1:654436/76‑1 (MQ=255)

                                                   ggcGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                   <  2:5370484/74‑1 (MQ=255)

                                                   gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                   >  1:243/1‑76 (MQ=255)

                                                   gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                   <  1:506104/76‑1 (MQ=255)

                                                   gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                   >  1:4709907/1‑76 (MQ=255)

                                                   gTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACaa                   >  1:3659961/1‑76 (MQ=255)

                                                    tCGAAGGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                  <  2:5483222/76‑1 (MQ=255)

                                                     cGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAcc                 >  1:3648740/1‑76 (MQ=255)

                                                      gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                <  1:3514412/76‑1 (MQ=255)

                                                        atgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCg              <  1:4000080/76‑1 (MQ=255)

                                                         tgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGc             >  1:2281791/1‑76 (MQ=255)

                                                          gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCt            <  2:2117541/76‑1 (MQ=255)

                                                           aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt           >  1:1991154/1‑76 (MQ=255)

                                                           aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt           <  1:5474748/76‑1 (MQ=255)

                                                           aatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCtt           <  1:4669764/76‑1 (MQ=255)

                                                            atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta          <  2:2907070/76‑1 (MQ=255)

                                                            atgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTta          <  2:779974/76‑1 (MQ=255)

                                                            atgGCTAAAACATCGAAACACAGGACGGCGCTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCGCGCtca          >  2:2667164/1‑74 (MQ=255)

                                                             tgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtat         >  1:3231249/1‑76 (MQ=255)

                                                                 tACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACa     >  1:1755427/1‑76 (MQ=255)

                                                                  acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaa    <  2:3187139/76‑1 (MQ=255)

                                                                  acaacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaa    >  2:3199089/1‑76 (MQ=255)

                                                                   caacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa   >  2:3235033/1‑76 (MQ=255)

                                                                   caacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACaaa   <  1:1989522/76‑1 (MQ=255)

                                                                    aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc  <  1:5520076/76‑1 (MQ=255)

GTTACCGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAC  >  NC_000913/4544097‑4544240

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`