BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,596,118	A→T	6.4%	intergenic (‑15/+14)	leuP ← / ← leuQ	tRNA‑Leu/tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,596,118	0	A	T	6.4%	1124.0 / 8.3	390	intergenic (‑15/+14)	leuP/leuQ	tRNA‑Leu/tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base A (277/81); new base T (27/2); total (304/86)
Fisher's exact test for biased strand distribution p-value = 5.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 406 total aligned reads displayed.

AAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGA  >  CP009273/4596053‑4596185

aaGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGACTGAGAGAAatata                                                                    >  1:13531681/1‑67 (MQ=14)

aaGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                                                    <  1:26440609/68‑1 (MQ=255)

  gAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATTTCa                                                                   >  1:38986762/1‑67 (MQ=255)

  gAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                                                   >  1:5329961/1‑67 (MQ=255)

  gAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCa                                                                   >  1:21873365/1‑67 (MQ=255)

    acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGAGTGCGAGCAATTTCAACg                                                                >  1:44854089/1‑68 (MQ=255)

    acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCTGTGCGAGCAATATCAAc                                                                 <  1:35297043/67‑1 (MQ=255)

    acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACg                                                                >  1:1450624/1‑68 (MQ=255)

    acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGCCAGGAATATCaaag                                                                >  1:32715565/1‑66 (MQ=14)

      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGGCAGAGAAATATCAAAGTg                                                              >  1:12003510/1‑68 (MQ=14)

      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCGAGGACTATCTACGTg                                                              >  1:16843910/1‑68 (MQ=17)

      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGCGAGGAATTTCAAAGTg                                                              >  1:29764965/1‑68 (MQ=17)

       cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                                             <  1:11052252/68‑1 (MQ=255)

       cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTgg                                                             <  1:29047070/68‑1 (MQ=255)

       cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGAAATATCAACGTgg                                                             <  1:27170445/68‑1 (MQ=255)

         tacacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                                           <  1:34244885/67‑1 (MQ=255)

         aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                                           <  1:33580053/68‑1 (MQ=255)

         aaCACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtt                                                           <  1:23832093/68‑1 (MQ=255)

           cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGAAATATaaaa                                                                 >  1:49032535/1‑59 (MQ=17)

            acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCGACGGCGAGCAATATCAACGTGgttgtt                                                        >  1:38725017/1‑68 (MQ=18)

            acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGttttt                                                        >  1:15154466/1‑68 (MQ=255)

             ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGCGAGCAATTTCAACGTGGTTTTTg                                                       >  1:35906928/1‑68 (MQ=255)

               tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                                     <  1:27877280/68‑1 (MQ=255)

               tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                                     <  1:3229214/68‑1 (MQ=255)

               tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                                     <  1:23260555/68‑1 (MQ=255)

               tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                                                     <  1:28672449/68‑1 (MQ=255)

                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTTAGAGCTATTTCAACGTGTTTGTTGGTg                                                    >  1:13958632/1‑68 (MQ=17)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCTATTTCAAAGTGGTTTTTGGTGCg                                                  >  1:16667237/1‑68 (MQ=21)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                                  >  1:37390794/1‑68 (MQ=255)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                                  <  1:20324924/68‑1 (MQ=255)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                                                  <  1:4546782/68‑1 (MQ=255)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGAAATATATAGGTGGTTTTTgg                                                      >  1:50522604/1‑64 (MQ=18)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGGGGGAGATATTTAAAAGTGTTTTTTGGTGCg                                                  >  1:16081883/1‑68 (MQ=17)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGGGATCTATTTCAACGTGTTTTTTGGTGCg                                                  >  1:799854/1‑68 (MQ=18)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGCGAGCAATATCAACGTGGTTGTTGGTGCg                                                  >  1:31580729/1‑68 (MQ=255)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAAGTGCGAGCAATATAAACGTGGTTTTTGGTggg                                                  >  1:14881095/1‑66 (MQ=25)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAACTGTGAGGAATTTCAACGTGGTTGTTGGTGCg                                                  >  1:42156233/1‑68 (MQ=18)

                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCAACTGGCAGCAATATAAAAGAGGTTTTTGGTGCg                                                  >  1:43709450/1‑68 (MQ=17)

                     tAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑g                                               <  1:4694325/67‑1 (MQ=255)

                      aGCGCTTCTACCAATTCCTCCACCTTAGCCCATTGCTAGCAATATCAACGTGTTTTTTGGTGCGA‑ggg                                             <  1:16389461/68‑1 (MQ=37)

                        cgcgTCTACCAATTCCGCCCCCTTCGCCCAGTGGGAGCAATTTCAAAGTGGTGGTTGGTGCGA‑ggggg                                           >  1:22135096/1‑68 (MQ=21)

                        cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGGAATATAAACGTGGTTGTTGGTGCGA‑ggggg                                           >  1:39829719/1‑68 (MQ=37)

                        cgcgTCTACCAATTCCGCCACCTTCGCCAAGTGGGAGCAATATCAACGAGTTTTTTGCTGCGA‑ggggg                                           >  1:6470056/1‑68 (MQ=25)

                                 cAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGCCTTGa                                  >  1:18051863/1‑68 (MQ=255)

                                   aTTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGaa                                 >  1:5539350/1‑67 (MQ=255)

                                      ccGCCACCTTCGCCCCGTGCGATCAATACCAACGTGGTGTTTGGTGCGAGGGGGGGGGACTTGAAccc                              <  1:31094816/68‑1 (MQ=25)

                                          ccccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                         <  1:30086005/66‑1 (MQ=255)

                                           aCCTTCGCCCAGTGCTAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGt                        <  1:20903839/68‑1 (MQ=255)

                                           aCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                         >  1:12128503/1‑67 (MQ=255)

                                           aCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                         >  1:11515365/1‑67 (MQ=255)

                                           aCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                         >  1:6149929/1‑67 (MQ=255)

                                            ccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTc                       <  1:2488466/68‑1 (MQ=255)

                                             cTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACt                                     <  1:10481810/53‑1 (MQ=255)

                                              ttCGCCCAGTGCGAGCAATATAAACGGGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCg                     <  1:33947122/68‑1 (MQ=255)

                                                  cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                 >  1:13180075/1‑68 (MQ=255)

                                                  cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                 >  1:15558047/1‑68 (MQ=255)

                                                  cccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                 >  1:18947031/1‑68 (MQ=255)

                                                  ccaaGTGCGAGCAATATCAAAGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                 >  1:2192844/4‑68 (MQ=255)

                                                   ccAGTTAGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGTGGGGGAATTGAACCCCCACGTCAGTAAgt                >  1:21624270/1‑67 (MQ=21)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:42816972/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:26913957/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:12335599/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:37737812/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:45882408/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:17322853/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:46199999/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:8725855/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:7896130/1‑68 (MQ=255)

                                                   ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:699116/1‑68 (MQ=255)

                                                    cAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                >  1:9663917/1‑67 (MQ=255)

                                                       tGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAggcc              >  1:18426085/1‑64 (MQ=255)

                                                       tGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGac              >  1:19833383/1‑66 (MQ=255)

                                                        gCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGaca             <  1:25735100/66‑1 (MQ=255)

                                                        gCGAGCAATATAAACGTGGTTTTTGGGGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt           >  1:8381813/1‑68 (MQ=255)

                                                            gCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAac        >  1:7594866/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGACCCCCCACGTCCGTAATGACACTAaca       >  1:31413449/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACCCTAaca       >  1:50965845/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:50078367/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:17550558/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:18649847/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:22338621/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:25381155/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:28277298/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:36788785/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:41638174/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:40710559/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:42525121/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:40018012/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:39235173/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:38049258/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:37199113/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:48138744/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:34899671/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:31777903/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:47434831/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:47695573/1‑67 (MQ=255)

                                                             cAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACAAaca       >  1:41241530/1‑67 (MQ=255)

                                                               atatCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAaca       >  1:5817449/1‑65 (MQ=255)

                                                                 atCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGa  >  1:24449857/1‑68 (MQ=255)

                                                                 atCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACAATAACACCTGa  >  1:30494391/1‑68 (MQ=255)

AAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGA  >  CP009273/4596053‑4596185

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