BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,910	T→A	W30R (TGG→AGG)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,910	0	T	A	100.0%	708.1 / NA	188	W30R (TGG→AGG)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base A (104/84); total (104/84)

Only 100 of 253 total aligned reads displayed.

CGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCA  >  NC_000913/49848‑49969

cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAgg                                                           >  1:3287633/1‑65 (MQ=255)

cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAgg                                                           >  1:2409655/1‑65 (MQ=255)

cGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAgg                                                           >  1:1248658/1‑65 (MQ=255)

 ggTGGACCGTGTTTTCGGCTTGGAAATCGCAGTGCCGGGGAACCGGCCTGCCGATCGCGCCAGGt                                                          <  1:2505953/65‑1 (MQ=255)

 ggTAGATCGCGTTATCGGCATGGAAACCGCCATCCCGTGGAACCTGCCTGCCGATCGCGCCAGGt                                                          <  1:837173/65‑1 (MQ=255)

  gTAGTTCGCGTTATCGGCATGGAAACCGCCATGCGGTGGAACCTGCCTGCCGATCGCGCCAGGtt                                                         <  1:2388680/65‑1 (MQ=255)

  gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAggct                                                         >  1:3518377/1‑63 (MQ=255)

  gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGtt                                                         >  1:3827428/1‑65 (MQ=255)

     gATCGCGTTAGCGGCATGGAAACCGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTaa                                                      <  1:318036/65‑1 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGCTCGCGCCAGGTTTaaa                                                     >  1:163914/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaaa                                                     >  1:3492682/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaaa                                                     >  1:1050876/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaaa                                                     >  1:322312/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaaa                                                     >  1:601218/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaaa                                                     >  1:1907026/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaaa                                                     >  1:3683874/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTGTaaa                                                     >  1:2347046/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGGTTaaa                                                     >  1:3155574/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGATTTaaa                                                     >  1:1167225/1‑65 (MQ=255)

      aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTCCCGATCGCGCCAGGTTTaaa                                                     >  1:547936/1‑65 (MQ=255)

       tcacgTTATCGGCATGGAAAACCCCATCCCGGGGACCCTGCCTGCCGATCGCGCCAGGTTTAAAc                                                    <  1:140260/62‑1 (MQ=255)

         gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACgg                                                  >  1:2843870/1‑64 (MQ=255)

         gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGACGGCGCCAGGCTTAAACGc                                                  >  1:207969/1‑65 (MQ=255)

         gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCCAGCGCGCCAGGCTTAAACgt                                                  >  1:1284512/1‑64 (MQ=255)

         gcgTTATCGGCATGGAAAACGCCATGCCGGGGAACACGCCTGCCGATCGCGCCAGGTTTAAAAGc                                                  >  1:2200453/1‑65 (MQ=255)

              atgggCATGGAAAACGCCATGCCGGGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:607673/62‑1 (MQ=255)

              accGGCATGGAAACCGCCATGCCGGGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:475573/63‑1 (MQ=255)

              aTCGGCATGGAAAGCGCATTGCCTTGAACCCTCCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:2774754/65‑1 (MQ=255)

              aTCGGCATGGAAACCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:401157/65‑1 (MQ=255)

              aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCCGCCGATCGCGCCCGGGTTAAACGCAacac                                             >  1:2463156/1‑65 (MQ=255)

              aTCGGCATGGAAAACCCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:679388/65‑1 (MQ=255)

              aTCGGCATGAAAAGCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:699960/65‑1 (MQ=255)

              aTCGGCAGGGAAACCGCCATGCGGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:2609523/65‑1 (MQ=255)

              aTCGCCATGGAAACCGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                             <  1:1557471/65‑1 (MQ=255)

               tCGGCATGGAAACCGCCATGCGGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:74699/65‑1 (MQ=255)

               tCGGCATGGAAACCGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:3882716/65‑1 (MQ=255)

               tCGGCATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:412153/65‑1 (MQ=255)

               tCGGCATGAAAACCGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:1688619/65‑1 (MQ=255)

               tCGGCATCGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:3066693/65‑1 (MQ=255)

               gcGGAATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:1679609/64‑1 (MQ=255)

               ccGGCGTGGAAACCGCCATGCGGTGACCCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACAcc                                            <  1:3455623/64‑1 (MQ=255)

                 ggCATGGAAAAGGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCtt                                          >  1:2750208/1‑65 (MQ=255)

                 ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCtt                                          >  1:2142343/1‑65 (MQ=255)

                 ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCtt                                          >  1:1099286/1‑65 (MQ=255)

                 ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCtc                                          >  1:2037456/1‑64 (MQ=255)

                 ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGGTTAAACGCAACACCtt                                          >  1:2449454/1‑65 (MQ=255)

                   cagggAAAACGCCACGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                        <  1:286877/62‑1 (MQ=37)

                   cATGGAAACCGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                        <  1:2193540/65‑1 (MQ=255)

                   cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                        <  1:1057616/65‑1 (MQ=255)

                   cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                        <  1:1181251/65‑1 (MQ=255)

                   cATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                        <  1:2766373/65‑1 (MQ=255)

                   aaTGGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                        <  1:3006272/64‑1 (MQ=255)

                     tGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaact                                      >  1:3605956/1‑63 (MQ=255)

                        aaaaCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaaataaa                                   >  1:1359394/1‑65 (MQ=255)

                         aaaCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAAc                                  >  1:1480056/1‑65 (MQ=255)

                          aaCGCCAGCCCGTGAAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAAcc                                 <  1:2100594/65‑1 (MQ=255)

                           aCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAAccc                                >  1:2587636/1‑65 (MQ=255)

                            cGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCg                               <  1:2756225/65‑1 (MQ=255)

                             gCCATGCCGTGGAACCTGCTTGCCCATCGCGCCAGGTTTCAACGCAACACCTTAAATAAACCCGt                              >  1:787706/1‑65 (MQ=255)

                             gCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGt                              >  1:731217/1‑65 (MQ=255)

                             gCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGAAACACCTTAAATAAACCCgc                              >  1:3840513/1‑64 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGCTTa                         >  1:2525306/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:2916679/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:27963/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:2293474/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:2266417/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:3246889/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:3318777/1‑65 (MQ=255)

                                  gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTa                         >  1:3446144/1‑65 (MQ=255)

                                    cGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATg                       >  1:668930/1‑65 (MQ=255)

                                         aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGCTTATGGgccg                  >  1:1200332/1‑65 (MQ=255)

                                         aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                  >  1:3361319/1‑65 (MQ=255)

                                         aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                  >  1:1140528/1‑65 (MQ=255)

                                         aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                  >  1:541175/1‑65 (MQ=255)

                                         aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                  >  1:1503424/1‑65 (MQ=255)

                                         aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccg                  >  1:3748527/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:704519/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:1028176/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:3237255/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:1965010/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:1892339/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:1849115/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgc                 >  1:1404011/1‑65 (MQ=255)

                                          aCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCAGTGATTATGGgccgc                 >  1:1062427/1‑65 (MQ=255)

                                           cctgcctgCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                >  1:922049/1‑65 (MQ=255)

                                           cctgcctgCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                >  1:3212312/1‑65 (MQ=255)

                                           cctgcctgCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                >  1:1369923/1‑65 (MQ=255)

                                           cctgcctgCCGATCACGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgccgcc                >  1:2972376/1‑65 (MQ=255)

                                           ccggcctgccGCTCGCGCCAGTTTTAAACGCAACCCCTTAAATAACCCCGTAATTATGGgccgcc                <  1:1827478/62‑1 (MQ=255)

                                            ctgcctgcCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa               >  1:2335769/1‑65 (MQ=255)

                                            ctgcctgcCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCa               >  1:2115712/1‑65 (MQ=255)

                                            ctgcctgcCGATCGCGCCAGGTTTAAACGCAACAACTTAAATAAACCCGTGATTCTGGgccgccc               >  1:3770228/1‑64 (MQ=255)

                                                ctgccGATCGCGCCAGGTTAAAACGCACCCCCTTAAAAAACCCCGTGATTATGGGCCGCCATAcc           <  1:198242/65‑1 (MQ=255)

                                                cggccGATCGCGCAAGTTTTAAACGCAACACCTAAAATAACCCCGTGATTATGGGCCGCCATAcc           <  1:971124/63‑1 (MQ=255)

                                                  gccGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGACATACCTg         >  1:1020001/1‑65 (MQ=255)

                                                    ctaTCGCGCCAGGTTTAAACGCAACACCTAAAAAAAACCCGTGATTATGGGCCGCCATACCTggg       <  1:3682327/63‑1 (MQ=255)

                                                    cGATCGCGCCAGGTTTAAACGCAACACCTTAAATAACCCGGTGATTATGGGCCGCCATACCTggg       <  1:362854/65‑1 (MQ=255)

                                                      aTCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCGGGGATTATGGGCCGCAATACCTGGGaa     >  1:961259/1‑65 (MQ=255)

                                                        cgcGCCAGGTTTAAACGCACCACCTAAAAAAACCCCGTTATTATGGCCCGCCATACCTGGGaatc   <  1:1590160/63‑1 (MQ=255)

                                                         gcGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGaatca  >  1:3865795/2‑65 (MQ=255)

CGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAATCA  >  NC_000913/49848‑49969

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