BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,985,309	T→C	45.7%	L266P (CTG→CCG)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,985,309	0	T	C	45.7%	314.4 / inf	720	L266P (CTG→CCG)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand): ref base T (249/141); new base C (201/127); total (452/268)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 785 total aligned reads displayed.

TTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCG  >  CP009273/3985242‑3985375

  tGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgc                                                                  <  1:21991257/68‑1 (MQ=255)

   ggTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAccgct                                                                 <  1:12557349/68‑1 (MQ=255)

     tGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgct                                                                 >  1:30053152/1‑66 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACActgctgc                                                               <  1:25927494/67‑1 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAccgctgc                                                               <  1:18185136/67‑1 (MQ=255)

      gCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACAccgct                                                                 >  1:9349794/1‑65 (MQ=255)

           ccTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGaa                                                          >  1:16064388/1‑67 (MQ=255)

            gtttGGCATCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTAATGAAAACACCGCTGCTGGAAg                                                         <  1:31046998/66‑1 (MQ=255)

            cTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGc                                                        <  1:7114276/68‑1 (MQ=255)

            cTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAg                                                         <  1:19871480/67‑1 (MQ=255)

            cTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGc                                                        <  1:6996917/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                    <  1:24972455/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    <  1:3938285/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    <  1:4383860/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    <  1:30485108/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    <  1:26925654/68‑1 (MQ=255)

                ggCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    <  1:21785476/68‑1 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                    >  1:9127682/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                    >  1:23963444/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAt                                                    >  1:15863623/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    >  1:17438088/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    >  1:23112250/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    >  1:29376162/1‑67 (MQ=255)

                 gCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAt                                                    >  1:30816740/1‑67 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTAtt                                                   <  1:10291677/67‑1 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAtt                                                   <  1:25338189/67‑1 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAtt                                                   >  1:23348277/1‑67 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAtt                                                   <  1:7615827/67‑1 (MQ=255)

                  cAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTAtt                                                   <  1:6560606/67‑1 (MQ=255)

                      tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTcc                                                 <  1:26681872/65‑1 (MQ=255)

                      tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTgg                                              <  1:16386455/68‑1 (MQ=255)

                      tctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTgg                                              <  1:27797944/68‑1 (MQ=255)

                       ctACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTggg                                             <  1:2097698/68‑1 (MQ=255)

                           aaGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAAt                                          >  1:10829627/1‑67 (MQ=255)

                            tgagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAc                                        <  1:25702096/67‑1 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:31721438/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:11948772/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:14765585/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:11379529/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:15561685/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:31383290/1‑67 (MQ=255)

                            agagTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATa                                         >  1:4865274/1‑67 (MQ=255)

                              agTATCGATTCCCCATAGAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:22182914/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:13691375/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:18868121/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:15303104/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:16655697/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGACTATTCCTGGGAATAccc                                      >  1:14979354/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:13016730/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:24685022/1‑68 (MQ=255)

                              agTATCGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATAccc                                      >  1:16936977/1‑68 (MQ=255)

                                 aTCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGa                                    >  1:27562956/1‑67 (MQ=255)

                                 aTCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGa                                    <  1:26347570/67‑1 (MQ=255)

                                  tCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAc                                  <  1:17230629/68‑1 (MQ=255)

                                  tCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAAc                                  <  1:24283214/68‑1 (MQ=255)

                                   cGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGacc                                  >  1:10541452/1‑65 (MQ=255)

                                   cGATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAg                                                         <  1:5616211/44‑1 (MQ=255)

                                    gATTCCCCATACAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAcccc                                >  1:18480374/1‑68 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCa                               <  1:8300627/67‑1 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCa                               <  1:27987331/67‑1 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCa                               <  1:13268322/67‑1 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCa                               <  1:20477445/67‑1 (MQ=255)

                                      ttCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCAc                              <  1:21243008/68‑1 (MQ=255)

                                        ccccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACg                             >  1:8080478/1‑67 (MQ=255)

                                         cccATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACg                             >  1:13397949/1‑66 (MQ=255)

                                            aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgc                        <  1:8589824/68‑1 (MQ=255)

                                            aTACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTct                          >  1:28880080/1‑66 (MQ=255)

                                             tACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTct                          >  1:21602764/1‑65 (MQ=255)

                                              aCAAAGCGGTACTGAAACCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCCCGTctgctg                      >  1:22654444/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGaa                                           >  1:7897602/1‑47 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTACTGGGAATACCCGAACCCACGTctgctg                      >  1:20763482/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                      >  1:4161592/1‑68 (MQ=255)

                                              aCAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                      >  1:11675926/1‑68 (MQ=255)

                                               cAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                      >  1:4662186/1‑67 (MQ=255)

                                               cAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTctgctg                      >  1:8930670/1‑67 (MQ=255)

                                               cAAAGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTgg                     <  1:28333327/68‑1 (MQ=255)

                                                aaaGCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGc                    >  1:29099174/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                 >  1:28630444/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                 >  1:29695550/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                 >  1:10614921/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                 >  1:19168231/1‑68 (MQ=255)

                                                   gCGGTACTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGaa                 >  1:19378286/1‑68 (MQ=255)

                                                     ggTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAg               <  1:17892321/68‑1 (MQ=255)

                                                        aCTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat             <  1:12346562/67‑1 (MQ=255)

                                                        aCTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat             <  1:10033587/67‑1 (MQ=255)

                                                        aCTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat             <  1:2012820/67‑1 (MQ=255)

                                                        aCTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat             <  1:5187234/67‑1 (MQ=255)

                                                        aCTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat             <  1:11111228/67‑1 (MQ=255)

                                                        aCTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGat             <  1:14417333/67‑1 (MQ=255)

                                                         cTGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGatat           <  1:30307245/68‑1 (MQ=255)

                                                          tGAAAACACCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTATGCTGGCGAAAGATATc          >  1:23642160/1‑68 (MQ=255)

                                                             aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCaa        >  1:14134603/1‑67 (MQ=255)

                                                             aaaCACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCaa        >  1:27296553/1‑67 (MQ=255)

                                                               acacTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa     >  1:22361086/1‑68 (MQ=255)

                                                               acacCGCTGCTGGAAGCCTATTCCTGGGAATACCCGACCCCACGTCTGCTGGCGAAAGATATCAAACa     >  1:3324429/1‑68 (MQ=255)

                                                               acacCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa     >  1:19620997/1‑68 (MQ=255)

                                                               acacCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa     >  1:3615182/1‑68 (MQ=255)

                                                               acacCGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACa     >  1:26155853/1‑68 (MQ=255)

                                                                  ctgctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCg  <  1:6124389/68‑1 (MQ=255)

                                                                  ctgctgctgGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCg  <  1:9567805/68‑1 (MQ=255)

TTTGGTGCCAGCCTTTGGCAGCTCTACAAGAGTATCGATTCCCCATACAAAGCGGTACTGAAAACACTGCTGCTGGAAGCCTATTCCTGGGAATACCCGAACCCACGTCTGCTGGCGAAAGATATCAAACAGCG  >  CP009273/3985242‑3985375

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