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%	306.8 / NA	114	I502I (ATT→ATC)	fimD	FimD
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (69/45); total (69/45)

Only 100 of 114 total aligned reads displayed.

CGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACG  >  NC_000913/4544102‑4544241

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

cGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAacac                                                                  <  1:2136598/76‑1 (MQ=255)

 gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAAcaca                                                                 <  1:1575033/76‑1 (MQ=255)

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

 gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAAcaca                                                                 <  1:2987965/76‑1 (MQ=255)

 gTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAAcaca                                                                 <  2:3355044/76‑1 (MQ=255)

  ttcttCGAACAGCTGATCTTTTAAGTTCGCTGATACAACATACAGTCGACTGAATGGCTACAACATCGAAACACAg                                                                <  1:3147072/73‑1 (MQ=255)

  ttATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAg                                                                <  2:1858236/76‑1 (MQ=255)

    aTTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGa                                                              >  2:3643111/1‑76 (MQ=255)

      tCGACCAGCGGCTATTTTAATTTCGCTGATACACCACACACTCGCATGAATGGCTACAACATCGACTCACAGGACg                                                            <  1:2116615/76‑1 (MQ=255)

      tCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                            >  1:4758504/1‑76 (MQ=255)

      tCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACg                                                            >  1:3955854/1‑76 (MQ=255)

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

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

           cAGCGGATATTTAAATTTCGCCGAGACAACATACAGCCGAATGAATAGAGACAACAGCGAAACACACGGCGGAggt                                                       >  2:3479069/1‑74 (MQ=255)

            accGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTa                                                      >  2:2176815/3‑76 (MQ=255)

            aGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTa                                                      >  1:4843906/1‑76 (MQ=255)

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

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

              cGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTAtt                                                    >  2:572418/1‑76 (MQ=255)

               ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTc                                                   <  2:3287258/76‑1 (MQ=255)

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

               ggATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACACCGAACCCCAGGCCGGATTTTTTc                                                   >  2:319732/1‑76 (MQ=255)

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

                  tatTTTAATTTCGTTGGTACAACATACAGTCGAATGAATGGCTGCAACATCGAAACACAGGACGGAGTTATTCAga                                                <  2:4486915/76‑2 (MQ=255)

                  tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACCACATCGAAACACAGGACGGCGTTATTCAgg                                                >  2:1266819/1‑76 (MQ=255)

                  tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAgg                                                >  2:1236320/1‑76 (MQ=255)

                  tatTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAgg                                                <  1:2990110/76‑1 (MQ=255)

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

                    ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                              <  2:1139927/76‑1 (MQ=255)

                    ttttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGtt                                              >  2:4058381/1‑76 (MQ=255)

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

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

                     tttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGATa                                             >  2:655414/1‑76 (MQ=255)

                      ttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                            <  1:1963777/76‑1 (MQ=255)

                      ttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                            <  1:4706963/76‑1 (MQ=255)

                      ttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                            >  1:4805795/1‑76 (MQ=255)

                      ttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                            >  2:1382491/1‑76 (MQ=255)

                      ttAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTaa                                            >  1:2290737/1‑76 (MQ=255)

                       tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAg                                           <  1:3338414/76‑1 (MQ=255)

                       tAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACCGAGTTATTCAGGTTAAg                                           >  1:4462956/1‑76 (MQ=255)

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

                         aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                         >  1:3116726/1‑76 (MQ=255)

                         aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                         >  1:4139470/1‑76 (MQ=255)

                         aTTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGcc                                         <  1:3792202/76‑1 (MQ=255)

                          tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                        >  1:1905908/1‑76 (MQ=255)

                          tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                        >  1:1852833/1‑76 (MQ=255)

                          tttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCg                                        >  1:171864/1‑76 (MQ=255)

                           ttCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGa                                       >  1:4020236/1‑76 (MQ=255)

                            tCGCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAATTATTCAGGTTAAGCCGaa                                      >  2:3630174/1‑76 (MQ=255)

                              gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAAt                                    >  1:3155351/1‑76 (MQ=255)

                              gCTGATACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGGTATTCAGGTTAAGCCGAAAt                                    >  1:838803/1‑76 (MQ=255)

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

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

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

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

                               cTGATACAACATACAGTCGAATGAATGGCTAAAACATCGAAAAACAGGACGGAGTTATTCAGGTTAAGACGAAAtt                                   >  1:1601636/1‑76 (MQ=255)

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

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

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

                                  aTACAACATACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTCTTCAGGTTAAGCCGAAATTCAc                                >  1:4642170/1‑76 (MQ=255)

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

                                    acaacaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCg                              <  1:3875765/76‑1 (MQ=255)

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

                                       acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                           <  2:749868/76‑1 (MQ=255)

                                       acaTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACt                           <  2:2709882/76‑1 (MQ=255)

                                         aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAt                         <  1:1470286/76‑1 (MQ=255)

                                         aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAt                         <  1:4689094/76‑1 (MQ=255)

                                         aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAt                         >  1:2537622/1‑76 (MQ=255)

                                         aTACAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAt                         <  1:4606956/76‑1 (MQ=255)

                                           aCAGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTa                       >  2:1136991/1‑76 (MQ=255)

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

                                             aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACa                     >  1:803204/1‑76 (MQ=255)

                                             aGTCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACa                     >  2:2575621/1‑76 (MQ=255)

                                               tCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                   >  2:2586984/1‑76 (MQ=255)

                                               tCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                   >  2:4777317/1‑76 (MQ=255)

                                               tCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                   <  1:4184691/76‑1 (MQ=255)

                                               tCGAATGAATGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAc                   <  1:3542143/76‑1 (MQ=255)

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

                                                cGAATGAATGGCTACAACATCGAAACACAGGACGGAGCTATTCAGGTTAAGCCGAAATTCACCGACTATTACAAcc                  >  2:731713/1‑76 (MQ=255)

                                                 gaatgaatGGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCt                 >  1:4684898/1‑76 (MQ=255)

                                                  aatgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                <  2:1158859/76‑1 (MQ=255)

                                                  aatgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTc                >  1:879789/1‑76 (MQ=255)

                                                   atgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTATAACCTCg               >  1:1363181/1‑76 (MQ=255)

                                                   atgaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCg               >  1:1326567/1‑76 (MQ=255)

                                                     gaatgGCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCt             >  2:3690386/1‑76 (MQ=255)

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

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

                                                        tgGCTACAACATCTAAACGCAGGCCGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTAGTAGAACCTCGCTtat          <  2:4257163/76‑1 (MQ=255)

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

                                                         ggCTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTtata         <  1:270961/76‑1 (MQ=255)

                                                         agCTACAACATCGAAACACAGGACGGAGTTATTCATTTTAAGCCGCAATTCACCGACTATTACAACCTCGCTtata         <  1:24828/75‑1 (MQ=255)

                                                          gCTACAACATCGAAACACAGGACGGAGTTATTAAGGTTACGCCGAAATTCCCCGCCTATTACACACTCGCTTCTaa        >  2:1359698/1‑76 (MQ=255)

                                                           cTACAACATCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAAc       >  2:3434274/1‑76 (MQ=255)

                                                            tACAACATCGAAACACAGGACGGAGTTATTCAGGCTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATTACa      >  1:1916799/1‑76 (MQ=255)

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

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

                                                               aacaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAc   >  1:4828027/1‑76 (MQ=255)

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

                                                                acaTCGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAAcg  >  2:278901/1‑76 (MQ=255)

CGTTATTCGACCAGCGGATATTTTAATTTCGCTGATACAACATACAGTCGAATGAATGGCTACAACATTGAAACACAGGACGGAGTTATTCAGGTTAAGCCGAAATTCACCGACTATTACAACCTCGCTTATAACAAACG  >  NC_000913/4544102‑4544241

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