BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,614,238	G→T	L133I (CTA→ATA)	yjjI ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,238	0	G	T	99.3%	472.4 / ‑6.0	149	L133I (CTA→ATA)	yjjI	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base T (79/69); total (80/69)
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.93e-01

Only 100 of 149 total aligned reads displayed.

AGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTA  >  NC_000913/4614166‑4614309

aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAt                                                                      <  1:3464393/76‑1 (MQ=255)

aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAt                                                                      >  1:3262477/1‑76 (MQ=255)

aGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCCACCATCGATCTCGCCCTGTGTTATAAt                                                                      >  2:804947/1‑76 (MQ=255)

  gCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAAttt                                                                    >  2:3937415/1‑75 (MQ=255)

   cggcTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                   <  2:1506552/73‑1 (MQ=255)

   cGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                   >  2:732798/1‑76 (MQ=255)

   cGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCt                                                                   <  2:1481687/76‑1 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCTAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                  <  2:2562861/76‑1 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                  <  1:702395/76‑1 (MQ=255)

    gTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTa                                                                  >  2:3118185/1‑76 (MQ=255)

       tGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACa                                                               >  2:2082073/1‑76 (MQ=255)

       tGGCAGGGTTCTGTCGAGGGAACGCCAGAAACGTTTTATTTGAACATCGATCTAGTCTTTTGTTATAATTCTAACa                                                               >  1:1404504/1‑76 (MQ=255)

        ggCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACAt                                                              >  2:2934701/1‑76 (MQ=255)

         gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                             >  1:4052235/1‑76 (MQ=255)

         gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                             >  2:4305469/1‑76 (MQ=255)

         gCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATa                                                             <  2:647080/76‑1 (MQ=255)

          cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                            <  2:2003276/76‑1 (MQ=255)

          cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                            >  2:1124275/1‑76 (MQ=255)

          cAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATAc                                                            <  1:1651754/76‑1 (MQ=255)

           aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                           >  1:645847/1‑76 (MQ=255)

           aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                           >  1:4448629/1‑76 (MQ=255)

           aGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACg                                                           <  1:1216671/76‑1 (MQ=255)

              gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGtt                                                        <  1:1819057/76‑1 (MQ=255)

              gTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATCCTAACATACGGtt                                                        >  2:1987570/1‑76 (MQ=255)

                tCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGc                                                      >  1:264076/1‑76 (MQ=255)

                    tCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                  <  1:1052747/76‑1 (MQ=255)

                    gcGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaaca                                                  <  2:4162797/75‑1 (MQ=255)

                     cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                 <  1:4166137/76‑1 (MQ=255)

                     cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                 <  2:3029021/76‑1 (MQ=255)

                     cGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGcaacaa                                                 >  1:192123/1‑76 (MQ=255)

                      gAGGTAACGCCAGAGACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                >  1:2957360/1‑76 (MQ=255)

                      gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                <  2:4111816/76‑1 (MQ=255)

                      gAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCaacaac                                                >  1:3200177/1‑76 (MQ=255)

                       aGGTAACGCCAGAAACGTTTTATTCGAACATTGATCCCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                               >  1:2975181/1‑76 (MQ=255)

                       aGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                               >  2:4119796/1‑76 (MQ=255)

                       aGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACg                                               >  1:2436108/1‑76 (MQ=255)

                          tAACGCCAGCAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                            <  1:2532731/76‑1 (MQ=255)

                           aaCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                            >  2:4406696/1‑75 (MQ=255)

                           aaCGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCAt                                            <  1:4406696/75‑1 (MQ=255)

                             cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                         <  2:1460918/76‑1 (MQ=255)

                             cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                         >  1:1495250/1‑76 (MQ=255)

                             cGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCa                                         >  2:3247213/1‑76 (MQ=255)

                               cgaGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGt                                       <  2:1726350/74‑1 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                      >  2:4778007/1‑76 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                      <  1:1179748/76‑1 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                      <  2:2664597/76‑1 (MQ=255)

                                cAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGtt                                      >  1:190299/1‑76 (MQ=255)

                                  gAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGc                                    >  1:2368718/1‑76 (MQ=255)

                                   aaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcc                                   >  2:50722/1‑76 (MQ=255)

                                   aaaCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcc                                   >  2:3921487/1‑76 (MQ=255)

                                    aaCGTTTTATTCGAACATCGATCTCGTCTTGTGTGATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGccc                                  <  1:1580848/76‑1 (MQ=255)

                                     aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAGCGCATCCAGTTGcccc                                 <  2:1708993/76‑1 (MQ=255)

                                     aCGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGcccc                                 >  2:332326/1‑76 (MQ=255)

                                      cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa                                >  2:829382/1‑76 (MQ=255)

                                      cGTTTTATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCa                                <  1:4408758/76‑1 (MQ=255)

                                        ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                              >  2:1231074/1‑76 (MQ=255)

                                        ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                              >  1:2267772/1‑76 (MQ=255)

                                        ttttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgg                              >  1:505594/1‑76 (MQ=255)

                                         tttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAgga                             >  1:475967/1‑75 (MQ=255)

                                         tttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGt                             >  1:2983325/1‑76 (MQ=255)

                                          ttATTCTAACCTCGTTCTCGTCTTGTGTTATAATTCTAACATACTGTTTCAACACCGCATCGAGTTGCCCCAGGTa                            <  1:2878317/76‑1 (MQ=255)

                                          ttATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTa                            >  1:909600/1‑76 (MQ=255)

                                           tATTCTAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAg                           <  1:1833197/76‑1 (MQ=255)

                                           tATTCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAg                           <  2:1931068/76‑1 (MQ=255)

                                              tCGAACATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGAcc                        <  1:1987054/76‑1 (MQ=255)

                                                 aaCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGc                     <  2:3204530/76‑1 (MQ=255)

                                                  aCATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCa                    >  2:1392756/1‑76 (MQ=255)

                                                   cATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                   >  2:153349/1‑76 (MQ=255)

                                                   cATCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCAt                   >  1:2508455/1‑76 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                  >  2:1423940/1‑76 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                  >  1:559596/1‑76 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                  <  2:1229163/76‑1 (MQ=255)

                                                    aTCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATc                  <  1:2571432/76‑1 (MQ=255)

                                                     tCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                 >  2:3406211/1‑76 (MQ=255)

                                                     tCGATCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCg                 >  2:1757500/1‑76 (MQ=255)

                                                     tCGATCACGTATTGTGGTAGAATTCTAACATACGGCCCCTACAACGCAACCAGTGGCCCCAGGCAGACCCGCATCg                 >  2:4490730/1‑76 (MQ=255)

                                                        aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg              >  2:227011/1‑76 (MQ=255)

                                                        aTCTCGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtg              <  1:3728609/76‑1 (MQ=255)

                                                         tctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgt             >  1:2689910/1‑76 (MQ=255)

                                                          ctcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGAtgtg            <  2:3262477/76‑1 (MQ=255)

                                                           tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa           <  2:3761784/76‑1 (MQ=255)

                                                           tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa           <  1:73369/76‑1 (MQ=255)

                                                           tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa           >  1:1097883/1‑76 (MQ=255)

                                                           tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa           >  1:1869533/1‑76 (MQ=255)

                                                           tcGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGa           <  1:3717239/76‑1 (MQ=255)

                                                            cGTTTTGTGTTATAATTTTAACATGCGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGCCCGGCATCGATGTGAc          <  2:1034688/76‑1 (MQ=255)

                                                            cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc          <  1:4233327/76‑1 (MQ=255)

                                                            cGTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAc          >  1:4123630/1‑76 (MQ=255)

                                                             gTCTTGTGTTATAATTCTAACATACGGTTGCAACAAGGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc         <  1:1074050/76‑1 (MQ=255)

                                                             gTCTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGAcc         >  2:2814047/1‑76 (MQ=255)

                                                               cTTGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGa       <  1:1636810/76‑1 (MQ=255)

                                                                ttGTGTTATAATTCTAACATACGTTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATGGATTTGACCgcc      <  1:1564033/76‑3 (MQ=255)

                                                                ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc      <  1:4087699/76‑1 (MQ=255)

                                                                ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc      <  1:4530502/76‑1 (MQ=255)

                                                                ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc      >  2:733191/1‑76 (MQ=255)

                                                                ttGTGTTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGAc      <  2:3605248/76‑1 (MQ=255)

                                                                  gtgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACgg    >  2:4134116/1‑76 (MQ=255)

                                                                   tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGt   >  1:3031111/1‑76 (MQ=255)

                                                                   tgtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGt   <  2:1609715/76‑1 (MQ=255)

                                                                    gtTATAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGGCCGACGGTa  <  1:3822199/76‑1 (MQ=255)

AGGCGTCTGGCAGGGTTCTGTCGAGGTAACGCCAGAAACGTTTTATTCGAACATCGATCTCGTCTTGTGTTAGAATTCTAACATACGGTTGCAACAACGCATCCAGTTGCCCCAGGTAGACCGGCATCGATGTGACCGACGGTA  >  NC_000913/4614166‑4614309

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