BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,759,336	C→T	92.0%	T45I (ACA→ATA)	yfjJ →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,759,336	0	C	T	92.0%	282.2 / 12.3	100	T45I (ACA→ATA)	yfjJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (5/3); new base T (50/42); total (55/45)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAC  >  W3110S.gb/2759270‑2759399

atgaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                             <  1:103305/71‑1 (MQ=255)

atgaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                             <  1:934775/71‑1 (MQ=255)

atgaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                             <  1:1203267/71‑1 (MQ=255)

atgaCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaa                                                             <  1:1733580/71‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:1772475/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:2865886/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:3073502/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:1357615/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:3659478/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:1999097/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:1990263/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:819478/70‑1 (MQ=255)

   aCTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:959738/70‑1 (MQ=255)

    cTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAata                                                         <  1:253349/71‑1 (MQ=255)

       ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGaaaa                                                           >  1:2263309/1‑66 (MQ=255)

       ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATg                                                       >  1:518566/1‑70 (MQ=255)

       ccTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATg                                                       >  1:2580776/1‑70 (MQ=255)

       ccTAAACCAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                       >  1:1479119/1‑70 (MQ=255)

           aaCGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGa                                                     <  1:1979846/68‑1 (MQ=255)

           aaCGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGa                                                     <  1:2112976/68‑1 (MQ=255)

           aaCGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGa                                                     <  1:1428278/68‑1 (MQ=255)

           aaCGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGa                                                     <  1:1356801/68‑1 (MQ=255)

                aCATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCAt                                             <  1:3692319/71‑1 (MQ=255)

                   tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                <  1:3585723/65‑1 (MQ=255)

                   tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                <  1:415828/65‑1 (MQ=255)

                   tCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAAcc                                                <  1:979485/65‑1 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCa                                              >  1:2885514/1‑65 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:289139/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:2550173/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:2314033/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:746424/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:1805192/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:1831859/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:3400136/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:821093/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:399528/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCttt                                        >  1:1511301/1‑71 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCtt                                         >  1:1262834/1‑70 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCtt                                         >  1:2883610/1‑70 (MQ=255)

                     ttAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCtt                                         >  1:1871708/1‑70 (MQ=255)

                       aaCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTc                                       >  1:2705621/1‑70 (MQ=255)

                       aaCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTc                                       >  1:1246400/1‑70 (MQ=255)

                        aCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                      <  1:1568853/70‑1 (MQ=255)

                        aCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                      <  1:3715850/70‑1 (MQ=255)

                        aCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                      <  1:308445/70‑1 (MQ=255)

                        aCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                      <  1:1672175/70‑1 (MQ=255)

                        aCTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                      <  1:759492/70‑1 (MQ=255)

                         cTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAata                                                         >  1:1895260/1‑50 (MQ=255)

                         cTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAata                                                         >  1:1277310/1‑50 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           >  1:3132311/1‑41 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:3684310/41‑1 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           >  1:3664420/1‑41 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:1022218/41‑1 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:393566/41‑1 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           <  1:524438/41‑1 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           >  1:913366/1‑41 (MQ=255)

                                aTTCTACGTGTTGATCTGCATGACCCAATTGATATAGaaaa                                                           >  1:282408/1‑41 (MQ=255)

                                 ttCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                            <  1:576560/71‑1 (MQ=255)

                                 ttCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                            <  1:3467420/71‑1 (MQ=255)

                                  tCTACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                       <  1:3211492/43‑1 (MQ=255)

                                  tCTAAGTGTTGATCTGCATGACCCAATTGATATAGAAAATATg                                                       <  1:3597025/43‑1 (MQ=255)

                                    tACGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGa                                                     >  1:3353360/1‑43 (MQ=255)

                                     aCGTGTTGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCa                                      <  1:3242454/57‑1 (MQ=255)

                                       gtgtTGATCTGCATGACCCAATTGATATAGAAAATATg                                                       <  1:2630276/38‑1 (MQ=255)

                                       gtgtTGATCTGCATGACCCAATTGATATAGAAAATATg                                                       <  1:592041/38‑1 (MQ=255)

                                         gtTGATCTGCATGACCCAATTGATATAGAAAATATg                                                       <  1:1824811/36‑1 (MQ=255)

                                          ttGATCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGt                   >  1:2016271/1‑71 (MQ=255)

                                            gATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:1811853/1‑39 (MQ=255)

                                            gATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:941276/1‑39 (MQ=255)

                                            gATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:3265599/1‑39 (MQ=255)

                                            gATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:3191136/1‑39 (MQ=255)

                                            gATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:1464071/1‑39 (MQ=255)

                                            gATCTGCATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:383744/1‑39 (MQ=255)

                                            gATCTGAATGACCCAATTGATATAGAAAATATGGATAAc                                                 >  1:3195442/1‑39 (MQ=38)

                                             aTCTGCATGACCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGc                 <  1:2590707/70‑1 (MQ=255)

                                                      aCCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAg                                >  1:1618786/1‑46 (MQ=255)

                                                      aCCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAg                                >  1:1834947/1‑46 (MQ=255)

                                                      aCCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAg                                >  1:958974/1‑46 (MQ=255)

                                                      aCCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAg                                >  1:253033/1‑46 (MQ=255)

                                                      aCCCAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAg                                >  1:3563769/1‑46 (MQ=255)

                                                         cAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCttt    <  1:1378535/71‑1 (MQ=255)

                                                         cAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCttt    <  1:1998824/71‑1 (MQ=255)

                                                         cAATTGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCttt    <  1:1233053/71‑1 (MQ=255)

                                                         cAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCttt    <  1:1069767/71‑1 (MQ=255)

                                                           aTTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTa   <  1:1171546/70‑1 (MQ=255)

                                                           aTTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:541081/1‑71 (MQ=255)

                                                           aTTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:1887283/1‑71 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGtt                            >  1:3579828/1‑43 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCt      >  1:2854807/1‑65 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCt      >  1:3164961/1‑65 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCt      >  1:282275/1‑65 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCt      >  1:790511/1‑65 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTa   >  1:1965909/1‑68 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:1627116/1‑69 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:1877809/1‑69 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:1245989/1‑69 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:458820/1‑69 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAc  >  1:271466/1‑69 (MQ=255)

                                                             tGATATAGAAAATATGGATAACCCATTATTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCt      >  1:2771093/1‑65 (MQ=255)

                                                              gATATAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCg        >  1:2542677/1‑62 (MQ=255)

ATGACTCCCTAAACGAACATCTTAACTCAATGATTCTACGTGTTGATCTGCATGACCCAATTGATACAGAAAATATGGATAACCCATTCTTTCAACCCAGGGTTGACTCTGGTGCTATATCTCGCTTTAC  >  W3110S.gb/2759270‑2759399

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