BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,722,918	T→A	N247Y (AAC→TAC)	b1647 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,722,918	0	T	A	99.1%	309.3 / NA	107	N247Y (AAC→TAC)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (66/40); minor base C (1/0); total (67/40)
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.98e-01

Only 100 of 111 total aligned reads displayed.

GCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGG  >  NC_000913/1722845‑1722988

gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTaaa                                                                      >  1:3846323/1‑76 (MQ=255)

gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTaaa                                                                      >  2:2244787/1‑76 (MQ=255)

 cccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                     >  2:3461886/1‑76 (MQ=255)

 cccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                     <  1:4125968/76‑1 (MQ=255)

  ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTCAacc                                                                    >  2:5166705/1‑76 (MQ=255)

   cGTGCGCGCTCAATTTGACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAacca                                                                   <  2:1008458/76‑1 (MQ=255)

   cGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGAGGGTAAACGTAATACCCAGGCGCAAaccc                                                                   >  2:333273/1‑75 (MQ=255)

    gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                  >  1:4902072/1‑76 (MQ=255)

    gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccac                                                                  <  1:3206384/76‑1 (MQ=255)

     tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                 <  1:1190881/76‑1 (MQ=255)

     tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAAGCGCAAaccacc                                                                 >  2:1982625/1‑76 (MQ=255)

      aCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                <  1:605936/76‑1 (MQ=255)

       cgcgCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTg                                                               >  2:1050382/1‑76 (MQ=255)

        gcgcTCAATTTTACCTGAACAGATCATTGGCAGCGGCTGCGATGGTAAACGTAATAACCAGGCTTCAACCACCTGt                                                              >  2:3581755/1‑76 (MQ=255)

         cgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGtt                                                             >  1:346042/1‑76 (MQ=255)

          gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                            <  1:3776681/76‑1 (MQ=255)

          gcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTc                                                            <  2:2969653/76‑1 (MQ=255)

            tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                          >  1:4619288/1‑76 (MQ=255)

            tCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                          <  1:1838107/76‑1 (MQ=255)

            tCAATTTTACCTGAACAGATAATTGGCAGCGGCTGCGATGGTCAACGTAATACCCAGGCGTAAACCACCTGTTCaa                                                          >  1:4617210/1‑76 (MQ=255)

             cAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                         <  1:5525324/76‑1 (MQ=255)

              aaTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                        <  2:1928871/76‑1 (MQ=255)

               aTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCg                                                       >  2:1919543/1‑76 (MQ=255)

                ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                      >  2:1356195/1‑76 (MQ=255)

                ttttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                      >  1:5366165/1‑76 (MQ=255)

                 tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGCAATACCCAGGCGGAAACCACCCGTTCAATCGac                                                     >  2:2884373/1‑75 (MQ=255)

                  ttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCCCCTGTTCAATCGAGc                                                    >  2:5499253/1‑76 (MQ=255)

                   tACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                   >  2:3693174/1‑76 (MQ=255)

                   tACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGcc                                                   <  2:5035442/76‑1 (MQ=255)

                     ccTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcg                                                 <  1:233424/76‑1 (MQ=255)

                      cTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgc                                                >  1:3819380/1‑76 (MQ=255)

                       tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                      <  2:259847/69‑1 (MQ=255)

                       tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGa                                                      >  1:259847/1‑69 (MQ=255)

                        gAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                              <  1:5411994/76‑1 (MQ=255)

                        gAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGt                                              >  2:4018219/1‑76 (MQ=255)

                        caaCCGATAATTGGCAGCAGCTGCGTCGGTAAACGTAATCACCACGCGACAACAACCGAGTCAGTCGACCCggcgc                                              >  2:2018803/2‑75 (MQ=255)

                         aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCcgcgca                                             >  1:2703237/1‑74 (MQ=255)

                         aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                             >  2:1450257/1‑76 (MQ=255)

                           ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                           >  1:2774203/1‑76 (MQ=255)

                           ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                           >  1:5734966/1‑76 (MQ=255)

                           ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                           >  1:5782477/1‑76 (MQ=255)

                            cGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTaa                                          <  1:3461886/76‑1 (MQ=255)

                             gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAc                                         >  1:5395688/1‑76 (MQ=255)

                              aTAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAct                                        <  2:1539806/76‑1 (MQ=255)

                              aTAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAAct                                        >  1:2814476/1‑76 (MQ=255)

                               tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCTTTTAActc                                       <  1:4023300/76‑1 (MQ=255)

                               tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActc                                       <  1:2150513/76‑1 (MQ=255)

                               tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActc                                       <  1:2931206/76‑1 (MQ=255)

                                aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                      >  2:3054072/1‑76 (MQ=255)

                                aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActca                                      <  2:5395687/76‑2 (MQ=255)

                                 aTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                     >  1:5122931/1‑76 (MQ=255)

                                 aTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcct                                     <  2:3274913/76‑1 (MQ=255)

                                   tGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGTCCGCCTTTAATCCctct                                   >  1:5547013/1‑76 (MQ=255)

                                    ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                  >  2:2089359/1‑76 (MQ=255)

                                    ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                  <  2:2047548/76‑1 (MQ=255)

                                    ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                  >  2:4500371/1‑76 (MQ=255)

                                    ggCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTg                                  >  1:3376462/1‑76 (MQ=255)

                                     gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCCGc                                 >  2:5554085/1‑76 (MQ=255)

                                       aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccc                               >  1:5532697/1‑75 (MQ=255)

                                       aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccc                               >  2:5016728/1‑75 (MQ=255)

                                       aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCCCGCTTCACTCCTCTGccc                               >  2:5766849/1‑75 (MQ=255)

                                         cGGCTGCGATGGGAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacc                             >  1:5386699/1‑76 (MQ=255)

                                          ggCTGCGCTGGTAAGCGTAATTCCCAGGCGTAAGCCACCTGTTCACTCGAGCCCGCGTTTAACTCCTCTGccccca                            <  1:2998105/76‑1 (MQ=255)

                                            atGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCCCcaca                          >  2:2072741/2‑76 (MQ=255)

                                             tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTCCCACcccac                         >  2:61123/1‑75 (MQ=255)

                                              gCGATGGTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                        <  1:61123/76‑1 (MQ=255)

                                              gCGATGGTAAACGTAATACCCAGGCGTAACCCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                        <  1:1666172/76‑1 (MQ=255)

                                              gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                        <  1:4358399/76‑1 (MQ=255)

                                              gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGc                        >  2:4842890/1‑76 (MQ=255)

                                               cGAGGGAAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGcc                       <  2:4232392/76‑1 (MQ=255)

                                                gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCa                      >  1:4185275/1‑76 (MQ=255)

                                                gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCa                      >  1:2972894/1‑76 (MQ=255)

                                                gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCa                      >  1:1552112/1‑76 (MQ=255)

                                                 aTGGTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAg                     <  1:2312741/76‑1 (MQ=255)

                                                 aTGGTAAACGAAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAg                     <  1:5556462/76‑1 (MQ=255)

                                                  tGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGt                    >  2:3354166/1‑76 (MQ=255)

                                                   ggTAACCGTAAGACCCAGGCGTAACCCCCCTGTTCAACCGAGCCCGCGTTTAACTCCGCTGCCCCCACAGCCAGtt                   <  1:1750924/76‑1 (MQ=255)

                                                   ggTAAACGTAATCCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                   <  1:3844724/76‑1 (MQ=255)

                                                   ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                   >  1:2142099/1‑76 (MQ=255)

                                                   ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtt                   >  2:3533680/1‑76 (MQ=255)

                                                   ggTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGtc                   >  2:3290509/1‑75 (MQ=255)

                                                    gTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtc                  >  2:2263530/1‑76 (MQ=255)

                                                    gTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtc                  >  2:1810791/1‑76 (MQ=255)

                                                      aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcat                >  1:4249966/1‑76 (MQ=255)

                                                       aaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatc               >  1:791091/1‑76 (MQ=255)

                                                        aCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCTCGCGTTTAACTCCTCTGCCACCACAGCCCGCtcatca              >  1:2985774/1‑76 (MQ=255)

                                                        aCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatca              >  1:906423/1‑76 (MQ=255)

                                                         cgaaaTCCCCGGGCGTAACCCCCCTGTTCAATCAAGCCCGCGTTTAACCCCTCTGCCACCCCAGCCAGTTCATCAc             <  2:2203757/73‑1 (MQ=255)

                                                          gTAATACCCACGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGCTCATCACg            >  2:4519089/1‑76 (MQ=255)

                                                            aaTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCa          >  1:5702504/1‑76 (MQ=255)

                                                             aTACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAg         >  2:467011/1‑76 (MQ=255)

                                                              tccccGGGCGTAAACCCCCTGTTCAATCGAGCCCGCTTTTACCCCCTCTGCCCCCCCAGCCAGTTCATCACGCAGc        <  1:3385446/74‑1 (MQ=255)

                                                               aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg       >  1:3954397/1‑76 (MQ=255)

                                                               aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCg       >  2:2094665/1‑76 (MQ=255)

                                                               aCCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGATCATCACGCAGCg       >  1:5519877/1‑76 (MQ=255)

                                                                 ccAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGc     >  2:1908755/1‑76 (MQ=255)

                                                                  cAGGCGAAACCCCCCTTTTCAATCGACCCCGCTTTTACCTCCTCTGCCACCCCAGCCAGTTCATCACGCAGCGGCt    <  1:5078292/76‑1 (MQ=255)

                                                                   aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCCGCGGccg   >  2:1881262/1‑74 (MQ=255)

                                                                   aGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTg   >  2:4002022/1‑76 (MQ=255)

                                                                    ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTgg  >  2:2221625/1‑76 (MQ=255)

GCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGG  >  NC_000913/1722845‑1722988

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