BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,998,290	+G	coding (379/486 nt)	b4566 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,998,289	1	.	G	98.6%	380.6 / ‑1.3	141	coding (380/486 nt)	b4566	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/0); new base G (76/63); total (78/63)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

AATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCA‑GACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGT  >  NC_000913/3998214‑3998358

aaTCTCTTTACTGGTTTTGCATGCATCCAGACGAACGATGCCGACTTTATCTCTCGCGACATCAAACACGC‑CTCCa                                                                        >  1:707436/1‑76 (MQ=255)

aaTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCa                                                                        >  1:3515621/1‑76 (MQ=255)

 aTCTCTTTGCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCA‑g                                                                      <  2:689167/76‑1 (MQ=255)

 aTCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCA‑g                                                                      <  2:3667083/76‑1 (MQ=255)

  tctgTTTGCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CTCCA‑gg                                                                     <  2:5524798/76‑2 (MQ=255)

  tctcTTTACTGTTTTGGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCA‑gg                                                                     <  2:2819449/76‑2 (MQ=255)

        tACTGTTTTTGAAGGCATCCCGCCGAACGCTGCCGCGTTTACCTCCCGCTACATAAAGCACGC‑CTCCAGGACGGAt                                                                <  1:3281618/76‑1 (MQ=255)

        tACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGAt                                                                <  2:5249811/76‑1 (MQ=255)

         ncTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGAAATAACACACGA‑CGACAGGACGGATa                                                               >  1:3991294/2‑76 (MQ=255)

         aCTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATa                                                               >  2:579177/1‑76 (MQ=255)

          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATaa                                                              <  2:135063/76‑1 (MQ=255)

          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATaa                                                              <  1:4061967/76‑1 (MQ=255)

          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATaa                                                              <  1:1810781/76‑1 (MQ=255)

          cTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGA‑CGCCAGGACGGATaa                                                              <  1:983727/76‑1 (MQ=255)

            ggTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGAGAAAg                                                            >  2:1916986/1‑76 (MQ=255)

             gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGa                                                           <  1:4560040/76‑1 (MQ=255)

             gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGa                                                           >  1:2328683/1‑76 (MQ=255)

             gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGa                                                           >  1:6621/1‑76 (MQ=255)

             gTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGa                                                           <  2:5427035/76‑1 (MQ=255)

              ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGaa                                                          >  2:1850708/1‑76 (MQ=255)

              ttttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGaa                                                          >  2:3877756/1‑76 (MQ=255)

               tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGaaa                                                         <  1:5542701/76‑1 (MQ=255)

               tttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGaaa                                                         >  2:5453001/1‑76 (MQ=255)

               gttGAAGGCAGCCAGACGAACGATGCCGACGTTATCTACAACGACATAAAACAAGCACGCAG‑CACGGAgaaagaaa                                                         >  2:3871245/2‑76 (MQ=255)

                ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAAt                                                        <  2:4930871/76‑1 (MQ=255)

                ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAAt                                                        <  2:4927655/76‑1 (MQ=255)

                ttGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAAt                                                        >  1:1583257/1‑76 (MQ=255)

                 tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAAtt                                                       <  2:4424559/76‑1 (MQ=255)

                 tGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAAtt                                                       >  1:3112937/1‑76 (MQ=255)

                  gAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTc                                                      >  2:1206820/1‑76 (MQ=255)

                    aGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTccc                                                    <  1:3716785/76‑1 (MQ=255)

                     ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCt                                                   <  1:4546456/76‑1 (MQ=255)

                     ggCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCt                                                   <  1:2987412/76‑1 (MQ=255)

                      gCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCtt                                                  <  2:4501758/76‑1 (MQ=255)

                       cATCCAGACGAACGATGCCGACTTTAGCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTc                                                 <  2:2580436/76‑1 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGa                                                           >  2:5657339/1‑65 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGa                                                           <  1:5657341/65‑1 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGCAATTCCCTTCa                                                >  2:2084309/1‑76 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCa                                                >  2:3484781/1‑76 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCa                                                >  2:3143669/1‑76 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCGTCa                                                >  2:2630362/1‑76 (MQ=255)

                        aTCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAAAAGCC‑CCCCG‑GGCCGGTAAAGGAATTCCCCTCAc                                               >  2:4038424/1‑76 (MQ=255)

                         tCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCAc                                               <  1:416905/76‑1 (MQ=255)

                           cAGACGATCGATGCCGACTTTATCTCCCGCGCCATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCt                                             <  1:1026280/76‑1 (MQ=255)

                            aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCCCCTg                                            >  2:5743287/1‑76 (MQ=255)

                            aGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTg                                            <  1:3009512/76‑1 (MQ=255)

                              acgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGaa                                          <  1:5642546/76‑1 (MQ=255)

                               cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                         >  1:1533916/1‑76 (MQ=255)

                               cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAg                                         <  1:4810514/76‑1 (MQ=255)

                               cgaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGACAAAGAAATTCCCCTCACCTCaac                                         >  2:567558/1‑75 (MQ=255)

                                gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAg                                                            >  2:2424882/1‑56 (MQ=255)

                                gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGc                                        >  2:3761179/1‑76 (MQ=255)

                                gaacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGc                                        <  1:5234282/76‑1 (MQ=255)

                                 aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                       >  1:71247/1‑76 (MQ=255)

                                 aacgaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCg                                       >  2:4325944/1‑76 (MQ=255)

                                    gaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                    <  1:3179242/76‑1 (MQ=255)

                                    gaTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCgaga                                    >  2:5692816/1‑76 (MQ=255)

                                     aTGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAt                                   >  2:1152021/1‑76 (MQ=255)

                                      tGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAtt                                  >  1:2527479/1‑76 (MQ=255)

                                       gcagACGTTAGCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAATGAAATTCCCTTCACCTGAAGCGAGAatt                                 >  1:5412786/4‑76 (MQ=255)

                                       gCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAttt                                 >  2:1752626/1‑76 (MQ=255)

                                       gCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGAttt                                 <  2:4215105/76‑1 (MQ=255)

                                         cGACTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTaa                               >  2:4499544/1‑76 (MQ=255)

                                            cgttATCTCCCGCGACATAAAACACGG‑CGCCAGGGCGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGc                            <  2:5332721/74‑1 (MQ=255)

                                            cTTTATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGACGCGAGATTTAACGc                            >  2:1660937/1‑76 (MQ=255)

                                             tttATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGcc                           >  2:2706985/1‑76 (MQ=255)

                                              ttATCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCa                          >  1:4386645/1‑76 (MQ=255)

                                                aTCTCCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGATATTTAACGCCAcc                        >  2:2492055/1‑76 (MQ=255)

                                                aTCTCACGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCAcc                        <  1:1242054/76‑1 (MQ=255)

                                                 tctcCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCt                       >  2:2676322/1‑76 (MQ=255)

                                                  ctcCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                      >  1:1952242/1‑76 (MQ=255)

                                                  ctcCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                      >  2:5696014/1‑76 (MQ=255)

                                                  ctcCCGCGACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                      <  2:4105344/76‑1 (MQ=255)

                                                  ctcCCGCGACATAAAAAACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCtt                      <  2:3991293/76‑1 (MQ=255)

                                                       gcgACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                 >  2:3564435/1‑76 (MQ=255)

                                                       gcgACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgct                 <  1:1214994/76‑1 (MQ=255)

                                                        cgACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgctg                >  1:2957751/1‑76 (MQ=255)

                                                        cgACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgctg                <  1:579177/76‑1 (MQ=255)

                                                        cgACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgctg                >  1:932888/1‑76 (MQ=255)

                                                        cgACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTctgctg                >  2:4710959/1‑76 (MQ=255)

                                                         gACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGa               >  1:2574503/1‑76 (MQ=255)

                                                         gACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGa               >  2:318572/1‑76 (MQ=255)

                                                         gACATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGa               >  1:4298884/1‑76 (MQ=255)

                                                          aCATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAc              <  1:1547540/76‑1 (MQ=255)

                                                           cATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAct             <  2:151371/76‑2 (MQ=255)

                                                           cATAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGAcc             >  1:2368310/1‑76 (MQ=255)

                                                            aTAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGACGCGAGATGTAACGCCACCTTCTGCTGACCg            >  2:5177688/1‑76 (MQ=255)

                                                             tAAAACACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGt           >  1:1465971/1‑76 (MQ=255)

                                                               aaaCACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt         <  2:4393262/76‑1 (MQ=255)

                                                               aaaCACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGttt         >  1:1921563/1‑76 (MQ=255)

                                                                aaCACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt        >  2:5182482/1‑76 (MQ=255)

                                                                aaCACGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGtttt        <  1:3446180/76‑1 (MQ=255)

                                                                 acacGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc       >  1:3680975/1‑76 (MQ=255)

                                                                 acacGC‑CGCCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTc       >  2:1597189/1‑76 (MQ=255)

                                                                    cgc‑cgcCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt    <  1:554381/76‑1 (MQ=255)

                                                                    cgc‑cgcCAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGAt    <  1:1152021/76‑1 (MQ=255)

                                                                     gc‑cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATg   >  1:3834881/1‑76 (MQ=255)

                                                                      c‑cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGacgc  >  2:1541337/1‑73 (MQ=255)

                                                                      c‑cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGt  <  1:5696016/76‑1 (MQ=255)

                                                                      c‑cgccAGGACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCGGACCGGTTTCGATga  >  2:1129753/1‑75 (MQ=255)

AATCTCTTTACTGGTTTTGAAGGCATCCAGACGAACGATGCCGACTTTATCTCCCGCGACATAAAACACGC‑CGCCA‑GACGGATAAAGAAATTCCCTTCACCTGAAGCGAGATTTAACGCCACCTTCTGCTGACCGTTTTCGATGT  >  NC_000913/3998214‑3998358

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