BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,500,309	Δ1 bp	coding (333/363 nt)	b3372 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,500,308	0	C	.	99.2%	570.0 / ‑3.2	126	coding (332/363 nt)	b3372	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (1/0); new base . (63/62); total (64/62)
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.72e-01

Only 100 of 136 total aligned reads displayed.

GCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACT  >  NC_000913/3500233‑3500378

gcgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgc                                                                        >  1:3013288/1‑76 (MQ=255)

gcgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgc                                                                        >  2:787271/1‑76 (MQ=255)

gcgAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTAGACGCTGATTTccgc                                                                        >  2:952606/1‑76 (MQ=255)

  gAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcgg                                                                      >  2:4805395/1‑74 (MQ=255)

  gAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcgg                                                                      >  2:4689473/1‑74 (MQ=255)

  gAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcgg                                                                      >  2:1921459/1‑74 (MQ=255)

  gAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcgg                                                                      >  1:2615615/1‑74 (MQ=255)

   aaaGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTccgcggc                                                                     <  1:2189070/76‑4 (MQ=255)

     aGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGccc                                                                  >  1:3707714/1‑76 (MQ=255)

       cttCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCat                                                                <  1:3232605/74‑2 (MQ=255)

       ccTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAc                                                                <  1:4728920/76‑1 (MQ=255)

       ccTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAc                                                                >  1:1211836/1‑76 (MQ=255)

       ccTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCCGATTTCCG‑CGGCCCAc                                                                >  2:3762635/1‑76 (MQ=255)

        cTCGACATGATCAAGCTGGCGATGGCTATGGCCCCAGCGATGACCGCGGGTTATCCGCTGATTTCCG‑CGGACCAcg                                                               >  2:2656384/1‑76 (MQ=255)

        cTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcg                                                               >  2:2615291/1‑76 (MQ=255)

          cGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcgcg                                                             <  2:1049265/76‑1 (MQ=255)

          cGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCAcgcg                                                             <  2:4004049/76‑1 (MQ=255)

           gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                            <  1:2073024/76‑1 (MQ=255)

           gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                            <  1:2401126/76‑1 (MQ=255)

           gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                            <  1:1157363/76‑1 (MQ=255)

           gACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCgg                                                            <  2:4360152/76‑1 (MQ=255)

            aCATGATCAAGCTGGCGCTGGATATGGCAAAAGAGATGCACGCGGGTTATACGCTGATTTCCG‑CGTACCACGCggg                                                           >  1:1839861/1‑76 (MQ=255)

            aCATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCggg                                                           >  1:2279877/1‑76 (MQ=255)

            aCATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCggg                                                           <  2:3794124/76‑1 (MQ=255)

            aCATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCggg                                                           >  1:1802793/1‑76 (MQ=255)

              aTGATCAGGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCt                                                         <  1:4444929/76‑1 (MQ=255)

              aTGATCAAGCTTGCGATGGATATGGCAAAAGAGATGAACTCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCt                                                         <  1:2780490/76‑1 (MQ=255)

              aTGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCt                                                         <  1:1622309/76‑1 (MQ=255)

                 aTCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGACCACGCGGGCTAtc                                                      >  1:35376/1‑76 (MQ=255)

                  tCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtct                                                     <  2:1909808/76‑1 (MQ=255)

                   cAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtctc                                                    >  2:1597224/1‑76 (MQ=255)

                    aaGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTAtctcc                                                   >  1:61930/1‑75 (MQ=255)

                    aaGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCa                                                   >  1:1135688/1‑76 (MQ=255)

                      gCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCAcc                                                 >  1:1370146/1‑75 (MQ=255)

                      gCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACg                                                 <  2:3538219/76‑1 (MQ=255)

                       cTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGc                                                <  1:198063/76‑1 (MQ=255)

                            gATGGATATGGCAAAAGAGATGATCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCt                                           <  1:1543845/76‑1 (MQ=255)

                             aTGGATATGGCAAAAGAGATGACCTCTGGTTATACGCTGCTTTCCG‑CGGCCCACTCGGGCTATCTCACGCCACCTa                                          <  2:2970515/76‑1 (MQ=255)

                             aTGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTa                                          <  1:1046823/76‑1 (MQ=255)

                               ggATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAAt                                        >  1:3064245/1‑76 (MQ=255)

                                gATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATg                                       <  2:4396091/76‑1 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCCCCCAATGtt                                     >  2:4814276/1‑76 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                     >  1:1004921/1‑76 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                     <  2:2516701/76‑1 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                     <  2:3123580/76‑1 (MQ=255)

                                  tatGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                     <  1:1658236/76‑1 (MQ=255)

                                  tagggCAAAAGAGAGGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGtt                                     <  2:156684/73‑1 (MQ=255)

                                      gCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCCCGCGGGCTATCTCCCCCCCCCCACTGTTATcg                                 >  2:144955/1‑75 (MQ=255)

                                      gCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCCCCTAATGCTATCt                                 >  2:137169/1‑76 (MQ=255)

                                        aaaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTgg                               >  2:1067478/1‑76 (MQ=255)

                                         aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCAc                                             >  2:4035534/1‑61 (MQ=255)

                                         aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCAc                                             <  1:4035534/61‑1 (MQ=255)

                                         aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTggg                              <  2:3188972/76‑1 (MQ=255)

                                         aaaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGCTATCTggg                              >  1:1525763/1‑76 (MQ=255)

                                          aaGAGATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTgggg                             >  1:3120443/1‑76 (MQ=255)

                                           agagATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGa                            <  1:3274344/76‑1 (MQ=255)

                                           agagATGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGa                            >  1:1854010/1‑76 (MQ=255)

                                            gagaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                           >  1:2464654/1‑76 (MQ=255)

                                            gagaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                           <  2:935641/76‑1 (MQ=255)

                                            gagaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                           >  1:4676585/1‑76 (MQ=255)

                                            gagaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGAc                           >  1:1804756/1‑76 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTggggccg                          <  1:2123828/76‑4 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACg                          <  1:2896036/76‑1 (MQ=255)

                                             agaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACg                          <  1:3439554/76‑1 (MQ=255)

                                             agaTGAACGCGGGTGATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACg                          <  1:2904886/76‑1 (MQ=255)

                                              gaTGAACGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACgg                         <  2:2492986/76‑1 (MQ=255)

                                                  aaCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTgg                     <  2:2685558/76‑1 (MQ=255)

                                                  aaCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTgg                     <  2:4650281/76‑1 (MQ=255)

                                                  aaCGCGGGTTATACGCGGGTTCCCG‑CGGCCCTCGCGGGCTATCTCACGCCGACTAATGTTATCTGGGGACGGCTgg                     <  2:3099518/76‑1 (MQ=255)

                                                   aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGc                    <  2:2975374/76‑1 (MQ=255)

                                                   aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGc                    >  1:4599385/1‑76 (MQ=255)

                                                   aCGCGGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGc                    >  1:2015392/1‑76 (MQ=255)

                                                     gcgGGTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCag                  >  2:4077386/1‑76 (MQ=255)

                                                      cgGGGTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCGGGCaga                 >  2:951704/1‑76 (MQ=255)

                                                       gggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTCATGTTATCTGGGGACGGCTGGCagag                >  1:2251618/1‑76 (MQ=255)

                                                       gggTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGCCGGCTGGCagag                <  1:1067479/76‑1 (MQ=255)

                                                         gTTATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGaa              >  2:651689/1‑76 (MQ=255)

                                                          ttATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAc             >  2:1184026/1‑76 (MQ=255)

                                                          ttATACGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAc             >  1:2016208/1‑76 (MQ=255)

                                                           tataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAcc            <  2:3068334/76‑1 (MQ=255)

                                                           tataCGCTGATTTCC‑CCGGCCCTCGCGGGCTATCTCACGCCACCTAATGTTATCTGGGCACGGCCGGCAGCGAAcc            >  2:2133471/1‑76 (MQ=255)

                                                            ataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAAcca           >  2:1438962/1‑75 (MQ=255)

                                                            ataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCt           >  2:4465753/1‑76 (MQ=255)

                                                            ataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCt           >  2:1371460/1‑76 (MQ=255)

                                                            ataCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCt           >  1:946517/1‑76 (MQ=255)

                                                              aCGCTGCTTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa         <  2:4286019/76‑1 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa         >  1:3730246/1‑76 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa         >  2:4173363/1‑76 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa         <  2:1233540/76‑1 (MQ=255)

                                                              aCGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTaa         >  1:1662240/1‑76 (MQ=255)

                                                               cGCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAg        >  2:1382475/1‑76 (MQ=255)

                                                                gCTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGc       >  1:3404359/1‑76 (MQ=255)

                                                                 cTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTCAGCg      >  2:4751916/1‑76 (MQ=255)

                                                                 cTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCg      <  2:873802/76‑1 (MQ=255)

                                                                 cTGATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCg      <  1:431203/76‑1 (MQ=255)

                                                                  tGATTGCCG‑CGGCCCACGCGGGCTATCTCTCGCCACCTAATGTGATCTGGGCACGGCTGGCAGAGGACGTAAGCGa     >  1:2888200/1‑76 (MQ=255)

                                                                   gATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGaa    >  2:3044288/1‑76 (MQ=255)

                                                                   gATTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGaa    <  2:2439269/76‑1 (MQ=255)

                                                                    aTTTCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAAc   >  2:3733748/1‑76 (MQ=255)

                                                                     tttCCG‑CGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGCACCTAAGCGCACt  >  2:2363356/1‑76 (MQ=255)

GCGAAAGCCTCGACATGATCAAGCTGGCGATGGATATGGCAAAAGAGATGAACGCGGGTTATACGCTGATTTCCGCCGGCCCACGCGGGCTATCTCACGCCACCTAATGTTATCTGGGGACGGCTGGCAGAGAACCTAAGCGAACT  >  NC_000913/3500233‑3500378

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