BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,334	A→G	V30V (GTT→GTC)	b2983 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,334	0	A	G	98.7%	407.0 / NA	156	V30V (GTT→GTC)	b2983	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (63/90); minor base C (2/0); total (65/91)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 157 total aligned reads displayed.

GTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACAACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTT  >  NC_000913/3130261‑3130405

  aCGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACAt                                                                     >  2:2589227/1‑76 (MQ=255)

  aCGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACCACAt                                                                     >  1:113750/1‑76 (MQ=255)

   cGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACAtt                                                                    <  1:3652197/76‑1 (MQ=255)

    gCCAAACATGGCGGTCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTa                                                                   <  2:3343315/76‑1 (MQ=255)

     ccAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTAc                                                                  <  2:4370322/76‑1 (MQ=255)

     ccAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTAc                                                                  <  1:3937720/76‑1 (MQ=255)

       aaaCATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTg                                                                >  1:4347541/1‑76 (MQ=255)

       aaaCATGGCGGGCGGCATCCCTTTACCGGCACACGAGAGCGCCCACAGACCCCGCAAGGCGCTCCCGACATTACTg                                                                >  2:1670365/1‑76 (MQ=255)

          cATGGCGTGCGTCATCCATTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAg                                                             <  1:3735585/76‑1 (MQ=255)

          cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAg                                                             <  2:3652922/76‑1 (MQ=255)

          cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAg                                                             <  1:3727796/76‑1 (MQ=255)

          cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAg                                                             <  1:726299/76‑1 (MQ=255)

          cATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACCACCTTACTTGAg                                                             >  1:1063997/1‑76 (MQ=255)

           aTGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGc                                                            >  1:2831697/1‑76 (MQ=255)

           aTGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGc                                                            <  2:3407253/76‑1 (MQ=255)

           aTGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGc                                                            <  1:1324255/76‑1 (MQ=255)

            tGGCGGGCGTCAGCCCTTTCCCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGcc                                                           <  2:168481/76‑1 (MQ=255)

            cggcgggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACCGGAGcc                                                           >  2:3623725/2‑76 (MQ=255)

             ggggggcgTCATCCCTTGACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCt                                                          <  2:1411602/73‑1 (MQ=255)

              gcgggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTa                                                         <  2:890809/76‑1 (MQ=255)

              gcgggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTa                                                         <  1:2501294/76‑1 (MQ=255)

               cgggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAACGCGCTCACGACATTACTGGAGCCTaa                                                        >  2:1337437/1‑76 (MQ=255)

                 ggcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAAcc                                                      <  2:2977714/76‑1 (MQ=255)

                  gcgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCa                                                     >  1:785596/1‑76 (MQ=255)

                   cgTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAg                                                    >  1:4766281/1‑76 (MQ=255)

                    gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAgg                                                   <  1:2652509/76‑1 (MQ=255)

                    gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAgg                                                   >  1:2033609/1‑76 (MQ=255)

                    gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAgg                                                   >  1:152617/1‑76 (MQ=255)

                    gTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAgc                                                   >  1:248796/1‑75 (MQ=255)

                     tCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgc                                                  <  2:4599420/76‑1 (MQ=255)

                     gcATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgc                                                  <  2:3211427/75‑1 (MQ=255)

                      cATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGgcg                                                 <  1:4225144/76‑1 (MQ=255)

                        tCCCTTTACCGGCACACGATGGCGCCAACAGACCCAGCAAGGCGCTCACGGCATTACTGGAGCCTAACCAGGCGCt                                               <  1:2924955/76‑1 (MQ=255)

                        tCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCt                                               <  1:4638708/76‑1 (MQ=255)

                        tCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCt                                               <  1:4547536/76‑1 (MQ=255)

                        gcccTTTGCCGGCACGCGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCt                                               <  2:845989/75‑1 (MQ=255)

                         cccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACAAGGCGCTa                                              <  1:3084466/76‑1 (MQ=255)

                          ccTTTGCCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                             <  1:4695205/76‑1 (MQ=255)

                          ccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                             >  2:2456563/1‑76 (MQ=255)

                          ccTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAt                                             <  1:502611/76‑1 (MQ=255)

                           cTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTAtt                                            >  2:3603563/1‑76 (MQ=255)

                             ttACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGc                                          >  2:1126504/1‑76 (MQ=255)

                             ttACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGc                                          <  2:3239622/76‑1 (MQ=255)

                              tACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCg                                         >  1:2012251/1‑76 (MQ=255)

                               aCCGTCACACGAGAGCGCGAACAGCCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGa                                        <  1:4015172/76‑1 (MQ=255)

                                ccGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaa                                       <  2:2822405/76‑1 (MQ=255)

                                 cGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaat                                      <  2:3335203/76‑1 (MQ=255)

                                  ggCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaata                                     >  2:1974373/1‑76 (MQ=255)

                                   gCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataa                                    >  1:1801711/1‑76 (MQ=255)

                                   gCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataa                                    >  2:1462315/1‑76 (MQ=255)

                                    cacgCGAGAGCACCACCAGCCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGaataat                                   <  1:2099246/76‑1 (MQ=255)

                                      cacGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTCTTGCGAATAATgg                                 >  2:829492/1‑76 (MQ=255)

                                      cacGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATgg                                 >  1:4574108/1‑76 (MQ=255)

                                       acGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGt                                >  1:2286646/1‑76 (MQ=255)

                                        cGAGAGCGCCAACAGCCCCAGCCAGGCGCTCCAGACCTTACCGGCGAATAACCCTGCGCTCTTGCGACTAATGGTg                               >  2:1505856/1‑76 (MQ=255)

                                        cGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTg                               <  1:4324101/76‑1 (MQ=255)

                                        cGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTg                               <  1:4175801/76‑1 (MQ=255)

                                        cGAGAGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTg                               >  2:1975639/1‑76 (MQ=255)

                                         gagaGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCCGGCGCTATTGCGAATAATGGTGc                              >  2:3144268/1‑76 (MQ=255)

                                         gagaGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                              <  1:1423170/76‑1 (MQ=255)

                                         gagaGCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGc                              <  2:3944449/76‑1 (MQ=255)

                                          ggagCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCg                             <  2:4460178/75‑1 (MQ=255)

                                          agagCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCg                             <  2:3098012/76‑1 (MQ=255)

                                           gagCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGa                            <  2:4178890/76‑1 (MQ=255)

                                           gagCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGa                            <  2:4182294/76‑1 (MQ=255)

                                           gagCGCCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAAACAGGCGCCATTGTGGATAATTGTGCGa                            >  1:1944388/1‑76 (MQ=255)

                                             gcgcCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGaaa                          <  2:2327805/76‑1 (MQ=255)

                                              cgcCACCAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTACCCAGGCGCTATTGCGAATAATGGTGCGaacc                         <  1:586164/76‑3 (MQ=255)

                                              cgcCAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAc                         <  1:1609725/76‑1 (MQ=255)

                                                 cAACAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAAcgcg                      >  1:869837/1‑76 (MQ=255)

                                                   aCAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGcc                    <  1:1535832/76‑1 (MQ=255)

                                                    cAGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCa                   <  1:3373604/76‑1 (MQ=255)

                                                     aGACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCAt                  >  1:637479/1‑76 (MQ=255)

                                                      gcaccAGAAAGCCGCTCACGACATGGCTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATc                 <  1:1085622/73‑1 (MQ=255)

                                                      gACCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATc                 >  1:3905975/1‑76 (MQ=255)

                                                       aCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                >  1:2778530/1‑76 (MQ=255)

                                                       aCCCAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCt                >  1:3498984/1‑76 (MQ=255)

                                                        cccAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTg               <  2:694740/76‑1 (MQ=255)

                                                        cccAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTg               <  1:1588316/76‑1 (MQ=255)

                                                         ccAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGc              >  2:659849/1‑76 (MQ=255)

                                                          cAGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCa             <  2:1055541/76‑1 (MQ=255)

                                                           aGCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaa            <  2:25530/76‑1 (MQ=255)

                                                            gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaaa           <  2:1339465/76‑1 (MQ=255)

                                                            gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaaa           <  1:559830/76‑1 (MQ=255)

                                                            gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaaa           <  1:3371595/76‑1 (MQ=255)

                                                            gCAAGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCaaa           >  2:1726464/1‑76 (MQ=255)

                                                             cAAGGCGCTCACGACATTACAGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAc          <  1:3892364/76‑1 (MQ=255)

                                                              aaGGCGCTCACTACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAcc         <  1:3217433/76‑1 (MQ=255)

                                                              aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAcc         >  1:836948/1‑76 (MQ=255)

                                                              aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAcc         >  2:2208561/1‑76 (MQ=255)

                                                              aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAcc         >  1:3903009/1‑76 (MQ=255)

                                                              aaGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAAcc         >  1:786476/1‑76 (MQ=255)

                                                               aGGCGCTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCa        >  2:1513894/1‑76 (MQ=255)

                                                                 gcgcTCGCGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAAt      <  2:3074289/76‑1 (MQ=255)

                                                                 gcgcCCACGACATTCCTGGAGCCTAACCAGGCGCTATTGCGAATAATGCTGCGAAACGCGCCATCTGCAAACCAAt      <  1:3605075/76‑1 (MQ=255)

                                                                  cgcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATg     >  2:3754498/1‑76 (MQ=255)

                                                                   gcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCGAACCAATGt    <  2:1232615/76‑1 (MQ=255)

                                                                   gcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCCATGt    >  1:254337/1‑76 (MQ=255)

                                                                   gcTCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGt    >  2:3438843/1‑76 (MQ=255)

                                                                     tCACGACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGttt  >  1:128042/1‑76 (MQ=255)

GTACGCCAAACATGGCGGGCGTCATCCCTTTACCGGCACACGAGAGCGCCAACAGACCCAGCAAGGCGCTCACAACATTACTGGAGCCTAACCAGGCGCTATTGCGAATAATGGTGCGAAACGCGCCATCTGCAAACCAATGTTT  >  NC_000913/3130261‑3130405

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