BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,751,091	G→C	A426P (GCT→CCT)	recN →	RecN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,751,091	0	G	C	99.3%	373.3 / NA	143	A426P (GCT→CCT)	recN	RecN
Reads supporting (aligned to +/- strand): ref base G (0/1); major base C (73/68); minor base A (1/0); total (74/69)
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.99e-01

Only 100 of 143 total aligned reads displayed.

CATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGGCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGT  >  NC_000913/2751019‑2751162

gaTCGCCTGGGCTGTGGCGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCCGGTCAGCCCATGCAGCCTa                                                                      <  2:1125589/75‑1 (MQ=255)

cATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTa                                                                      <  1:3745468/76‑1 (MQ=255)

cATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTa                                                                      <  1:4506744/76‑1 (MQ=255)

cATCACCTGGGCGCTGACGGCGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCtc                                                                      >  2:5185292/1‑75 (MQ=255)

 aTCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTAt                                                                     >  1:5823872/1‑76 (MQ=255)

 aTCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTAt                                                                     >  1:16754/1‑76 (MQ=255)

  tCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTAtt                                                                    <  2:3098662/76‑1 (MQ=255)

  tCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTAtt                                                                    <  1:5760069/76‑1 (MQ=255)

   cACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCGCCAACCCAGGTCAGCCAATGCAGCCTATTg                                                                   <  1:2240782/76‑1 (MQ=255)

   cACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCCACCCAGGCCAGCCACTCCAGCCTCttc                                                                   >  2:620267/1‑75 (MQ=255)

   cACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTg                                                                   <  1:3015164/76‑1 (MQ=255)

   cACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTg                                                                   >  1:5562699/1‑76 (MQ=255)

    aCCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGc                                                                  >  1:3048391/1‑76 (MQ=255)

     cgtGGGCGCTGAGGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGcc                                                                 <  2:5103050/74‑1 (MQ=255)

      cTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAACGCAGCCTCTTGCCa                                                                >  2:4739079/1‑76 (MQ=255)

       tGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCaa                                                               <  1:1286713/76‑1 (MQ=255)

       tGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCaa                                                               >  2:2702907/1‑76 (MQ=255)

          gcgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGt                                                            <  1:1223827/76‑1 (MQ=255)

           cgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTc                                                           <  1:3663092/76‑1 (MQ=255)

           cgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTc                                                           <  1:2790398/76‑1 (MQ=255)

           cgcTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTc                                                           <  1:241722/76‑1 (MQ=255)

             cTGCCGTTGCGGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGc                                                         <  2:1527792/76‑1 (MQ=255)

             cTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGc                                                         >  1:4978520/1‑76 (MQ=255)

              tGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCa                                                        >  1:4555684/1‑76 (MQ=255)

                gcGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATc                                                      <  2:5375613/75‑1 (MQ=255)

                aCGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATc                                                      >  1:2562123/1‑76 (MQ=255)

                 cGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATcc                                                     <  1:187289/76‑1 (MQ=255)

                 cGGGGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATcc                                                     <  2:3683/76‑1 (MQ=255)

                  ggTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCg                                                    >  1:1130193/1‑76 (MQ=255)

                  ggTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCg                                                    <  1:353667/76‑1 (MQ=255)

                   gggCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCgg                                                   <  1:4170590/74‑1 (MQ=255)

                   gTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCgg                                                   <  2:5164226/76‑1 (MQ=255)

                    tGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggt                                                  >  2:3621552/1‑76 (MQ=255)

                    tGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggg                                                  >  1:102780/1‑75 (MQ=255)

                    tGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGCCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggt                                                  >  1:4851681/1‑76 (MQ=255)

                     gCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtg                                                 >  1:3083291/1‑76 (MQ=255)

                      ccGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtgg                                                >  1:404300/1‑76 (MQ=255)

                      ccGATCGGATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtgg                                                >  2:2912911/1‑76 (MQ=255)

                       cGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtggt                                               >  1:2013019/1‑76 (MQ=255)

                       cGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCggtggt                                               >  2:4561061/1‑76 (MQ=255)

                        gATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCCAAGTCGCATCCGgtggtg                                              >  2:4173219/1‑76 (MQ=255)

                         aTCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCCGCGGTGa                                             >  2:2201507/1‑76 (MQ=255)

                          tCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGaa                                            >  2:3522497/1‑76 (MQ=255)

                          tCGTATTGACTCTCGGGTAACCACCAACCCAGTCCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGaa                                            <  1:379356/76‑1 (MQ=255)

                           cGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGCGAAt                                           >  2:1508911/1‑76 (MQ=255)

                            gTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAAtt                                          >  2:5028623/1‑76 (MQ=255)

                            gTATTGAGTTTCGGGTAACCACCAACCCAGGGCAGCCAATGCAGCCTATTGCCAAAGTCGCATCAGGTGGTGAAtt                                          >  2:861342/1‑76 (MQ=255)

                             tATTGAGTGTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTg                                         <  2:1240693/76‑1 (MQ=255)

                              aTTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGt                                        >  1:2504308/1‑76 (MQ=255)

                              aTTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGt                                        >  1:2465052/1‑76 (MQ=255)

                                 gAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTccc                                     <  1:1322736/76‑1 (MQ=255)

                                   gTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGc                                   <  1:4655703/76‑1 (MQ=255)

                                   gTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGc                                   <  2:1540616/76‑1 (MQ=255)

                                   gTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGc                                   <  1:2792126/76‑1 (MQ=255)

                                      tCGGGTACCCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                <  1:303147/76‑1 (MQ=255)

                                      tCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                >  1:2235785/1‑76 (MQ=255)

                                      tCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                <  2:1914446/76‑1 (MQ=255)

                                      tCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                <  2:3644237/76‑1 (MQ=255)

                                      tCGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATc                                <  1:2912911/76‑1 (MQ=255)

                                       cGGGTACCCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCg                               <  2:2596807/76‑1 (MQ=255)

                                       cGGGTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCg                               <  1:5441883/76‑1 (MQ=255)

                                       cGGGGAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTccc                                     >  1:3998089/1‑70 (MQ=255)

                                        gggTAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGc                              <  1:2975491/76‑1 (MQ=255)

                                           tacccaccaGCCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACt                           <  1:5079707/73‑1 (MQ=255)

                                           tAACCACCAACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACt                           >  2:3402096/1‑76 (MQ=255)

                                              ccaccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGc                        >  1:1229327/1‑76 (MQ=255)

                                              ccaccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGc                        >  1:167797/1‑76 (MQ=255)

                                              ccaccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGc                        >  2:4761306/1‑76 (MQ=255)

                                               caccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCa                       >  1:4398295/1‑76 (MQ=255)

                                               caccaACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCa                       >  1:4407511/1‑76 (MQ=255)

                                                 ccaCCCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAAt                     <  2:2119027/76‑1 (MQ=255)

                                                  caACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATc                    <  1:3207908/76‑1 (MQ=255)

                                                  caACCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATc                    >  1:3469711/1‑76 (MQ=255)

                                                   aaCCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATcc                   <  1:1125575/76‑1 (MQ=255)

                                                   aaCCCAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATcc                   >  1:5516724/1‑76 (MQ=255)

                                                    ccccAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCa                  <  2:1992827/75‑1 (MQ=255)

                                                     cgcAGGTCAGCCAGTGCAGGCTATTGCCAGGGTGTGATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAg                 <  2:1032255/74‑1 (MQ=255)

                                                     cccAGGTCAGCGAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAg                 <  1:3307795/76‑1 (MQ=255)

                                                     cccAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAg                 >  2:830347/1‑76 (MQ=255)

                                                      ccAGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAgg                >  1:4573425/1‑76 (MQ=255)

                                                        aGGTGAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtc              <  2:5823871/76‑1 (MQ=255)

                                                        aGGTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtc              <  1:1631803/76‑1 (MQ=255)

                                                        aGGTCAGCCAATGCAGCCCATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtc              <  2:2465293/76‑1 (MQ=255)

                                                         ggTCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtca             >  1:4237684/1‑76 (MQ=255)

                                                           tCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatc           <  1:3321771/76‑1 (MQ=255)

                                                           tCAGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCAGCGCACTGGCAATCCAGGgcatc           >  1:3292350/1‑76 (MQ=255)

                                                           gcAGCGAAGGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGtcatc           <  2:16754/75‑1 (MQ=255)

                                                             aGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCCCTGGCAATCCAGGTCATCAc         >  2:1382166/1‑76 (MQ=255)

                                                             aGCCAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCAc         <  1:2618978/76‑1 (MQ=255)

                                                             aGCCAATGCAGCCTATTGCCAAAGGCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCATTCGCGGTCATCAc         >  1:5250208/1‑76 (MQ=255)

                                                               ccAATGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACgg       >  2:851419/1‑76 (MQ=255)

                                                                 aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcg     >  2:2756007/1‑76 (MQ=255)

                                                                 aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcg     >  2:4616124/1‑76 (MQ=255)

                                                                 aaTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcg     >  2:3851504/1‑76 (MQ=255)

                                                                  aTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACggagg    >  1:5821557/1‑73 (MQ=255)

                                                                  aTGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGgcgc    <  2:3240672/76‑1 (MQ=255)

                                                                   tGCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGGCACCCCGGcgcg   >  2:365885/1‑76 (MQ=255)

                                                                    gCAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGt  <  2:4207733/76‑1 (MQ=255)

                                                                    gCAGCCTATTGCCAAAATCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGt  >  2:4663803/1‑76 (MQ=255)

                                                                     cAGCCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAg                 <  2:872998/60‑1 (MQ=255)

CATCACCTGGGCGCTGACGGTGCCGATCGTATTGAGTTTCGGGTAACCACCAACCCAGGTCAGCCAATGCAGGCTATTGCCAAAGTCGCATCCGGTGGTGAATTGTCCCGCATCGCACTGGCAATCCAGGTCATCACGGCGCGT  >  NC_000913/2751019‑2751162

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