BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,940,935	C→A	5.5%	intergenic (+3/‑31)	metW → / → metV	tRNA‑Met/tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,940,935	0	C	A	5.5%	635.7 / 2.3	201	intergenic (+3/‑31)	metW/metV	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand): ref base C (25/158); new base A (4/10); total (32/169)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

CCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCT  >  CP009273/2940871‑2940999

ccTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                                                >  1:18480982/1‑68 (MQ=35)

 cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                                                <  1:24089025/67‑1 (MQ=35)

 cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAt                                                               >  1:40767878/1‑68 (MQ=14)

 cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCCCCGGTTCAAATCCGGCCCCCGCAACCAATCaaa                                                                <  1:20682996/67‑1 (MQ=255)

      aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATAAAAT‑Ttgat                                                         >  1:19788021/1‑68 (MQ=255)

      aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCATAT‑Atgat                                                         >  1:23434398/1‑68 (MQ=255)

      aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCATAG‑Ttgat                                                         >  1:33106071/1‑68 (MQ=255)

      aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑Ttgat                                                         >  1:15414112/1‑68 (MQ=255)

      aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑Ttgat                                                         >  1:15465722/1‑68 (MQ=255)

      aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAT‑Ttgga                                                         >  1:40365907/1‑66 (MQ=255)

         tcgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGaa                                                      <  1:39243040/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGCAg                                                     <  1:29187038/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:22096050/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:19330935/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:19594633/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:20838379/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:32766344/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:38800622/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:23015350/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:23140591/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:3043784/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:25151535/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:26753905/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:13523043/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:9195414/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:8054739/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:42949839/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:42796516/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:14927056/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:40175330/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:15335919/68‑1 (MQ=255)

          cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:16775108/68‑1 (MQ=255)

          cggcgGGCTCTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAg                                                     <  1:42609167/65‑1 (MQ=255)

             cgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTa                                                   <  1:37111400/67‑1 (MQ=255)

                gCTCATAACCCTAAGGTCGTCGTTTCAAATCAGGCACCCGAAACTAATGAAAT‑TTGATGAATTaaaa                                                <  1:15768376/67‑1 (MQ=14)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:40115638/67‑1 (MQ=255)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:3403556/67‑1 (MQ=255)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:6071299/67‑1 (MQ=255)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:8252676/67‑1 (MQ=255)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:28369973/67‑1 (MQ=255)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:1051649/67‑1 (MQ=255)

                gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTaaaa                                                <  1:2394159/67‑1 (MQ=255)

                 cTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAT‑TTTATGAAGTAAAAGc                                              >  1:28760/1‑68 (MQ=255)

                 cTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAg                                               <  1:1223231/67‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGc                                              <  1:11617109/67‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGc                                              <  1:12553183/67‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGc                                              <  1:20576316/67‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCa                                             <  1:10152628/68‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCa                                             <  1:37982197/68‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCa                                             <  1:28466639/68‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCa                                             <  1:35968600/68‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCa                                             <  1:21593858/68‑1 (MQ=255)

                  tCATAACCCGAAGGTCGTCGGCTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGc                                              <  1:35686321/67‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:3974964/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:15976568/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:18993057/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:36522023/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:26407607/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:1509367/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:39436420/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:6956220/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:29098629/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:7886412/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:14273604/68‑1 (MQ=255)

                    aTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGt                                           <  1:31893157/68‑1 (MQ=255)

                      aaCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTAc                                         <  1:24014/68‑1 (MQ=255)

                       aCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACg                                        <  1:30926282/68‑1 (MQ=255)

                         ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGt                                      <  1:31587122/68‑1 (MQ=255)

                         ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGt                                      <  1:5685913/68‑1 (MQ=255)

                         ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGt                                      <  1:31175065/68‑1 (MQ=255)

                         ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGt                                      <  1:10619699/68‑1 (MQ=255)

                          cGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑ATCAAGCACTAAAATCAGTACgggg                                     >  1:35246305/1‑66 (MQ=14)

                             aGGTCGTCGGTTCAAATCCGGCCCCCGAAACCAATAAAAT‑TTTATGTATTAAAAGCAGTACGGTGAc                                   <  1:28707497/67‑1 (MQ=14)

                              ggTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGAcgc                                 >  1:11870890/1‑68 (MQ=255)

                                tcgtcgGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGAcgc                                 <  1:27036793/66‑1 (MQ=255)

                                 cgtcgGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGCTGACGCggg                              <  1:3423224/68‑1 (MQ=255)

                                 cgtcgGTTCAAATCCGGCCCCCGCAACCAATAAAAT‑TTTATGAAGTAAAAGCAGTACGGTGACGCggg                              <  1:8521486/68‑1 (MQ=14)

                                  gtcgGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATTAAGTAAAAGCAGTACGGTGACGCggg                              <  1:5065845/67‑1 (MQ=255)

                                  gtcgGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCggg                              <  1:1685751/67‑1 (MQ=255)

                                       ttCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGCAGTAAAAGCAGTACGGTGACGCGGGGTGGa                         <  1:35626815/67‑1 (MQ=255)

                                         gcaaTCCCGCCCGCGCCACCAATCAAAT‑ATGATTGAGTAAAAGCATTACGGTGACGCGGGGTGGagca                      <  1:11393432/66‑1 (MQ=17)

                                           taTCCTGCCCCCGCAACCAATAAAAA‑TTGAATCACTAAAAGCAGTACGGTGACGCGGGGTGGagcagc                    <  1:17870255/67‑1 (MQ=11)

                                           aaaccGGCCCCCGCAACCAATCCACT‑TTGGTGCACTACAACCAGTCCCCGGACGCGGGGTGGagcagc                    <  1:4220812/65‑1 (MQ=12)

                                           aaaccGGCCCCCGCAACCAAAAAAATATTGCAGAATTCACGGC‑CTCCCGTGACGCGGGGTGGagcagc                    <  1:13742907/65‑1 (MQ=2)

                                           aaTCCGGCCCCCGCAACCAATCAAAT‑TTTGTGAAGTAAACGCAGTACGGTGACGCGGGGTGGagcagc                    <  1:29853664/68‑1 (MQ=14)

                                           aaTCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGagcagc                    <  1:6634612/68‑1 (MQ=255)

                                           aaTCCGGCCCCCGCAACCAATCAAAT‑ATGGTCCACTATACGCAGTACGGTGACGCGGGGTGGagcagc                    <  1:29002859/68‑1 (MQ=17)

                                           aaTCCCTCCCCCGCAACCAATAATAT‑ATGATCCATTAAACGCATTACCGTGACGCGGGGTGGagcagc                    <  1:30808864/68‑1 (MQ=9)

                                             tCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTTCGGTGACGCGGGGTGGagcag                     >  1:21847778/1‑65 (MQ=255)

                                             tCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGagcag                     >  1:27130987/1‑65 (MQ=255)

                                                     ccGCAACCAAACAAAT‑TTCGGTAACTAAAAGCACTACGCTGACGCGGGGTGGAGCAGCCTGGTAGCtc          <  1:32044841/68‑1 (MQ=14)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:27595209/1‑68 (MQ=255)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:25476177/1‑68 (MQ=255)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:43083109/1‑68 (MQ=255)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:37677537/1‑68 (MQ=255)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:34392413/1‑68 (MQ=255)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:7834108/1‑68 (MQ=255)

                                                          aCCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCgg     >  1:29552407/1‑68 (MQ=255)

                                                           ccAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCggg    <  1:26169341/68‑1 (MQ=255)

                                                             aaTCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCt  <  1:33690076/68‑1 (MQ=255)

CCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAT‑TTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCT  >  CP009273/2940871‑2940999

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