BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,905	T→G	L28R (CTC→CGC)	folA →	dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,905	0	T	G	99.5%	719.6 / ‑5.0	203	L28R (CTC→CGC)	folA	dihydrofolate reductase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (118/84); minor base A (1/0); total (119/84)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 255 total aligned reads displayed.

GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAA  >  NC_000913/49843‑49966

gTTGGCGGTAGATCGCGTTACCGGAATGAAAAACGCCATGCCGTGGAACCTGCCTGCCGATcgcg                                                             <  1:2026030/65‑1 (MQ=255)

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

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

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

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

       gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAgggt                                                      >  1:3566285/1‑63 (MQ=255)

       gTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAgggt                                                      >  1:832045/1‑63 (MQ=255)

         aGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTa                                                    >  1:3807140/1‑65 (MQ=255)

          taTCGCTTTATCGGCATGGAAACCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaa                                                   <  1:961317/64‑1 (MQ=255)

          gATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTaa                                                   >  1:2008216/1‑65 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

           aTCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGATTaaa                                                  >  1:3350876/1‑65 (MQ=255)

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

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

            tCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAAc                                                 >  1:2911475/1‑65 (MQ=255)

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

              gcgTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCAGCCGATCGCGCCAGGGTTAAACGc                                               >  1:1383892/1‑65 (MQ=255)

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

                gTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCaa                                             >  1:3736029/1‑65 (MQ=255)

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

                   aTCGGCATGGAAAACGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                          <  1:702623/65‑1 (MQ=255)

                   aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                          <  1:2827085/65‑1 (MQ=255)

                   aTCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAacac                                          <  1:232803/65‑1 (MQ=255)

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

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

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

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

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

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

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

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

                      ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCtg                                       >  1:977845/1‑64 (MQ=255)

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

                      ggCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGATTAAACGCAACACCtt                                       >  1:1348515/1‑65 (MQ=255)

                        cagggAAAACGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                     <  1:1682931/62‑1 (MQ=255)

                        cATGGAAACCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                     <  1:320618/65‑1 (MQ=255)

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

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

                        cATGGAAAACCCCTTGCGGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                     <  1:2928691/65‑1 (MQ=255)

                        cATGAAAACCGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaa                                     <  1:3136398/65‑1 (MQ=255)

                             aaacCGCCGTGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTtaaataaa                                <  1:2794944/65‑1 (MQ=255)

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

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

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

                              aaccGCCATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTCTAAACGCAACACCTTAAATAAAc                               <  1:691325/62‑1 (MQ=255)

                              aaccGCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAAc                               <  1:3877298/62‑1 (MQ=255)

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

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

                                accccATGCCGTGGACCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAAccc                             <  1:1419742/62‑1 (MQ=255)

                                 cGCCAGGCCGTGAAGCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCg                            <  1:3415754/65‑1 (MQ=255)

                                  gccgtgccgtGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGt                           <  1:2852364/65‑1 (MQ=255)

                                  gCCATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTCAATAAACCCGt                           >  1:1634641/1‑65 (MQ=255)

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

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

                                   ccATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTg                          >  1:3886957/1‑65 (MQ=255)

                                    cATGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGa                         <  1:1665733/65‑1 (MQ=255)

                                     aTGCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGAt                        <  1:881082/65‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                       gCCGTGGAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACACGTGATTa                      >  1:2774628/1‑65 (MQ=255)

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

                                            ggAACCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTCAAATAAACCCGTGATTATGGgc                 >  1:1808927/1‑65 (MQ=255)

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

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

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

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

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

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

                                              aaCCTGCCTGCCGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGgcag               >  1:438108/1‑63 (MQ=255)

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

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

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

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

                                                    cctgccGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAc         >  1:1955365/1‑65 (MQ=255)

                                                    cctgccGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATAc         >  1:3224482/1‑65 (MQ=255)

                                                    cctgccGATCGCGCCAGGCTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCCTaa         >  1:1971101/1‑64 (MQ=255)

                                                       gccGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTg      >  1:405255/1‑65 (MQ=255)

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

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

                                                         cGATCGCGCCAGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTggg    >  1:779477/1‑65 (MQ=255)

GTTAGCGGTAGATCGCGTTATCGGCATGGAAAACGCCATGCCGTGGAACCTGCCTGCCGATCTCGCCTGGTTTAAACGCAACACCTTAAATAAACCCGTGATTATGGGCCGCCATACCTGGGAA  >  NC_000913/49843‑49966

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