BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,717	+C	coding (86/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,717	1	.	C	99.1%	271.8 / NA	109	T29R (ACG→AGG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); major base C (44/64); minor base A (1/0); total (45/64)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 121 total aligned reads displayed.

ATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACG‑TGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCAT  >  NC_000913/4001642‑4001785

ataaataaAACCGAACAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                        <  2:2786815/76‑1 (MQ=255)

aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATTAAATTTAAGTACg                                                                        >  1:3655715/1‑76 (MQ=255)

aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATTAAATTTAAGTACg                                                                        >  1:3616523/1‑76 (MQ=255)

aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                        <  2:2777289/76‑1 (MQ=255)

aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                        >  2:2482096/1‑76 (MQ=255)

aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACg                                                                        <  1:5270483/76‑1 (MQ=255)

aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCGACTTGAAGTACg                                                                        >  2:1605812/1‑76 (MQ=255)

aTAAAGAAAAACGACCACATCACGCGATGCGCCAGGATTCCATCGGCGGGCACGTACAAAAACAACCTGAAGTACg                                                                        >  2:5548060/1‑76 (MQ=255)

 tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgc                                                                       <  1:1367667/76‑2 (MQ=255)

  aaaGAAAAACGACCAGATCACGCGATGCGTCATGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                      <  1:1333357/76‑3 (MQ=255)

  aaaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                      <  2:771493/76‑3 (MQ=255)

      aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                  <  2:2770242/76‑1 (MQ=255)

      aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                  <  2:2860037/76‑1 (MQ=255)

      aaaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                  <  2:2223006/76‑1 (MQ=255)

        aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                <  2:807828/76‑1 (MQ=255)

        aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                >  1:852018/1‑76 (MQ=255)

        aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                <  1:4044707/76‑1 (MQ=255)

         aaCGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt‑a                                                              <  2:2166455/76‑1 (MQ=255)

          aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt‑at                                                             <  2:4979855/76‑1 (MQ=255)

          aCGACCAGATCACGCGAGGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTAAGATGGCGCtgt                                                              >  1:5060836/1‑74 (MQ=255)

           cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt‑ata                                                            <  1:449935/76‑1 (MQ=255)

            gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATAc                                                           >  1:303964/1‑76 (MQ=255)

            gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATAc                                                           <  2:4389034/76‑1 (MQ=255)

            gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCCGGCGCT‑ATAc                                                           >  2:4223947/1‑76 (MQ=255)

            gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAAATTGAAGTACGCTGGCGCT‑ATAc                                                           >  1:4485637/1‑76 (MQ=255)

              ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATAccc                                                         <  2:4619488/76‑1 (MQ=255)

              ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATAccc                                                         <  2:4590633/76‑1 (MQ=255)

               cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt‑a                                                              <  1:2888150/70‑1 (MQ=255)

               cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt‑a                                                              >  2:2888150/1‑70 (MQ=255)

               cAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATAcccc                                                        >  1:4746885/1‑76 (MQ=255)

                aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCa                                                       <  2:3950324/76‑1 (MQ=255)

                aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCa                                                       <  2:5453848/76‑1 (MQ=255)

                 gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATAccccca                                                      >  2:2015122/1‑74 (MQ=255)

                  aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCaaa                                                     <  1:564732/76‑1 (MQ=255)

                  aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCaaa                                                     <  1:3861130/76‑1 (MQ=255)

                   tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCaaat                                                    <  2:1030591/76‑1 (MQ=255)

                    cACGCGATGCGTCAGGATTTCATCGGTGGGCACGTAGTAAATCTACTTGAAGTACGCTGGCGCT‑ATACCCCaaata                                                   <  1:5311238/76‑1 (MQ=255)

                    cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCaaata                                                   >  2:1371796/1‑76 (MQ=255)

                     aCGCGATGCGTCAGGATTTCATCGGCGGGAGCCGCGGAAATCAACTTTAAGTACGCTGGCGCTAATACCCaaaaaa                                                   >  1:5583338/1‑74 (MQ=255)

                      cgcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCaaataaa                                                 >  1:2831939/1‑76 (MQ=255)

                       gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATGCCCCaaataaaa                                                <  2:5327243/76‑1 (MQ=255)

                       gcgATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCTGGCGCT‑CTACCCCcaagaaaa                                                >  2:5792460/1‑76 (MQ=255)

                         gATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCaaataaaata                                              <  1:4948827/76‑1 (MQ=255)

                          aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAaataaaataa                                             <  2:5115440/76‑1 (MQ=255)

                           tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAg                                            <  1:1111424/76‑1 (MQ=255)

                            gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGc                                           <  1:4833960/76‑1 (MQ=255)

                              gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCgg                                         <  2:3105843/76‑1 (MQ=255)

                              gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCgg                                         <  1:3068723/76‑1 (MQ=255)

                               tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGc                                        <  1:4701714/76‑1 (MQ=255)

                               tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGc                                        >  2:1735101/1‑76 (MQ=255)

                                cAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGTCGCT‑ATACCCCAAATAAACTAAGCGGCa                                       <  1:4451420/76‑1 (MQ=255)

                                 aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCCAATAAAATAAACCGCaa                                      >  2:3718339/1‑76 (MQ=255)

                                 aGGATTTCATCGGCGGGCACGCACTAAATCAACTTGAAGGACGCTGGCGCT‑ATACCCCCAATAAAACAAGCGGCaa                                      >  2:1573648/1‑76 (MQ=255)

                                  ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAag                                     <  2:3839422/76‑1 (MQ=255)

                                  ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAag                                     <  1:2068230/76‑1 (MQ=255)

                                   gATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAaga                                    <  2:2772129/76‑1 (MQ=255)

                                    aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAagag                                   <  1:4109571/76‑1 (MQ=255)

                                    aTTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAagag                                   >  2:2540645/1‑76 (MQ=255)

                                     tttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGc                                  <  1:5528514/76‑1 (MQ=255)

                                     tttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGc                                  <  2:5805995/76‑1 (MQ=255)

                                     tttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGc                                  <  1:1394223/76‑1 (MQ=255)

                                       tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGa                                <  1:1284759/76‑1 (MQ=255)

                                       tCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGa                                <  1:1760358/76‑1 (MQ=255)

                                        cATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAg                               <  2:3655715/76‑1 (MQ=255)

                                        cATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAg                               <  2:3616523/76‑1 (MQ=255)

                                          tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTa                             <  2:5060835/76‑1 (MQ=255)

                                              ggggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                         <  2:4782634/75‑1 (MQ=255)

                                              cggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAtcc                         >  2:5582/1‑74 (MQ=255)

                                              cggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                         >  1:5330092/1‑76 (MQ=255)

                                              cggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                         <  1:4882668/76‑1 (MQ=255)

                                                 gCACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcc                      >  1:3742674/1‑76 (MQ=255)

                                                  cACGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                     >  2:4164210/1‑76 (MQ=255)

                                                   aCGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATACGCGGCAAGAGCGAGTAATACGCCCt                    >  2:4640377/1‑76 (MQ=255)

                                                   aCGTAGTAAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                    >  2:4583597/1‑76 (MQ=255)

                                                      tagtaAATCAGCTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcg                 <  2:4575757/76‑2 (MQ=255)

                                                      tagtaAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcc                 >  1:99048/1‑76 (MQ=255)

                                                       agtaAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                >  2:2955373/1‑76 (MQ=255)

                                                        gtaAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc               <  1:3210087/76‑1 (MQ=255)

                                                        gtaAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc               <  1:3389213/76‑1 (MQ=255)

                                                        gtaAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc               >  2:487713/1‑76 (MQ=255)

                                                         taAATCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAACGGCAAGAGCGAGTAATACGCCCTCCcggg              >  2:5432802/1‑74 (MQ=255)

                                                          aaaTCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt             >  2:815439/1‑76 (MQ=255)

                                                           aaTCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt            <  1:2210692/76‑1 (MQ=255)

                                                           aaTCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt            <  2:2990043/76‑1 (MQ=255)

                                                            aTCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt           >  1:2211432/1‑76 (MQ=255)

                                                            aTCAACTTAAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt           <  1:3282078/76‑1 (MQ=255)

                                                             tCAACTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg          <  2:5519336/76‑1 (MQ=255)

                                                                aCTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgttt       >  1:5188504/1‑76 (MQ=255)

                                                                 cTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                   <  2:4103460/63‑1 (MQ=255)

                                                                 cTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTg                   >  1:4103461/1‑63 (MQ=255)

                                                                 cTTGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt      <  1:1751686/76‑1 (MQ=255)

                                                                   tGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc    <  2:3647268/76‑1 (MQ=255)

                                                                   tGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc    <  1:2739142/76‑1 (MQ=255)

                                                                   tGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc    <  1:3718339/76‑1 (MQ=255)

                                                                   tGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc    <  2:1858337/76‑1 (MQ=255)

                                                                   tGAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGCTTGTTTTGc    >  2:3448368/1‑76 (MQ=255)

                                                                    gAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca   >  2:684295/1‑76 (MQ=255)

                                                                    gAAGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca   <  2:2723288/76‑1 (MQ=255)

                                                                     aaGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcat  >  1:721711/1‑76 (MQ=255)

                                                                      aGTACGCTGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTg     >  1:9115/1‑72 (MQ=255)

ATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACG‑TGGCGCT‑ATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCAT  >  NC_000913/4001642‑4001785

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