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.2%	374.6 / ‑5.7	127	N247Y (AAC→TAC)	b1647	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base A (81/45); minor base C (1/0); total (82/45)
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.25e-01

Only 100 of 129 total aligned reads displayed.

GCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGG  >  NC_000913/1722845‑1722988

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

gCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGGAAACGTAAAACCCAGGCGTaaa                                                                      >  2:480708/1‑76 (MQ=255)

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

 accGTAAGCGCTAAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAac                                                                     >  1:1213664/2‑76 (MQ=255)

  ccGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGGAAACGTAATACCCAGGCGTAAacc                                                                    >  1:534413/1‑76 (MQ=255)

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

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

    gTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGGAAACGTAATACCCAGGCCTAAaccac                                                                  >  2:4778955/1‑76 (MQ=255)

     tACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                 >  1:4017906/1‑76 (MQ=255)

     tACGCGCTCAATTTCACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAaccacc                                                                 >  1:3563333/1‑76 (MQ=255)

      aCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                >  1:2329354/1‑76 (MQ=255)

      aCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                >  1:549433/1‑76 (MQ=255)

      aCGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCt                                                                >  2:4246323/1‑76 (MQ=255)

       cgcgCTCCATGTCACCTGATCCGATAATTTGCAGCGGCTGCGATTGTAGACGTAATGCCCAGGCGTAAACCGCCTg                                                               <  2:2654049/76‑1 (MQ=255)

        gcgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGt                                                              >  1:2236791/1‑76 (MQ=255)

         cgcTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGtt                                                             <  2:3130869/76‑1 (MQ=255)

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

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

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

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

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

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

             cAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAAt                                                         >  1:3097656/1‑76 (MQ=255)

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

              aaTTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATc                                                        >  2:3649550/1‑76 (MQ=255)

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

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

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

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

                 tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                     >  1:270671/1‑76 (MQ=255)

                 tttACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAg                                                     <  1:2170459/76‑1 (MQ=255)

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

                    aCCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGccc                                                  <  1:1803819/76‑1 (MQ=255)

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

                       tGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAAGCGTAAACCACCTGTTCAAACGCGCCcgcg                                               >  2:3535392/1‑76 (MQ=255)

                        gAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGCGCCCGCGt                                              >  2:393828/1‑76 (MQ=255)

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

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

                        gAACCGATAATTGGCAGCGGCGGACATCGTAAACCTAATACCCAGGCGTAAACCCCCTGATCCACCGCGCCCGCGt                                              >  2:4710844/1‑76 (MQ=255)

                         aaCCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGtt                                             <  2:3302372/76‑1 (MQ=255)

                           ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                           <  2:1279068/76‑1 (MQ=255)

                           ccGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTa                                           <  2:1006705/76‑1 (MQ=255)

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

                            cGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTac                                          >  2:4446159/1‑75 (MQ=255)

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

                             gATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTCCAATCGAGCCCGCGTTTAAc                                         >  1:3647641/1‑76 (MQ=255)

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

                               tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActc                                       >  1:4218513/1‑76 (MQ=255)

                               tAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCCGTTCAATCCAGCCCGCGCTTAAatc                                       >  2:1365845/1‑76 (MQ=255)

                                aaTTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcc                                      <  2:4149241/76‑1 (MQ=255)

                                  ttGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAActcctc                                    >  2:3067807/1‑76 (MQ=255)

                                   tGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGCAAACCAACTGTTCAATCCAGCCCGCGTTTAACTCctct                                   >  2:3748845/1‑76 (MQ=255)

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

                                     gCAGCGGCTGCGATGGTAAACTTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGc                                 <  1:2669850/76‑1 (MQ=255)

                                     gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGCTTAACTCCTCTGc                                 >  1:404055/1‑76 (MQ=255)

                                     gCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGAGTTTAACTCCTCTGc                                 >  2:2630563/1‑76 (MQ=255)

                                       ggCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                               <  1:1777038/75‑1 (MQ=255)

                                       aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAACCCACCTGTTCCATCGAGCCCGCGTTTAACTCCTCTGcca                               <  2:2505435/76‑1 (MQ=255)

                                       aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                               >  1:2676136/1‑76 (MQ=255)

                                       aGCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGcca                               <  2:2479517/76‑1 (MQ=255)

                                        gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccac                              >  1:1189073/1‑76 (MQ=255)

                                        gCGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccac                              >  1:1255017/1‑76 (MQ=255)

                                         cGGCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacc                             <  2:1699568/76‑1 (MQ=255)

                                         cGGCTGCGATGGTAAACGTAACACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacc                             <  1:3405050/76‑1 (MQ=255)

                                          ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacca                            >  2:2858178/1‑76 (MQ=255)

                                          ggCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGccacca                            >  1:688493/1‑76 (MQ=255)

                                           gCTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCcaccac                           >  2:2373606/1‑76 (MQ=255)

                                            cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                          <  1:4209094/76‑1 (MQ=255)

                                            cTGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACcaca                          >  1:4174792/1‑76 (MQ=255)

                                             tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCCCcacac                         >  2:767648/1‑75 (MQ=255)

                                             tGCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAg                         >  2:3205829/1‑76 (MQ=255)

                                              tcGATGGTAAACGTAATCCCCAGGCGTCAACCACCGTTTCAATCGAGCTCGCGTTTACCTCATCTGCCACcacacc                        >  1:4576003/2‑74 (MQ=255)

                                              gCGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCCCCACAGc                        >  2:553542/1‑76 (MQ=255)

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

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

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

                                               cGATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGcc                       >  2:2600355/1‑76 (MQ=255)

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

                                                gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCa                      <  2:2492325/76‑1 (MQ=255)

                                                gATGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCa                      <  2:4060994/76‑1 (MQ=255)

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

                                                gATGGCAAACGTAATACCCACGCGGAAACCACCTGCACAATCGCGCCCTGCTTTCACACCTCTGCCACCACATCCa                      >  2:4611595/1‑76 (MQ=255)

                                                 aTGGTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAg                     <  2:3508878/76‑1 (MQ=255)

                                                   tgTAAACGACACACCCGGGCGAAAACCGCCTGTTCACTCGCGCCCGCGTTTAACCCCTCTGCCACCACAGCCAGtt                   <  1:1230455/75‑1 (MQ=255)

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

                                                    gTAAACGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCCCAGCCAGTtc                  >  2:1956353/1‑76 (MQ=255)

                                                    atAAACGAAATCCCCGGGCGAAACCCACCTGTTCTGTCGAGCCCGCGTTTAACCCCTCTCCCACCCCAGCCAGTtc                  <  1:2087509/75‑1 (MQ=255)

                                                     gaaaCGTAATCCCGAGGCGAAAACCCCCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCCCCCCAGCCAGTtca                 <  2:501778/75‑1 (MQ=255)

                                                      aaaCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcac                >  2:1031795/1‑75 (MQ=255)

                                                        aCGTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTtcatcc              >  2:3748003/1‑75 (MQ=255)

                                                         cGTAATACCCAGGCGTAACCCCCCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCAc             <  1:2137256/76‑1 (MQ=255)

                                                          gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACg            >  1:3701803/1‑76 (MQ=255)

                                                          gTAATACCCAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACg            >  2:2701671/1‑76 (MQ=255)

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

                                                              cccccAGGCGAAACCCCCCTGTTAAACCGACCCCGCTTTTAACCCCTCTCCCCCCCCAGCCAGTTCATCACGCAGc        <  2:3814781/74‑1 (MQ=255)

                                                                cccGGGCGTAACCCCCCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCCCAGCCAGTTCATCACGCAGCgg      <  1:3141129/76‑1 (MQ=255)

                                                                cccAGGCGTAAACAACCTGTTCAATCGAGCCCGCCTTTAACTCATCTGCCACCACAGTCATGTCAACGCGCAGcag      >  1:3779678/1‑74 (MQ=255)

                                                                  cAGGCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCt    >  2:1045739/1‑76 (MQ=255)

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

                                                                    ggCGTAAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCAAAGCCAGTTCATCACGCAGCGGCTgg  >  2:3886753/1‑76 (MQ=255)

GCCCGTACGCGCTCAATTTTACCTGAACCGATAATTGGCAGCGGCTGCGATGGTAAACGTAATACCCAGGCGTTAACCACCTGTTCAATCGAGCCCGCGTTTAACTCCTCTGCCACCACAGCCAGTTCATCACGCAGCGGCTGG  >  NC_000913/1722845‑1722988

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